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Summary. — Measurements of angular distributions of secondary par- 
ticles in intezactions of high energy primaries of cosmic radiation with the 
nuclei of photographic emulsion have been performed. A systematic 
increase in anisotropy of angular distribution in C.M. with increasing 
primary energy has been stated. It was found that in many high energy 
jets the number of secondaries in the vicinity of 7/2 C.M. angle in smaller 
than that predicted by the theories of Heisenberg or Landau. For ex- 
plaining these experimental results a model has been introduced in which 
the secondaries are emitted isotropically from two centres moving in 
opposite direction in the C.M. system. Some experimental arguments 
have been given for this isotropy. The relation between the velocity of 
two centres, the inelasticity coefficient, multiplicity and transversal mo- 
menta of secondaries has been discussed. 


1. — Introduction. 


TAKAGI (1) proposed a model of multiple meson production in nuclear col- 
lisions of very high energy based on the assumption that the generated particles 
are emitted independently from two centres moving in opposite directions 
in the C.M. system. 


(4) S. TAKAGI: Progr. Theor. Phys., 7, 123 (1952). 
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TAKAGI assumed that the secondary particles are emitted isotropically in 
the frames associated with each centre and as a result the particles are con- 
centrated in the forward and backward directions in the C.M. system. 

KRAUSHAAR and MARKS (2) suggested a similar model in which the centres 
are assumed to be highly excited nucleons. 

BHABHA (*) suggested a model based on the assumption of a spatial distri- 
bution of mass in the nucleons. This results in the emission of secondary par- 
ticles from distinct regions which move in opposite directions in the C.M. 
system. 

All authors mentioned above used the theories of FERMI (*) or HEISEN- 
BERG (5) for describing the phenomena of emission of secondaries from both 
centres. 

Koga (5) examined Takagi’s model in comparison with Heisenberg’s and 
Landau’s (7) theories. This comparison concerned mainly the problem of 
multiplicity of secondary particles in high energy jets and their composition. 

In a recently published paper (*) the authors used Takagi’s model for de- 
scribing the phenomena observed in high energy interactions. In that paper 
—which we shall further denote as «I »—-the authors analysed the angular 
distributions of secondary particles using the method of DULLER and WAL- 
KER (°) and LINDERN (!°). As a result of this analysis we observed an increase 
of anisotropy in angular distributions in the C.M. system with increasing energy 
of the primary. It was found that the model under consideration— two 
centres »—model as we called it in «I»—gave a consistent description of the 
observed experimental facts. For jets of the highest energies in which the 
narrow and the diffuse cone can be easily distinguished, the distribution ob- 
served in the cones taken separately was not in contradiction with an isotropic 
one in the centres own system. We found that the y i.e. the y-values of the 
emitting centres in C.M. system were much smaller than the corresponding 
values for nucleons. 

We related the factor y with the inelasticity coefficient K by means of 
the formula (1) given in «I» and formula (8) in the Sect. 6 of the present 


) 
) H. J. BHABHA: Proc. Roy. Soc. London, 219, 293 (1953). 
) E. FERMI: Phys. Rev., 81, 683 (1951). 

5) W. HEISENBERG: Zeits. Phys., 188, 65 (1952). 

) Z. KoBA: Progr. Theor. Phys., 17, 288 (1957). 

) 8. Z. BELEN'KIJ and L. Lanpav: Suppl. Nuovo Cimento, 8, 15 (1956). 

) P. Crox, T. CoGHEN, J. GIERULA, R. Hoznyxsxi, A. JuRAK, M. Mrgsowicz. 
IEWSKA, O. STANISZ and J. PERNEGR: Nuovo Cimento, 8, 166 (1958). 

N. M. DuLLER and W. D. WALKER: Phys. Rev., 93, 215 (1954). 

L. von LINDERN: Nuovo Cimento, 5, 491 (1957). 
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paper. The comparison of this formula with a very similar one given recently 
by EDWARDS et al. (11) will be given later. 

A similar attempt to describe some events of high energy collisions known 
from the literature was then given by G. Cocconi (2). 

The paper «I» was based on experimental material from our laboratories 
only. The present paper gives the results based on our extended material 
and on data from published papers (!*) and private communications (*). 


2. — Experimental procedure. 


From the whole available material we selected for analysis only inter- 
actions with primary particles: p, n, «, and with low excitation of the target 
nucleus (N, < 5). 

In order to obtain sufficient statistical significance in the energy evaluation 
from the angular distribution we restricted the events by the condition n, > 6. 

We also used information about angular distributions of neutral pions in 
the few cases in which an analysis of electron-photon cascades was performed. 

The events fullfilling these conditions were divided into 4 groups according 
to their primary energy (in eV) 


| TOUS, <= 10" 38 jets 374 tracks 
| CO a, 1058 42 jets 611 tracks 
| 102 < B,< 108 29 jets 475 tracks 
| LOS, 10 jets 143 tracks 


In separating the events into particular energy groups we used the for- 
mula (14) 


(1) Ig y. = lg cotg 0, 


(1) B. Epwarps, J. Losty, D. H. PERKINS, K. PINKAU and J. Reynoups: Phil. 
Mag., 3, 237 (1958). 

(12) G. CoccoNI: private communication. 

(13) M. ScHEIN, R. G. GLASSER and D. M. HASKIN: Nuovo Cimento, 2, 647 (1955); 
A. DEBENEDETTI, C. M. GareLLI, L. TALLONE and M. Vigone: Nuovo Cimento, 4, 
1142 (1956); M. Kosuipa and M. F. KAPLON: Phys. Rev., 97, 193 (1955); V. D. HOPPER, 
S. Biswas and J. F. DarBY: Phys. Rev., 84, 457 (1951); E. G. Boos, A. H. WINNICKIJ, 
7. S. TAKIBAEW and I. Ja. Casnrxow: Zurn. Exp. Theor. Fyz., 34, 622 (1958). 

(*) We are very much indebted to Dr. D. H. PERKINS of the H. H. Wills Labo- 
ratory, Bristol and to Dr. K. LANIUS of the Kern-physikalisches Institut der Deutschen 
Akademie der Wissenschaften zu Berlin for kindly sending us the angular distributions 


‘of their material. 


(4) C. CASTAGNOLI, G. CORTINI, D. Moreno, C. FRANZINETTI and A. MANFREDINI: 
Nuovo Cimento, 10, 1539 (1953). 
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(2) E, = 2y,—1, 
y, is therefore the Lorentz factor of a system for which.there is a symmetry 
in the forward and backward directions. In the case of nucleon-nucleon col- 
lision, y, is the energy of each nucleon in the C.M. system in rest mass units. 
It may be remarked that the analyzed material consists at least in the 
high energy groups of jets with relatively small numbers of shower particles n,. 
If we tentatively hold to the tunnel-model of nucleon-nucleus collision, and 
evaluate the average length of the tunnel expressed in the number of nucleons / 
in it, we get for the highest energy group in average /<1 and in the group 
102 eV+10"% eV J = 2 (see for instance (!)). 
The angular distribution was expressed in co-ordinates (°) 


F(0,) 
(8) x = log (y. tg 0,); ETES TEE 


where F(0,) is the fraction of tracks with an angle less than 9, with the jet 
axis. It follows from the transformation formula 


Trg peli el 
RER GE Cone 


where 9 and 0, are the angles in the C.M. and laboratory system respectively, 
that in case fp,» fh, & 1 the angular distributions in such co-ordinates give 
angular distributions in the C.M. system. In these co-ordinates the isotropic 
distribution in the C.M. system is represented as a line of slope 2 and the 
anisotropic distribution of HEISENBERG (!*) as a line of slope 1. The distri- 
bution given by the theory of LANDAU (7) is also an approximately straight 
line whose slope decreases slightly with increasing primary energy. The Landau 
distribution for the primary energy interval (101? :101%) eV is very similar to 
the Heisenberg distribution. 

LINDERN (1) used a = 1g (y, tg 0,) for plotting differential distributions, and 
showed that the typical angular distributions in the C.M. system e.g. iso- 
tropic, Heisenberg’s anisotropic, ete., could be well approximated by Gaussian 
curves, and the dispersion o can be taken as a measure of anisotropy. For 


instance, for isotropy we get o=0.36 and for Heisenberg’s anisotropic dis- 
tribution o = 0.7. 


(5) F. D. RoEsLER and C. B. A. McCusker: Nuovo Cimento, 10, 127 (1953). 


(5) K. SyMANZIK: in Kosmische Strahlung, 2nd ed. (Berlin-Géttingen-Heidelberg, 
1953). 
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3. — The dependence of anisotropy of angular distribution on the primary energy. 


In Figs. 1. 2. 3 and 4 we see the experimental differential angular distri- 
bution dN/dlg (y. tg0,) vs. lg (y.t20,), i.e. angular distributions in the C.M. 


= 
one tra 


ck 
=A 0 1 log(y_tg 6) 


Fig. 1. — Differential angular distribution 
for jets in the energy interval 101 eV to 
10%. eV. The continous curve represents 
the best fit of a Gaussian curve to the 
experimental data with o = 0.46. 


- 
one trach 
-1 0 


1 logly.tg@) 


Fig. 2. — Differential angular distribution 
for jets in the energy interval 1011 eV to 
102 eV. The continous curve represents 
the best fit of a Gaussian curve to the 
experimental data with o = 0.55. 


system. These figures correspond, to 4 energy intervals. As mentioned be- 
fore the o of these distributions may be considered as a measure of anisotropy. 
The values of o for different energy intervals are given in Table I. 


= 
one track 


#9 0 1 log(y tg CA 


Fig. 3. — Differential angular distribution 
for jets in the energy interval 101? eV to 
1013 eV. The continous curve represents 
the best fit of a Gaussian'eurve to the ex- 
perimental data with o = 0.66. 


= 
one track 


-1 0 1 


log(Y, tg 9, 


Fig. 4. — Differential angular distribution 

for jets with energy higher than 101 eV. 

The continous curve represents the best 

fit of a Gaussian curve to the experimental 
data with o = 0.96. 


We see from Table I the general increase of the anisotropy—expressed 
by o,,,—with increasing energy. The experimental material used in this in- 
vestigation is sufficiently large to give significant differences in o-values of 
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TABLE I. 
| oo Vinten oo | jou | 1012 | 1013 | 1014 
baci LE | | | 

FRS | 0.49 | 063 | 077 | 0.96 
Tins | 0.46 | 0.55 0.66 0.96 | | 
= | 1.25 | 1.50 2.0 | 4.0 | | 
“AG st 39%, | 8% | 12% | à 

| | 


neighbouring energy intervals. The levels of significance of these differences 
corresponding to the Kolmogorov-Smirnov test (17) are given in the last row 
of Table I and in Fig. 5. 


La 
10 aa 
10" ev- 10? ev CS 
10° ev- 10" ev more than 10 ev 
/ | 102 ev-10)3 ev 
05 if 


more than 10° eV 


10 2e - 107 ev 


0° eV- 10"! eV 


log(¥,t9 0) 
7 0 4 2 


Fig. 5. — Integral angular distributions for four energy intervals of jets. The maximum 
differences corresponding to the Kolmogorow-Smirnow test between neighbouring 
distributions are marked. 


(17) A. Kozmocorow: Giorn. dell’Ist. Ital. degli Attuari, 4, 1 (1933); H. SmIRNOW: 


Recueil Math., N.S. 6, 1 (1939); see also F. J. Massey: Amer. Statist. Assoc. Journ. 
(March 1951), p. 68. 
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The general dependence of anisotropy i.e. the o-value of the distribution, 
on the primary energy is roughly consistent with the prediction of Landau’s 
theory (*) as seen from the table. Besides we see that for lower energy inter- 
vals (1028-1041 and 10!!--10!) the experimental histograms are well approx- 
imated by Gaussian curves. However, for the energy interval (1072 --1018)-eV 
the agreement of the experimental histogram with the Gaussian curve of corres- 
ponding o is rather a poor one. There is a characteriste lack of tracks in the 
neighbourhood of the C.M. angle x/2. For energies higher than 10! eV the 
errors in the evaluation of y, and poor statistics make the comparison with 
any theoretical distribution difficult. 


4. — An attempt to analyze angular distributions on the basis of the «two- 
centre » model. 


In «I» the authors drew attention to the fact that in a number of high 
energy jets we observe a characteristic angular distribution which were shown 
for co-ordinates x, y (Fig. 7 in «I», Fig. 2 of Coccont’s paper (#2) and Fig. 6 
of the present paper). If these distributions significantly deviate from a straight 
line, as represented by the Heisenberg or Landau distributions, then we can 
say that the observed distribution is marked by a lack of tracks in the di- 
rections close to the angle 7/2 in the C.M. system. In the differential distri- 
bution dN/dæ vs. x this fact appears in the deviation from the Gaussian 
curve namely the dip in the vicinity æ = 0, which seems to be present in Fig. 3 
and 4. The authors attempted to discuss the characteristic angular distri- 
butions of the type shown in Fig. 6 on the basis of the two centre model. 

Let us assume (as in «I ») that after the collision of the nucleons the emis- 
sion takes place from two centres independently which move in opposite di- 
rections with the same velocity (in the simplest case) in the C.M. system and 
for which the corresponding value of y is y. Let us assume that in the system 
associated with the centre the emission of particles is isotropic and let us as- 
sume that the particles are emitted in this system with a uniform energy E_: 
These assumptions lead in a simple way to an angular distribution which, in 
the æ, y co-ordinates have the characteristic shape seen in Fig. 6. In this 
figure are represented several events with distinctly differing narrow and dif- 
fuse cones. Besides the events given in «I» we have included in this paper 
the additional events Nos. 18 and 26 from our material and events P5 and 
P20 published in the paper of EDWARDS et al. (11). 

For each event y, calculated from the formula lg y. = lg cotg 0, and », 
calculated from the individual cones by evaluations y, and y, for two emitting 
centres in the L system, in an analogous way to that introduced in the paper 
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of CIOK et al. (18) are given in the figure. We have (see Appendix) 


(4) ve = ie + VAI VA—14-1) & pin. 
From the values y, and y, we also calculate the value of 7 from the formula 


(5) = pra 
2y, 


In addition, the values of y (denoted as 7, in Fig. 6) calculated from the co- 
efficients of inelasticity (formula (8)) are given for the events for which these 
coefficients were directly measured. The solid lines in Fig. 6 represent distri- 
butions calculated on the basis of the «two-centre model» for values of y 
and y given in the figure. 


5. — Angular distribution of secondary particles in the system of the emitting 
centre. 


In the preceding Section we gave a comparison of the observed angular 
distribution for several events with the distribution calculated under the as- 
sumption of the « two-centre » model. Here we assume that emission in the 
system of the given centre is isotropic. There exists the possibility of checking 
this assumption experimentally. In the events where the cones are distinctly 
separated we can investigate the distributions in the individual cones treated 
as separate jets. 

In «I» an analysis was made in this direction. Calculations of y, and > 
from the formula 


(6) log ya = log cotg 0, , 


were made for individual forward and backward cones of our jets Nos. 155, 115 
and for the star «S». Then the total integral angular distribution for this 
set of forward and backward « jets » was computed. Fig. 6 in «I» shows that 
there is no contradiction with the assumption of isotropy of the distribution 
if it is assumed that the particles in the individual cones are emitted by two 
centres moving with the Lorentz factors y, and y, in the L system. The results 
of such analysis are shown in Fig. 7 on the basis of more extensive material. 


(8) P. Crox, M. Danysz, J. GIERULA, A. Jurax, M. Mrgsowicz, J. PERNEGR, 
J. VRANA and W. WoLTER: Nuovo Cimento, 6, 1409 (1957). 


8711 


750 P. CIOK, T. COGHEN, J. GIERULA, R. 


Pcs VITO 


IIOLYNSKI, A. JURAK, M. MIESOWICZ, ETC. 


Represented here are the integral angular distributions for narrow and dif- 
fuse cones separately. We believe that the results presented here confirm se 
opinion expressed in «I» that both distributions can be regarded as isotropic. 


12 


Narrow cones Diffuse cones 


08 0.8 


04 04 


tog(y .tg@,) 


log(y,199,) 


04 08 04 08 
-04 
a) Pa L) 
12 
1.6 
Fig. 7. — Integral angular distributions for narrow (a) and diffuse (b) cones treated 


separately as independent jets. Here are presented no. 115 (18), « S »-star (18), «5, P20, 
P2054, P10 and P2484 (11) and private communication from Dr. PERKINS. 


6. — The problem of transversal momentum in the « two-centre » model. 


Now we would like to consider the consequences, aside from the shape of 
the angular distribution, which result from the adoption of the « two-centre » 
model. In the work of the Bristol group (1) there was introduced a new de- 
finition of the coefficient of inelasticity Æ' from which the following relation 
results : 


(cosec 6), = 


— 
EN | 
— 


4 (4) 0 angle in C.M.; 
3 d p, transversal momentum. 


As we have shown in «I» ((Eq. (1)) an analogous relation results from our 
model as an expression of the principle of the conservation of energy: the 
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total energy of the centres moving in the C.M. System with a y-value equal 
to y and emitting secondary particles with a mean energy E, (total energy 
of a particle in the rest system of the centre) is equal to the energy available 
in the C.M. multiplied by the factor K (K, coefficient of inelasticity) : 


(8) 1.50, y = 2Ky,. 


In our approach the measure of anisotropy introduced in the Bristol paper 
—(cosec 0),,—finds physical expression in y. The mean y-values for each 
energy interval were determined and given in Table I (see also Appendix). 

If we tentatively adhere to this model, we can, from the above relation 
determine: 


a) the value E, assuming that the upper limit of the quantity 
nyly. (= 2K/1.5E,) for the set of all observed jets corresponds to the highest 
possible value of the inelasticity coefficient, i.e. K —1 (say HL, = (H,),_,); 


b) assuming this value of (H,),_, as the energy of the emitted particles 
also for collisions with AK < 1 we can establish the relation between the quan- 
tities n,7/y. and K, namely: 


(9) de His (12) = 


We can check this relation for jets in which the inelasticity coefficient was 
determined by direct measurements (see the numerical data in Fig. 6). 

Let us make this analysis, in effect, analogous to that of the Bristol group, 
for the entire experimental material at our disposal. The results are shown 
in Fig. 8. The quantity n,7/y. for jets in the energy range from 10! eV to 
1014 eV are shown as a function of y,. The quantity y was determined from 
the angular distributions on the basis of the «two-centre » model (see Ap- 
pendix). For energies below 10% eV only the mean for the two groups of jets 
in the energy intervals 10% eV = 101eV and 10! eV-- 10! eV are given. 

It should be stressed that value of (£,),_1 = 0.4 GeV estimated in this 
way is in good agreement with the measurements of the transversal mo- 
menta (1:11). 

As may be seen from Fig. 8 the material investigated here characterized 
by a small excitation of the target nucleus (N,<5) confirms the conclusion given 


(19) O. Minakawa, Y. Nisuimura, M. TZUZUuKI, H. YAMANOUCHI, H. Aizu, H. HASE- 
GAWA, Y. Isuu, S. TOKUNAGA, Y. Fusimoro, S. HASEGAWA, J. NIisHIMURA, K. NIU, 
K. NisuixawA, K. IMAEDA and M. Kazuno: I. N. S. J.-7, preprint (Tokyo, March 1958). 
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originally by the Bristol group (11) about the decrease of the mean pis of K with 
increasing energy. From the smallness of K for very high energies and from 
the fact that 7 is much smaller than y, it follows that it would be difficult to 
identify the centres of emission with two «excited nucleons ». 


Ps scale nsŸ -value accepted n 


Us Y as. corresponding to K=1 scale 


0.5 


02 


5 10 20 50 100 7c 


Fig. 8. — Analysis of the formula 1.5n,E,y = 2Ky,, and the energy dependence of Æ 

based on the « two-centre » model. © correspond to jets with energy higher than 107? eV. 

© correspond to groups of jet energy smaller than 10!2eV. Numbers in the vicinity 
of some points denote experimental K-values. Other details: see text. 


7. — Conclusions. 


1) In interactions of nucleons of cosmic radiation with the nuclei of 
photographic emulsion (with N, <5) in the limits of energy 10! eV = 1014 eV 
the angular distributions of secondary particles in the C.M. approximated by 
Gaussian curves show an increase in anisotropy with energy in agreement with 
the predictions of the Landau theory. 


2) It was found, however, that in many high energy jets the angular 
distribution quite markedly departs from the angular distribution predicted 
by the theories of Landau and Heisenberg in their present form. This is also 


reflected in the total differential distribution by a dip in the vicinity of the 
angle 7/2 in the C.M. 


3) The «two-centre » model that has been introduced by TAKAGI seems 
to describe the shape of the angular distribution in the C.M. with better agree- 
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ment with experiment. In this model it is assumed that the secondary par- 
ticles are emitted isotropically with the same energy from two centres moving 
in the C.M. with y considerably less than the y of the nucleons. 


4) The « two-centre » model permits the factor 7 of the centres in C.M. 
to be related to the coefficient of inelasticity, multiplicity, energy of secondary 
particles in the system of the centre and the energy of the primary nucleon, 
as (cosec 0),, was in the paper of Bristol group. From this relation and from 
the observed angular distributions the values of the transversal momentum 
obtained are in agreement with experiment. 


* OK OK 


The authors wish to thank Dr. A. ZAWADZKI and Dr. P. ZIELINSKI for their 
comments on the statistical analysis of the experimental material. 


APPENDIX 


Formulae (4) and (5) can be derived in the following manner: Let y,, y: 
be the Lorentz factors of the two emitting centres in L. If we assume the 
angles 0, and 6,,, formed with the shower axis by the two centres respectively 
64, = 9. = 9 and the Lorentz factors y of the two centres in the C.M. equal, 
than the usual transformation formulae from L. to C.M. read 


(10) 


or inversely 


1 = pe + VyE 1 V1, 


(11) ne 
Va = Vy, — VyÈ—1 y? — 1. 


Comparing these equations and eliminating y we immediately receive for the 
velocity of C.M. in L 


vi — V2 _Vyi=1+ V1 
nn vai nth? 


and for y = 1/(V 1— f?) after some simple calculations, 


(13) ve = tp + Vyi 1 Vyi 1 +1) À pp 
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The right hand side of (12) can be directly obtained from energy and momentum 
conservation in L. 7 can be obtained either directly from (10) or by adding 
equations (11) 


pa V1 + Ya 
2ye 


(14) 7 


The value y can be also evaluated if we know 


(15) di V > [log vo TA y 2 (e; = log ye)? 


If we assume isotropic emission from two centres, as in the model under con- 
sideration we get two groups of values, each having o = 0.36 with respect 
to its mean value: — log y, or — log y, respectively. Let us denote the dis- 
tance between these two values by 26. From equations (11) follows: 


7 + Vyi—1 V1 
(16) 20 = og 4 = log Ye? ARE SS e: rw 2 log 27, 
Va yy -Vy-1V7?-1 
(for large y. and not too small y). 
For o with respect to — log y. = — 4(log y, + log y») we obtain immediately: 
(17) o = V0.36 + d? . 


This relation enables us to evaluate y in a very simple manner even if the 
narrow and diffuse cones can not be easily distinguished. 


RIASSUNTO (*) 


Si sono eseguite misure di distribuzioni angolari di particelle secondarie nelle inte- 
razioni di primari della radiazione cosmica coi nuclei dell’emulsione fotografica. Si è 
constatato un aumento sistematico dell’anisotropia della distribuzione nel sistema del 
centro di massa coll’aumento dell’energia primaria. Si è trovato che in molti jets di 
alta energia il numero dei secondari in prossimità dell’angolo 7/2 nel sistema del centro 
di massa è minore di quello previsto dalle teorie di Heisenberg e di Landau. Per spie- 
gare tali risultati sperimentali, è stato introdotto un modello nel quale i secondari 
sono emessi isotropicamente da due centri che si muovono in direzioni opposte nel 
sistema del centro di massa. Si espongono alcuni argomenti sperimentali a sostegno 
di tale isotropia. Si è discussa la relazione tra la velocità dei due centri, il coefficiente 
di inelasticità, la molteplicità e i momenti trasversali dei secondari. 


(*) Traduzione a cura della Redazione. 
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Possible Consequences of the Isobar and Fire Ball Models. 


Z. KOBA 
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Summary. — Two possible consequences, characteristic of the isobar 
model and the fire ball model of the multiple meson production, are 
discussed. One is the azimuthal asymmetry of secondary particles, and 
the other is the apparent anomaly of the inelasticity in meson-nucleon 
collision. 


1. — Introduction. 


Recently a model for multiple particle production at extremely high energy 
has been proposed by Nrv (!), the Polish group (2), and Cocconi (*). Ac- 
cording to them, the collision of two energetic nucleons gives rise to two «fire 


(1) K. Niu: Utyusen-Kenkyu (Mimeographed circular in Japanese for cosmic ray 
studies), 3, 85 (1958); K. Nru: Lecture at the semi-annual meeting of the Physical 
Society of Japan (May 1958). 

(2) P. Crox, T. COGHEN, J. GreruLa, R. HozyXski, A. JURAK, M. MIRSOWICZ, 
T. SANIEWSKA, 0. STANISZ and J. PERNEGR: preprint. We are grateful to these authors 
for sending us their work before publication. 

(3) G. Cocconi: preprint. We are grateful to Prof. S. Hayakawa for showing us 


this work. 
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balls », each of which separates from the incident nucleon, flies out, and after 
a very short time disintegrates into a number of mesons. 

This «fire ball model» (*) shares with the isobar model (+), which was 
presented by one of us (S. T.) (*) and later worked out by KRAUSHAAR and 
Marks (5), the assumption of a hypothetical metastable state in the inter- 
mediate stage. In this note we should like to point out two possible conse- 
quences of such models, which might be detected experimentally. 


2. — Azimuthal asymmetry of secondary particles. 


In these models it is conceivable that the azimuthal distribution of se- 
condary particles can sometimes show more asymmetry than is expected from 
simple statistical fluctuation. 

Suppose that partial waves of very high angular momenta are effective in 
the collision process due to certain properties of interaction, and a larger part 
of this angular momentum is assigned to the 
intrinsic angular momentum (or spin, so to 
speak) of isobars or fire balls. 

Let the isobar or the fire ball have linear 


pei momentum, p*, whose direction can be slightly 
Sd a different from that of the incident particle, p. 
cala Thus the spin of the isobar or fire ball will be 

P 


quantized in the direction of p x p*; this implies 
that a particular direction with respect to the 
azimuthal angle does exist. That is to say, the 
emission in the plane defined by p and p* will 
be most favoured when the perpendicular com- 
ponent of the spin of the isobar or the fire ball happens to be large (*) (Fig. 1). 

Detailed calculation of possible azimuthal asymmetry has already been 
carried out by KRAUSHAAR and MARKS (5). The point we should like to 
emphasize here is the essential role in the above argument played by the inter- 
mediate stage. One might be tempted to imagine the nucleon-nucleon col- 
lision as two classical balls, passing by each other at a certain distance o (Fig. 2). 


Target diagram 


Fig. 1. — Emission with 
azimuthal asymmetry. 


) Abbreviated hereafter as FB-model. 

) Abbreviated hereafter as I-model. 

4) S. TAKAGI: Prog. Theor. Phys., 7, 123 (1952). 

5) W. L. KRAUSHAAR and L. J. MARKS: Phys. Rev., 83, 326 (1954). 
) 


(*) One can imagine emission from a spinning ball: most of the emitted particles 
will be found in the plane of equator. 
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Then the system would have an angular momentum of in the direction per- 
pendicular to the line defined by the two parallel tracks of incident balls, 
so that the system would again have a preferential 
direction with respect to the azimuthal angle, quite 
irrespective of the I model or FB-model. But such 
a reasoning is not correct, because here the quantum- P 
mechanical system of two colliding nucleons is de- 
scribed by a plane wave (*) which has an axial 
symmetry around the direction of incidence; in 
other words, the system is represented by an equal- 
weight superposition of all the states which are obtained by rotating the 
Fig. 2 around the vector p. 

Thus it is obvious that the actual existence of a metastable state, which 
through the interaction with external fields defines the momentum p* and 
consequently the plane px p*, is a necessary condition for the detection of 
azimuthal asymmetry. The situation is analogous to the polarization exper- 
iment by a double scattering. 

In this way we conclude that, if we find more azimuthal asymmetry than 
statistical fluctuations, then it would make the models with intermediate meta- 
stable states rather plausible, while the reverse argument does not hold, of 
course. 


pie 


Fig. 2. — Classical picture 
of collision. 


3. — Apparent inelasticity of the meson-nucleon collision. 


Many of the secondary jet-showers are regarded as induced by mesons. In 
such cases there is no reason to believe that the collision will be symmetric. 
But experimentally most of the secondary showers exhibit symmetric angular 
distribution of emitted particles in suitably chosen reference systems. In the 
following we shall show that, according to the I- or FB-model, an asymmetric 
meson-nucleon collision can give rise to apparently symmetric distribution of 
secondary particles, and moreover, this circumstance can lead to an erroneous 
(sometimes unacceptable) estimation of the inelasticity. 

It is plausible to imagine that a meson-nucleon collision will result in emis- 
sion of not two, but one isobar or fire ball, which then disintegrates. In the 
I-model the isobar will stop in the center of mass system (+), if the incident 
meson loses its whole energy by the collision; but, in general, when the incident 
meson loses only a fraction of its original energy and goes on proceeding in 


an 


(*) In the actual condition the incident beam is located by the silver grain, which 
is infinitely large compared with the target, so that it can be described by a plane wave. 
(+) Abbreviated hereafter as c.m.s. 
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the initial direction, the resulting isobar will move in the opposite direction 
to the incident meson in the c.m.s. In the FB-model we can imagine, as far as 
kinematical relations are concerned, that the above mentioned isobar sepa- 
rates further into a nucleon and a fire ball, the nucleon proceeding faster than 
the fire ball in the opposite direction to the incident meson. Thus the fire ball 
will be moving slower (*) in the c.m.s. than the isobar, when the inelasticity 
has a given value (cf. Fig. 3). 


I-model FB-model 
N 
Before collision ca e Rife A 
After collision + co <e- 4 =e 
I H N FB H 


Fig. 3. — Partially inelastic collision. 


Let us denote the Lorentz factor of the laboratory system (+) with respect 
to the c.m:s. by y:, and that of the isobar or the fire ball in the c.m.s. by 7 


Tv: 
Then the Lorentz factor, y,,, of the isobar or the fire ball in the Ls. is given by 


(1) Van = YoVe TV — Vy —al ’ 


where the minus and the plus signs correspond respectively to the directions 
of motion of the isobar or the fire ball in the c.m.s. being the opposite and the 
same as the incident meson. 

Now it is obvious that, when we estimate the incident energy by the usual 
analysis of the angular distribution of secondary particles, the formula: 


1% 
(2) — logy =— > logtg@,, 


Vs i=1 
ns: number of shower particles, 0;: angle of emission of the i-th particle in 
c.m.s., will not give y=y,, but y=}, since the secondary particles are 
symmetrically emitted in the rest system of the isobar or the fire ball. (Here 
we neglect the effect of possible inclusion of the incident meson and the 
target nucleon in the summation of (2), since N,>>1.) 
Thus the estimation of the incident energy in the 1.s. through the relation 


(3) B~2M,, 


(gi In algebraic sense. Sometimes the direction of motion in c.m.s. can be reversed 
(i.e. if can move in the direction of the incident meson). 
(+) Abbreviated hereafter as Ls. 
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M being the mass of the target, will lead to an apparent value 
(4) E, = 2My? 


Consequently we shall obtain an apparent inelasticity: 


(5) E (2), 


where 7 represents the «true » inelasticity and is defined by the ratio of energy 
loss to the total energy of the incident meson in the 1.s. (*). 
In the Appendix I and II it is proved that 


(6) Map = co , 


where n’ stands for the rate of energy loss of the incident nucleon in the c.m.s. 
When the nucleon moves with the excited meson cloud (I-model), 7’ is put 
equal to 1. 

Thus in the I-model the apparent inelasticity always turns out one, even 
if the incident meson loses only a small fraction of its energy. If by chance 
the energy of the incident meson, which has thus retained a large energy, is 
measured and taken into the sum of the energy of the secondary particles 
then it will lead to an unacceptable result: Ya, >1. The situation in the 
FB-model model is more or less similar. 

It is suggested, therefore, that we should be rather careful in the analysis 
of secondary jets and examine whether the above-mentioned effect really 


exists. 
kK *K 


In conclusion, we should like to express our gratitude to the Kyoto, Osaka 
and Tokyo members of the research group of extremely high energy pheno- 
mena for their discussions. 


APPENDIX I 


Let the incident meson and nucleon lose fractions 7 and »' of their energies 
in the c.m.s. respectively. The energy, E*, and momentum, p*, of the fire 
ball are determined by the energy and momentum conservation : 


(A.1) 2My,=(1—n)My:+(1-n')My.+E", 
(A.2) 0=(1—n)My;— (1—7') My. — v* - 


(*) The inelasticity in the c.m.s. is (7 + n')/2. 
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These relations are valid when (1—n)y.> 1, (1 —#7')y. > 1, but they hold also 
when 7 = lor 7'— 1, provided that the nucleon moving with the cloud (isobar) 
is regarded as if it had disappeared and, instead, #* and p* are identified with 
energy and momentum of the isobar (*). We shall use, therefore, the above 
relations for all values of 0 <7 <1, 0<7'<1. ' 

The Lorentz factor of the isobar or the fire ball in the c.m.s. is given by 


E* 
AD (E*? con p*} ’ 


(A.3) 


or, with the help of (A.1) and (A.2), 


a 2) 


Now, by (1) 


(A.5) Yor =VV:FVR_1V_-1ar{nFVy=-1} ify >1, 


(+ corresponding to n < 7’). 
We have from (A.4) 


(A.6) Vyi—T= 


— 
=" 

~ 

~ 


bol 
= Sd 
Is] 3/3] 

| 

È ast 
AL 

eb 

IN 


to} 


3 


so that, from (A.5) and (A.6) 


(A.7) PESATA 
Va Va 0-3 
whence 
(A 8) Hg EN (2) n 
Var 


APPENDIX II 


The relation (6) can be easily understood by i ini i 

TI i y imagining a hypothetical 
collision of a hypothetical meson with the nucleon, which Richie, es rise 
to a fire ball identical with that of the real collision, but in which both particles 
lose the same fraction 7’ of the incident energy (Fig. 4) 


*) In other word : 
oe words (A.1) and (A.2) hold, as far as we can neglect the effect of rest 
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Let the hypothetical meson have an incident energy (7/7')My. (instead 
of My, of the real meson) in the original c.m.s. of the real collision. In this 
reference system the hypothetical meson and the nucleon lose the energy 7 My, 
and y'My, respectively, giving rise to a fire ball identical with the real one. 


Real c.m.s. Hypothetical c.m.s. 
2 My, My Ram 220 
n n N 
Before collision —o> «e da Rene 
After collision Pra +) > e () o> 
My ny)? My, VE am 1 (1-y)M 
My(1—7') th Ve À 7')My, 7 | n') My. 


Fig. 4. — Hypothetical collision. 


The rate of energy loss of the incident particles is n' for both of them, and 
as far as we can neglect the rest mass of nucleon and meson, this rate of 
energy loss is invariant with respect to the Lorentz transformation. Therefore 
in the c.m.s. of the hypothetical meson and the nucleon, the hypothetical colli- . 
sion will take place in a completely symmetric way, and the fire ball will be 
at rest and decay symmetrically. 

We can further show, that the Lorentz factor of this hypothetical c.m.s. 
with respect to the L.s. is identical with the ya, given by (A.7). For this purpose 
let us first prove that the product of the energies of two incident particles 
is independent of the Lorentz frame as far as we can neglect the effect of rest 
mass. Due to a Lorentz transformation the energy of one particle is multiplied 
by a factor J(1 +), and that of the other by a factor [(1—v), where 
JT = (1—v)t is the Lorentz factor of this transformation. Thus the product 
of both energies is multiplied by ri—v)(1 + v) =1, that is, remains 
invariant. Consequently the incident energy € of each particle (hypothetical 
meson and nucleon) in the hypothetical center of mass system is given by 
the geometrical mean of the incident energies in any system of reference: 


= Mye My: 


or 


(A.9) & = V2 My. . 


Evidently this implies that the Lorentz factor of the hypothetical c.m.s. 
with respect to the 1.8. is expressed as V n'/n:ye, Which is equal to ya, in (A.7). 
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Thus it is seen that, by analyzing the decay product of the fire ball with 
the assumption of a symmetric collision, we are led to the process of the 
hypothetical collision with the inelasticity 7’, and not to the real collision with 
the inelasticity n (or ((7 + 7')/2 in the c.m.s.). 


RIASSUNTO (*) 


Si discutono due possibili conseguenze, caratteristiche del modello isobarico e del 
modello a «fire ball » della produzione multipla dei mesoni. Una è l’asimmetria azimu- 
tale delle particelle secondarie, l’altra è l'apparente anomalia dell’anelasticità della 
collisione mesone-nucleone. 


(*) Traduzione a cura della Redazione, 
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Summary. — A homogenous sample of 540 7-meson decays found in 
emulsion stacks exposed to the Berkeley Bevatron has been studied. 
Most of the events have been located by track scanning, and the effects of 
geometrical and scanning biases have been minimized. The results are 
discussed in detail and then combined with those of a comparably large 
sample examined by BALDO-CEOLIN et al. (2). The most probable spin- 
parity assignment is 0° and the deviations of the positive and negative 
pion energy spectra from those predicted by the relativistic statistical 
factor, already noted by BALDO-CEOLIN et al., have been confirmed and 
estimated. 


(*) A preliminary report on this study was presented at the Padua- Venice Con- 
ference on Mesons and Recently Discovered Particles (1957). A: 

(+) Now at Dunsink Observatory, Institute for Advanced Studies, Dublin. 

(x) Now at the College of Science, Pilani, India. 
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1. — Introduction. 


The energy and angular distributions of the three pions emitted in the 
decay of the +-meson have been calculated by several authors (*) under various 
assumptions of spin and parity of the parent particle and of the angular mo- 
mentum configuration in the decay. The approximations and assumptions in 
the derivation of these curves are discussed fully in these references. The ex- 
perimental results agree best with the curves calculated for the spin-parity 
assignment (07). 

BALDO-CEOLIN et al. (2) have given a useful bibliography of the earlier work, 
which we shall refer to as A. 

The largest homogeneous sample of t-meson decays studied hitherto has 
been that of 419 events reported in A. In this work, attention was drawn 
to certain apparent discrepancies between the experimental results and the 
calculated curves. 

In the present work a second large homogeneous sample has been studied, 
with the aim, as far as possible, of avoiding selection bias. 

The sample of t-decay studied was comparable in size with that studied 
in A. It is a feature of the present experiment however, that almost all the 
events have been found by track-scanning, or modified track-scanning, in 
contrast to that of BALDO-CEOLIN et al., in which most of the events were 
found by area scanning. 

Of the events studied in the present work, 87 reported already by N. N. 
Biswas et al. (*) have been re-examined in the light of the criteria laid down 
in Sect. 2°1 and 2°2. 

The selection procedure, measurements and results are described in Sect. 2°1, 
2°2, and 2°3 respectively. 

In Sect. 3 the results are discussed separately, and combined with the 
results of the 419 events of A. 


() R. H. Datirz: Proc. Phys. Soc., 66, 710 (1953); Phil. Mag., 44, 1068 (1953); 
Phys. Rev., 94, 1046 (1954); 99, 915 (1955); Proc. Phys. Soc., 69, 527 (1956); E. FABRI: 
Nuovo Cimento, 11, 479 (1954); Suppl. Nuovo Cimento, 12, 205 (1954); G. Costa and 
L. TAFFARA: Nuovo Cimento, 3, 169, (1956). Y. EISENBERG, E. Lomon and $. Rosen- 
DORFF: Nuovo Cimento, 4, 610 (1956); I. S. SHAPIRO, E. I. DoLINSKY and A. P. MisHa- 
KOVA: Nucl. Phys., 3, 60 (1957). 

: (?) M. BaLDO-CEOLIN, A. BONETTI, W. D. B. GREENING, S. LIMENTANI, M. MERLIN 
and G. VANDERHAEGHE: Nuovo Cimento, 6, 84 (1957). 

(8) N. N. Biswas, L. CECCARELLI-FABBRICHESI, M. CECCARELLI, K. GOTTSTEIN, 

N. C. VARSHNEYA and P. WALOSCHEK: Nuovo Cimento, 3, 825 (1956). | 
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2. — Experimental procedure and results. 


21. Method of scanning. — All the events were found in 3 stacks of G5 
nuclear emulsion strips 600 ym thick exposed to K*-meson beams at the Be- 
vatron, the K; stack, the K; stack of the European Collaboration, and a 
third smaller stack exposed on behalf of W. F. FRy. 

289 events were found by track scanning in the course of experiments on 
the interactions of K* particles (4), 221 events by modified track scanning as 
described in (5), and 30 events by area scanning. 

Due attention has been paid to minimizing the effects of bias discussed 
in detail in Sect. 211 of A, and the geometrical criteria proposed by BRENE 
et al. (*) have been strictly satisfied. 

A decay configuration wicht might escape observation is that in which the 
negative x-meson is emitted in the direction of the primary t-meson with low 
energy ~ 1 or 2 MeV, and results in a zero prong capture star. 

In such a case, the change in grain density at the point of decay is liable 
to elude detection and the tracks of the 2 positive m-mesons, being fast and 
almost collinear, might be mistaken for a single crossing track. 

An approximate calculation predicts an expectation of one such event in 
the present sample. 


2°2. Acceptance criteria and analysis of events. — Only those t-meson decay 
events were selected for study of which the origins lay at distances greater 
than 8 mm from the borders of the stacks. In this case, as shown by BRENE 
et al. (8), at least two secondary pions should be brought to rest within the 
stack, making it possible to deduce unambiguously the signs of the charge of 
the three decay products. 39 events failed to satisfy this criterion and were 
excluded. For 499 events, the angles included between the directions of the 
emitted secondary particles were measured, and such further observations 
were made that the sign of the charge could be deduced for each of the three 
r-mesons. A stereographic projection chart was used, in the usual manner to 
reduce the projected angle and dip measurements to space angles. By use 


(4) N. N. Biswas, L. CECCARELLI-FABBRICHESI, M. CECCARELLI, K. GOTTSTEIN, 
N. C. VARSHNEYA and P. WaLoscHEK: Nuovo Cimento, 5, 123 (1957); M. BALDO- 
Coin, M. CRESTI, N. DALLAPORTA, M. Gricu, L. GUERRIERO, M. MERLIN, G. A. SA- 
LANDIN and G. ZAGO: Nuovo Cimento, 5, 402 (1957); B. BHOWMIK, D. Evans, S. NILSSON, 
D. J. PROWSE, F. ANDERSON, D. Krrrs, A. Kernan and J. LOSTY: Nuovo Cimento, 


6, 440 (1957). 
(5) G. ALEXANDER, R. H. W. Jounsron and C. O’CEALLAIGH: Nuovo Cimento, 


6, 478 (1957). : 
(6) N. BRENE, K. H. Hansen, J. E. Hooper and M. ScHARFF: Nuovo Cimento, 


4, 1059 (1956). 
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of the chart, the coplanarity of the three outgoing 7-mesons was checked for 
each acceptable event. In no case, could a departure from coplanarity be 
established which was significant, taking into account the uncertainty (typi- 
cally — 1°) in the individual measurements of the space angles. Range measu- 
rements were made on the secondary particles from 318 7-decays. In 258 
events, the energy of the 7--meson was found by direct range measurement, 
and in 60 of these, the range of one of the 7+-mesons was also measured. 
The energies of the two 7*-mesons were found by range measurement for 
60 events. 
For clarity this information is summarized in Table I. 


TABLE I 
No. of events Te Teli Tu | Angles 
| ; isti + | | 
| 1 Range Range Range | yes | 
| 60 No Range Range Range yes 
60 | Range Range No Range yes | 
197 Range No Range | No Range | yes | 
181 No Range No Range | No Range yes | 


499 Fully analysed | 
39 Excluded for reasons of geometry 
2 Could not be fully analysed, see text 


| 
: ha Net. | 
540 Total | 

Apart from the 499 events fully analysed, 2 additional events were investi- 
gated which satisfied the selection criteria, but for which the energy of the 
n particle could not be determined. In each of these events, one particle 
escaped from the stack and one particle was positively charged. In one event 
the third particle interacted, and in the other came to rest between pellicles. 
See Appendix II for a partial analysis of these events. The ranges were care- 
fully measured and the uncertainties in the measurement of length were esti- 
mated to lead to an uncertainty of — 0.7% in energy. This was neglected in 
comparison with the error in the energy estimate due to straggling (AH/H~2°%). 
When it was found impossible, or experimentally inconvenient, to measure 
directly the momentum of the x--meson, owing to particles crossing cut edges 
or going into parts of the stacks not to hand, the ranges of the two positively 
charged secondaries were measured. Their energies having been deduced, the 
energy of the m~-meson was derived either by subtraction from the Q-value 
of the 7-meson decay 74.9 MeV, or by closing the momentum triangle. In 
181 events the energies of the secondary particles were obtained from the 


ie) 
N 
n 
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momentum triangle defined by the space angles at decay, normalizing the 
energies to the Q-value. 

As a check, the energies of the 258 7 -mesons whose ranges were directly 
measured were also computed from the angles in this way. These energy esti- 
mates were in good agreement with those obtained from the range. 

The range energy curves published by BARKAS et al, (7) were used in all 
cases to deduce the energies of the mesons from their ranges. No decay in 
flight of a 7-meson was observed. 


2°3. Experimental results. — Fig. 1 shows the observed energy distributions 
of the negative and positive z-meson secondary to the 499 7-meson decays, 


32 


Tt” Energy spectrum 
499 Events 


28 


24 


20 


16 


lumber of events 


12 


n* Energy spectrum 
gop nt 


Number of events 


Fig. 1. 


(?) W. H. Barxas: U.C.R.L. Reports no. 3768, 3769. 
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plotted in class intervals of 2 MeV. Uncertainties in the energy estimates 
are greater for the higher than for the lower regions of the spectra since, in 
most cases, these energies were not estimated directly from the range of a 
single track, but depended rather on measurements of angles, or upon the 
measured ranges of two tracks of lower energy. Similarly, a large error is 
associated, with the calculated energy of a secondary particle which is emitted 
opposite a small angle. The uncertainties introduced in this energy region 
exceed those due to straggling, namely + 1.0 MeV at 48 MeV. 

Table I shows the uncertainties set out as a function of energy for the 
energies of the negative pions of the present sample. 


TABLE II. 


Buergy 


= È — 
saci 0+4/4—8 sa ~ 16/16 — 20/20 — 24/24 — 28/28 - 32 
(MeV) | 


es: 


Tei 
32- = 40 40-+ 4444 + 48] 


5: 


5% | 


| 
| 
| 
| 
| 
| 
| 


BEE 
2% 3% 


le |10%| 6% | 4% | 3% | 3% | 3% | 3% | 3% 


90/ 
2% 


3.7- Diseussion. 


31. Energy distribution. — The energy and angular distributions have been 
calculated by DALITZ and FABRI. They are based on the assumption of an 
interaction radius roughly equal. to the Compton wavelength of the +-meson 
and an interaction independent of the orbital momentum quantum numbers /, 1’ 
(defined in reference (1)). 

Coulomb effects and pion pion interactions in the final state are not con- 
sidered. 

Then, neglecting the possible influence of unknown partial or complete 
selection rules, the lowest order term in the expansion in terms of angular 
momentum eigenstates is by far the most important. The effects of this term 
on the energy distribution of the t-meson decay are displayed by these authors 
for the spin and parity combinations 0- 1+ 1- 2+ 3-. In Fig. 2 the experimental 
energy distribution for the negative pion is compared with the theoretical 
curves in the non-relativistic approximation. Also shown in Fig. 2 is the 
relativistic curve for the 0- configuration. Energy distributions calculated in 
this way show a maximum sensitivity to the value of the spin in the ex- 
treme parts of the spectrum, since in these regions they are less dependent 
on the simplifying assumptions made. 

Because of conservation of angular momentum the low energy negative 
pions must have a momentum dependence p" where n> 2 for +-mesons of 
odd spin and can be zero for a t-meson of even spin. Using the world data 
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available at the time, OREAR (8) showed that the results were consistent with 
n=0 and strongly inconsistent with n=2. This was done by calculating 
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Fig. 2. 


the probability that the low energy events occurred as they did, assuming 
1+ to be the true spin-parity, divided by the probability that they occurred 
as they did, assuming 07 
to” be the true spin-parity 


assignment. Fr 7% ue 
In the present work, the sole 
energy of the negative pion $ 
was found to be less than 30} 
10 MeV in 56 cases. The 3 3 
= 


probability of spin zero to 


spin one was found to be 104. 10 
We can infer that it is ex- ott , PI rr > Tv 
tremely unlikely that the pesa 

ig. 3. 


<-meson has spin 1 or any 
odd value. It does not seem 
possible to exclude spin 2 for the t-meson by means of a Dalitz type ana- 


lysis alone. 
It has been found that the energy spectra of the negative and positive 


(8) J. OREAR: Phys. Rev., 106, 834 (1957). 
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r-mesons are not identical and show departures from that calculated from the 
variation with energy of the relativistically computed statistical factor. This 
effect has been observed also by BALDO-CEOLIN et al., with their sample of 


T° Spectrum. 
899 Events 


Number of events 


24 
Fig. 4. 


40 MeV 48 


32 


400 decays. Since the present 
work was performed with 
criteria similar to theirs, the 
results of the two large sam- 
ples thus collected have been 
combined. The combined 
results of the 899 7-meson 
decays are displayed in Fig. 3 
(x--spectrum) and fig. 4, (x 
spectrum). In Table III, are 
given the coefficients of 
asymmetry and excess (7, 


and y, respectively) relative to the mean of the statistical factor for the spectra 
of Figs. 3 and 4, and for the statistical factor. 


TABLE III. 
Spectra of 499 Events 
mt Ca Theoretical 

Lo 145 148 144.34 

Ms hee 625 + 118.94 

La 42588 42681 41 850 

Va 0.040 7 — .3471 — 0.0686 + 0.3137 (*) 

Wey — 1.0044 — 1.0546 — 0.9913 + 0.1058 (*) | 

Spectra of 899 Events 
TT TU Theoretical 

Us 139.36 144.50 144.34 

Us — 73.88 534.23 — 118.94 

Ma 39 210 41 755 41 850 

y 0.0449 — 0.3075 — 0.0686 + 0.3137 (*) 

Va — 0.9815 — 1.0004 — 0.9913 + 0.1058 (i 
us = moment of order 2 with respect to the theoretical mean, for the statistical distribution ; 
#3 = moment of order 3 etc.; 
HA = moment of order 4 etc.; 
y. = coefficient of skewness; 
y: = coefficient of excess. 

(*) The errors given in coiumn 4 have been calculated on the assumption that the indi- 
vidual pion energies have an error of 1 Mey. 
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The table shows that the experimental distributions depart somewhat from 
that predicted by the relativistic statistical factor (r.s.f.). For the spectrum 
of the negative pions, the departure is towards the higher energies, and for 
the positive pions, is towards the lower energies. These tendencies can be 
more easily appreciated in Figs. 5 and 6 where in each diagram the r.s.f. is 


denoted by a horizontal straight line. Each ordinate is the ratio of the 
number of secondaries found experimentally with energies in a particular 
energy interval to the expected number of secondaries deduced from the same 
energy interval of the r.s.f. 

To the sets of ordinates thus obtained, the following straight lines were 
fitted, using the method of least squares 


Y(r-) =( 0.0105 + 0.001 64 )X ++ (0.7578 + 0.047 4), 


Y(x+) = (— 0.005 07 + 0.000 615)X + (1.1851 + 0.0165). 


The errors shown are the mean quadratic errors to be feared. It has been 
suggested by MITRA (*) that similar asymmetries in the energy spectra may 
be due to resonant pion-pion interactions, thus bringing about a certain mo- 
mentum dependence in the various transition amplitudes. However, as pointed 
out by DALITZ (1) very little is known about the strength of 7-7 forces. 

It is also possible that the discrepancy could be accounted for in terms of 
a difference between the non-relativistic and relativistic matrix elements, and 
that in extending the calculation to the relativistic case it is not a sufficiently 
good approximation to change only the statistical factor. 


(9) A. N. MITRA: Nucl. Phys., 6, 404 (1958). 
(0) R. H. DALITZ: Rep. Progr. Phys., 20, 163 (1957). 
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4. — Conclusions. 


A homogeneous sample of 540 7-decays has been studied under the strictest 
conditions of scanning and with every safeguard against appreciable bias. 
A comparison of the results with the curves calculated by several authors 
shows that the t-meson decays in a state of zero spin and odd parity. The 
sample has been combined with the events of BALDO-CEOLIN et al. and the slight 
departure from symmetry in the distribution of the decay energy among the 
three pions already noted by these authors has been confirmed. 


* * OK 
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APPENDIX I 


The following is a list of the energies of the negative pions from the 
499 +-meson decays used in this analysis. Those which have been derived 
from angular measurements are marked with a cross (+); those from the energies 
of the two positive pions with an asterisk (*); the remaining pions were traced 
to rest, and their energies deduced from the measured range. 


TABLE IV. 
IRR) 01 5.19+ 7.08 7.8* 9.3 10.45+ | 11.09+ | 12.62 1321 
Tora 3.61+ 5.2 Tel 7.99 9.3 10.50 11.38 12.75+ | 1314 
1.7 D 12 5.49 eni 8.3 9.33+ | 10.60+ | 11.75 12.76+ | 13.6 
2.0 3.8 5.65 7 ASS 8.4 9.5 10.7 11.8 12.85 1857 
2.054 3.82 5.8 7.16 8.6 9.8 10.8 11.8 12.90 157 
2.4 4.00+ 5.86+ 7.19+ 8.7 10.05 10.95 12.04+ | 12.9 13.89+ 
2.6 4.13 6.25 re) 8.9 10.2 11.0 12.10+ | 12.9 13.99 
2.85 4.14+ 6.32 7.4 8.9 10.20+ | 11.0 T2.137 12/994) 414987 
3.40 4,9 6.75 7.45 8.99+ | 10.3 RES PAT 13.0 14.30 
3.42 5.03+ 6.84+ Heel 9.2 10.38+ | 11.03+ | 12.6 13.04+ | 14.41+ 
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TABLE IV (continued). 
14.53+ | 15.4 16.2 | 16.85+ 1793S") 18377 ~|, 20:0 20:70+ 121.5 22.01* 
| 14.60 15.45+ | 16.2 16.9 18.0 18.90+ | 20.0 20.777 | 21-5 22.1 
14.86+ | 15.68+ | 16.24+ | 17.0 18.1 18.98+ | 20.22 | 20.8 21.555 | 22.13+ 
14.87 15.78+ | 16.3 Il 18.2 19.2 20.28+ | 21.06* | 21.56+ | 22.32+ 
15.0 15.8 16.98E 17134 4) 18.2 19.2 20.3 21.08+ | 21.6 22.5 
15.02+ | 15.9 16.46+ | 17.2 13.2* 19:7 +, 20.3 Des 21.6 22.60+ 
| 15.1 15.9* 16.6 17.6 18.27+ | 19:7 20.41+ | 21.3 21.6 22.63 
| 15.13+ | 15.9 16.7 17.64+ | 18.3 19.8 20.5 2140 21.8 DOTT 
| 15.2 L'6-13% | 16276" | .17.8 18.52 19.91+ | 20.5 21.44% | 21,94+ | 23.0 
| 15.4 | 16.2 16.80 | 17.80+ | 18.56 | 19.92+ | 20.6 21.5 21.99 | 23.07+ 
| 23.1 | 24.09 24.68+ | 25:0811#20-087 414971 27.9 28.43+ | 29.00+ | 29.80+ 
23.35 | 24.10 | 24.75+ | 25.92+ | 26.74+ | 27.2 27.9 28.44+ | 29.08+ | 29.81+ 
23.44+ | 24.3 24.8 25.95+ | 26.75* | 27.4 28.05 | 28.5 29.14+ | 29.85* 
23.500 242335 | 24.9 25.99+ | 26.8 27.50+ | 28.07+ | 28.5 29.29+ | 29.92+ 
| 23.50 | 24.42+ | 25.0* 26.0 26.89+ | 27.5 28.10+ | 28.54+ | 29.3 30.02+ 
23.53+ | 24.46+ | 25.0 26.18+ | 26.9 27.5 28.11+ | 28.66+ | 29.5 30.05 
23.6 24.46 25.2 26.30 | 26.9 27.59+ | 28.16+ | 28.70 | 29.5 30.24+ 
2307 24.48+ | 25.41+ | 26.4 27.0 27.69 | 28.18+ | 28.75+ | 29.5* 30.37+ 
23.8 24.50 | 25.5 | 26.5 25.0 SISI 28.3 28.8 29.52+ | 30.5 
| 24.0 24.5 25.5 26.5 27.05 | 27.8 28.4. 29.0*, | 29.6 30.5 
| 30.76+ | 31.5 32.03+ | 33.15+ | 33.99+ | 34.43* | 35.26+ | 35.8 36.49+ | 37.0* 
30.79+ | 31.5* 32.04+ | .33.2* 33.99 | 34.5 30.00" | 3028 36.5 Silk) 
30.8 31.58 * 532.2 Bors 34.04 | 34.5 90.9 35.84+ | 36.5 37.01 
30.8* 31.64+ | 32.39+ | 33.35+ | 34.1 34.6 en 90.0) 36:0") 37.3 
| 31.0 31.66+ | 32.4 39.5 34.17+ | 34.68* |.35.4 36.01+ | 36.5* 37.4 
31.0 DL. Gus aloes 30.0 34.19 34.87+ | 35.4 36.1 36.6 37.5 
31.04+ | 31.76+ | 33.0 33.5 34.19* | 35.0 35.41* | 36.10+ | 36:6 37.5* 
SO NS 77 300533770 224217 35.0 35.5 36.23+ | 36.74 | 37.5 
31.4 31.8 33.02+ | 33.8 34.82* | 3516+ | 35:75 | 36.4 36.98+ | 37.67+ 
3105 31.88 33:10" 133.9 34.4 202101035760 |), 36.40" Wo m0 SI 
37.76+ | 38.90+ | 39.7 40.62* | 41.6* 42.5* 43.56+ | 44.08* | 44.70* | 45.50+ 
37.9* 38.94* | 39.70 40.69* | 41.70+ | 42.5* 43.58+ | 44.09* | 44.87+ | 45.73+ 
37.94 | 39.0* 39.79+ | 41.0 41.83+ | 42.5* 43.62+ | 44.13+ | 44.9*° | 45.8 
38.29+ | 39.0 39.81* | 41.0 41.9 42.5 AS TOs Ne 44.0%) 45.0 4.6.1 
38.4 39.0* 39.95+ | 41.0 42.2 42.54+ | 43.81* | 44.2* | 45.02*| 46.2* 
38.41+ | 39.02+ | 39.99 41.2 AD. 2s) ADT 43.84+ | 44.34+ | 45.05* | 46.51+ 
38.5 39.09+ | 40.2 41259 ayn Aone 43.04+ | 43.9* 44.5* | 45.06+4 47.0* 
38.5 39.18+ | 40.46+ | 41.34+ | 42.4* 43.20 | 44.0* | 44.57* | 45.30+ | 47.0 
98.594 39.25) |040.5* 41.5 42.5 43.4 44.01 44.59+ | 45.39+ —- 
38.6 39.5% | 40.82% |-41.58+ |, 42.5* | 43.5 44.06+ | 44.62+ | 45.5* 7.79 
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APPENDIX II 


In the first event which cannot be fully analysed, one 7* has an energy 
15.4 MeV; another pion whose charge is unknown because its track ends 
between pellicles, has an energy 36.1 MeV; the third pion interacts. The energy 
of the negative pion is either 23.4 MeV or 36.1 MeV. 

One + from the second 7-decay has an energy 3.1 MeV and one pion 
interacts. 22.00 mm of the third track can be seen before it leaves the avail- 
able pellicles and the residual range indicated by ionization measurements 
is 0.54mm. Hence the energy of the negative pion is either 37.6 MeV or 34.2 MeV. 


RITA SS DENON) 


Si è studiato un campione omogeneo di 540 decadimenti 7 trovati in pacchi di 
emulsioni esposti al bevatrone di Berkeley. La maggior parte degli eventi è stata 
individuata osservando lungo la traccia e si sono resi minimi gli effetti dovuti a fattori 
geometrici e d’osservazione. Si discutono i risultati in dettaglio e si combinano poi 
con quelli forniti da un campione di volume comparabilmente grande studiato da 
BALDO-CEOLIN et al. (7). La più probabile assegnazione di spin-parità è 0- e abbiamo 
potuto confermare e valutare le deviazioni degli spettri dell'energia positiva e negativa 
dei pioni dai valori predetti dal fattore statistico relativistico, già osservate da BALDO- 
CEOLIN et al. 


(*) Traduzione a cura della Redazione. 
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Conversion Coefficient of the 87 keV Transition in el Die 


M. ©. Josm, B. N. SuBBA RAo and B. V. THOSAR 
Tata Institute of Fundamental Research, Bombay 


(ricevuto il 18 Agosto 1958) 


Summary. — The K-conversion coefficient of the 87 keV, H2 transition 
in }°Dy has been determined from the photoelectron spectrum and the 
internal conversion spectrum of 160Tb investigated with the intermediate 
image f-ray spectrometer. «* is found to be 2.0 + 0.2 which is large 
compared with the theoretical value. 


1. — Introduction. 


The decay of 1°Tb (7, = 73 d) has been studied by many workers (16) and 
there is good agreement regarding the low-lying excited states in 1 Dy. From 
these studies, it is well established that the first excited state is at 87 keV 
and that, in accordance with the systematics of even-even nuclei, this level 
has 2 + assignment for its spin and parity. So, the 87 keV transition (2+ > 0+) 
Should be purely £2. 

There is no accurate determination of «“ reported so far for this 87 keV 
transition and there seems to be disagreement among the values given by 
previous workers. Recently, McGowan et al. (7) found that experimental 


(1) S. B. Burson et al.; Phys. Rev., 94, 103, (1954). 

(?) V. KESHISHIAN et al.: Phys. Rev., 96, 1050 (1954). 

(3) G. BERTOLINI, M. BeTTONI and E. LAZZARINI: Nuovo Cimento, 3, 754, 1162 (1956). 
(4) M. A. CLARK et al.: BAPS. II, 4, 231 (1957). 

(5) O. NATHAN: Nucl. Phys., 4, 125 (1957). 

(9) S. OFFER: Nucl. Phys., 5, 331 (1957). 

(€) F. K. McGowan and P. H. Stenson: Phys. Rev., 107, 1674 (1957). 
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values of « for a few pure Æ2 transitions were consistently higher than the 
corresponding theoretical values. So, it was felt necessary to determine the 
conversion coefficient of this transition as accurately as possible. 


2. — Experimental methods. 


The Tb source was obtained by neutron irradiation in the Harwell pile. 
A small amount of the powder was used as source. The diameter of the 
source was 0.4 cem. Sufficient amount of perspex was used to absorb primary 
electrons from the source. Tin, deposited 
by vacuum evaporation on a perspex 
dise was used as converter for obtaining 
the photoelectron spectrum. The diameter 
of the converter was about 2cm and 
the source-converter distance was about 
0.4cm. Thickness of the converter was 
about 0.3 mg/em?. 

The photoelectron spectrum obtained 
with such an arrangement in the inter- 
mediate-image -ray spectrometer is shown 
in Fig. 1. K- and L-photo-electron line 
(ler Up, and #7, Lp,) of K,- and Kg-X-rays 
of Dysprosium and the K-photoline (y,) 
of the 87 keV y-ray are seen in it. Back- 
ground was determined in another run 
under identical conditions but without the 
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(N/Be «7) > 


Fig. 1. - The photoelectron spectrum 
of y-rays emitted in the decay of 


converter. 
The internal conversion spectrum of the 
low-energy portion, also observed in the 


160Tb using tin as radiator (low 


intermediate image 8-ray § 
ARA ge f-ray spectrometer 


using a thin source, is shown in Fig. 2. 
The most prominent lines were the X and 
(L+M) conversion lines of the 87 keV transition, represented by K, and 
(L1-+M,) in the figure. K, and (Z,+ M,) represent the conversion lines of the 
197 keV transition. The latter are shown on a larger scale in the inset. Other 
conversion lines could not be clearly brought out indicating thereby that they 
are very weak compared even with K,. 
The area under the photoelectron peak may be represented (8) by 


(1) P=k-Ix0, 


i 


(8) A. C. G. MircHELL: Beta- and Gamma-ray Spectroscopy (1955), p. 229. 
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ces I is the intensity of the incident electromagnetic radiation, 7 is its photo- 
electric absorption coefficient, ¢ is the thickness of the converter, C is a factor 
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Fig. 2. - Internal conversion electron spectrum (1#Tb). 


representing the effect of angular distribution of photoelectrons with respect 
to the incident photon momentum and k is a constant which depends on the 
transmission of the spectrometer and geometry of the source-converter assembly. 
Then the ratio of the areas under the photoelectron peaks due to X-ray and 
y-ray may be written, 

Px ma Ix tx Cx 


(2) PAPI 
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From this one can determine 7,/1,, using theoretical values for Tx and T,, 
provided C,/C, is known. The energy of X-rays of Dysprosium is soe 
50 keV and the energy of the y-ray is 87 keV. The corresponding separation 
of the maxima of intensity distributions of photoelectrons with angle is ex- 
pected theoretically to be about 10°. Further, such an angular distribution 
for low photoenergy is a fairly slowly varying function (°) of # around Pj 
for a given energy of photoelectrons. The measurement of the K conversion- 
coefficient (1°) of the 144 keV transition in ™'Pr;,,, by two independent methods, 
one based on B-spectrum and conversion line intensities and the other one 
being the same as the one adopted here, led to the same value within exper- 
imental errors though the separation between the maxima was about 20°. 
These considerations indicate that, for the experimental arrangement used 
here, CEE ales 


8. — Results and discussion. 


By adopting the procedure described above, /x/I,@, was found to be 
(2.2 + 0.2), where ©, is the K-fluorescence yield (!) of Dysprosium (@, = 0.93) 
But the X-ray intensity determined above is not due to the internal conver- 
sion of 87 keV transition alone. In order to account for this, the K-conversion 
line intensity of this transition relative to the total of all K-conversion line 
intensities occuring in the @--decay of 1Tb was used. From our observations 
of internal conversion spectrum, it was found that K,/(K,+4,) = 0.94. Using 
the conversion data of O. NATHAN (5), M. A. CLARK et al. (4) and M. P. Avo- 
TINA et al. (1%) also K,/K was found to be 0.91 + 0.03. Thus, the right value 
of 1 jlo, = aX — 2.0 + 0.2. 

The theoretical value due to Rose ('4), based on the «no penetration » 
model which takes into account all static effects, for the K-conversion coef- 
ficient of the 87 keV #2 transition in Dysprosium is xf = 1.53. The experimental 
value is clearly larger than this. This gives for the ratio «%/a*,,. = 1.29+0.13. 
McGowan et al. (7) in their study of the K-conversion coefficient of H2 transitions 
—the 80 keV transition in Er, the 84 keV transition in 1”Yb and the 89 keV 
transition in ‘*Hf—find this ratio to range from 1.13 to 1.23. We also find 


(*) C. M. Davisson and R. D. Evans: Rev. Mod. Phys., 24, 79 (1952). 
(1) M. C. JosHI, B. N. SusBA Rao and B. V. THosar: Nuovo Cimento, 9, 600 (1958). 
METE 


DAVISSON: Beta- and Gamma-ray Spectroscopy (1955) p. 836. 
I. BERGSTROM: Beta- and Gamma-ray Spectroscopy (1955) p. 630. 
M. P. Avorina et al.: Doklady. Akad. Nauk. SSRR, 119, 1127 (1958). 
M. E. Rose: Privately circulated tables. 
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the K-conversion coefficient of the 84 keV 2 transition in 190Yb, deter- 
mined by the method described in this work, to be nearly the same as that 
reported by McGowan et al. (7). These determinations indicate that exper- 
imental x are larger than the corresponding theoretical values by about 20%. 


* OK OK 


Our thanks are due to Mr. C. V. PANAT and Mr. S. D. BHAGWAT for their 
assistance during this work. 


RIASSUNTO. (*) 


Il coefficiente di conversione K della transizione H2 di 87 keV del 1®Dy è stato 
determinato dallo spettro dei fotoelettroni e lo spettro di conversione interna del 19°Th 
è stato studiato con lo spettrometro f a immagine intermedia. Si trova che af è 2.00.2, 
valore alto rispetto al valore teorico. 


(*) Traduzione a cura della Redazione. 
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The Effect of the Final State Interaction on the Decay of Hyper-*He (*). 


Yi SS IAN 440) 
University of Illinois - Urbana, Ill. 


(ricevuto il 18 Agosto 1958) 


Summary. — The effect of the final state interaction on the decay of 
hyper-5He has been studied. It is found that the strong interaction 
between the «-particle and the proton produced in this decay plays a 
significant role in determining the energy spectra and the angular distri- 
butions of the final particles. With the «-p force taken into account 
as determined by the experimental phase shifts, a simplified calculation 
indicates that the pion is emitted in the energy range from about 21 to 
30 MeV for approximately 90% of the events, and that the proton is 
scattered predominantly in the backward direction in the «-p center of 
mass system. These findings are in fair agreement with Hills analysis 
of recent experimental observations on the decay of hyper-5He. 


1. — Introduction. 


A recent study by HILL (!) has shown that the decay patterns are quite 
similar for all of the known cases of « uniquely defined » decays of hyper-*He. 
The energy of the pion produced in this decay shows a strong tendency to lie 
within the range approximately from 21 to 29 MeV, and for eight of the nine 
cases studied, the proton is projected in a backward direction in the center 
of mass of the «-p system. In view of these observed regularities, HILL there- 
fore suggested that the strong interaction between the «-particle and the proton 
must play a significant role in the decay of this hyper-nucleus. 

If the interactions among the final particles ‘He, p and ~~ produced in 
this decay are completely neglected, then the energy spectra and the angular 
distributions of those particles can be readily calculated. However, the results 
so obtained are quite inconsistent with the experimental observations. Instead 


(*) Assisted by a joint program of the U.S. Atomic Energy Commission and the 
U.S. Office of Naval Research. 


(*) Now at the Florida State University, Tallahassee, Florida. 
(1) R. D. Hiri: Nuovo Cimento, 8, 459 (1958). 
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of a prominent backward Scattering asymmetry, the theoretical calculation 
indicates that the proton should rather be scattered in a forward direction in 
the center of mass of the «-p system. In addition, the energy spectrum of 
the pion obtained in this manner does not show any distinct preference for the 
region above 21 MeV. These discrepancies, therefore, necessitate the consi- 
deration of the final state interactions in the decay process. 

According to the theoretical study of WATSON (2), three conditions would 
have to be met if the final state interaction is to seriously modify the cha- 
racteristics of a reaction. These conditions are a) that the primary reaction 
must be confined to a certain volume V of order a, where a is the range of 
the primary interaction, 6) that the relative energies of the particles whose 
interaction is under consideration must be comparatively low, and €) that 
the final state interaction must be strong and attractive. For the decay of 
the hyper-5He, the primary reaction, which is the free decay of a A°, takes 
place in a sphere with a radius of the order of 0.5-10-!* cm. In the «-p system, 
for which the interaction is most important, the relative energy of interest 
is around 10 MeV, corresponding to a wave number % of about 6-10! em-1. 
The magnitude of ka is therefore only about 0.3, which is small compared to 
unity. As for the third condition, it is well known that, although the s-state 
interaction in the «-p system may be repulsive, the p-state interaction is strongly 
attractive. Thus, all three conditions are satisfied; consequently, we expect 
the final state interaction to be significant. 

It is fortunate that the influence of the final state interaction is calculable 
in our case. The «-p force may be inferred from the experimentally determined 
phase shifts, which have been listed quite extensively for relative energies up 
to about 20 MeV (*%). The interaction between zm and « can be neglected, 
since in the decay of hyper-5He, the energy of the pion never exceeds 33.4 MeV. 


2. — Formulation of the problem. 

The matrix element U which determines the probability of the decay mode 
(1) ‘He, + ‘*He.+ p:-- 17 
may be written as 


ee d . i 7 o-K 
ey U=|yty"rexp —iK |r. =} + 4a) (s+ 2%] tide’ 


(2) K. M. Watson: Phys. Rev., 88, 1163 (1952). 

(8) C. L. CritcHFIELD and D. C. DoppER: Phys. Rev., 76, 602 (1949). 
(4) D. C. Dopper and J. L. GAMMEL: Phys. Rev., 88, 520 (1952). 

(5) K. W. Brockman: Phys. Rev., 108, 1000 (1957). 

(8) P. E. Hopeson: Adv. in Phys., 7, 1 (1958). 
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where o denotes the baryon spin vector, while s and p denote the amplitudes 
for the 1. =0 and J, =1 channels in free A-decay, respectively. The vector r 
is defined as {r,—r,), where r, and r, are the position vectors of the A-par- 
ticle and the «-core. AK is the momentum of the pion in the three-body decay 
mode of 5He,, and AK, is its value in the free decay of a A-particle. Finally, 
y,, y and y; represent the wave functions of the «-core, the x-p system and 
the hyper-nucleus 5He,, respectively. 

Since the binding energy of the A-particle in the hyper-"He is rather low, 
it is appropriate to assume the wave function to be in the form of a product; 
that is, 


(3) Vi = VM), 


where w\(r) represents the orbital motion of the A-particle relative to the 
a-core, and y, is the spin function. 

If we make a further assumption that the «-core is not distorted by the 
presence of the A-particle (7), then Eq. (2) may be simplified, yielding 


(4) feos | y* (7) exp[— 20K -r] ( a a anya 


with b equal to #. 
For the x-p system, the wave function must be so chosen that it represents 
a plane wave modified by an ingoing spherical wave (8). Taking into account 


the spin-orbit force, we may therefore write y(r) as 


(5) y(r) = 5 d'[r(21 +1)}?- 


1=0 
{C14 (04-3, M; 0, M)[H (kr) + exp [— 2187] h? (kr)] Vi à y+ 


+ a, M; 0, MAL (kr) + exp[— 2467 ]h? (kr)] Vi 13}, 


in which € is the Clebsch-Gordan coefficient and h, is the spherical Hankel 
function. #k is the relative momentum in the «-p system, Y is the eigen- 
function of the total angular momentum (*), and è, is the anomalous phase 
shift which has been determined extensively for various energies ranging up 
to around 20 MeV (**). In this work, we shall use only the terms involving 


(*) According to R. H. Dazrrz (to be published), the optimum compression of the 
a-particle core is only about 5% for the hyper-nucleus 5He. 

(3) G. BREIT and H. A. Berne: Phys. Rev., 98, 888 (1954). 

(°) J.M. BLATT and V. F. Weisskorr: Theoretical Nuclear Physies (New York), p. 790. 
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1-0 and /=1 in Eq. (5), since in the «-p System, the relative energy of 
interest does not exceed 20 MeV, at which the values of 0, for 1>2 are 
quite small. 

For the wave function y(r) shown in Eq. (5), we have not considered the 
repulsive Coulomb interaction, since a proper account of it by the modified 
Coulomb wave function would involve analytical difficulties in the integration 
procedure that follows. We, therefore, choose to use at this point the plane- 
wave approximation, keeping in mind that a modification has yet to be made 
after the integration process is completed. 

To find the transition probability of the decay mode (1), it is necessary 
to take the square of the matrix element U, and then average over the spin 
States and the outgoing directions of the pion. Using Eqs. (4) and (5), and 
with due consideration of the volume in phase space, we get for the transition 
probability per unit pion energy per unit solid angle 


a7 


(6) dE-dQ, 


2 
D (s+ pte) KL (+ LI, 
KA 
where D is a constant, and 
(7) 1, =| pi) exp[—WK-r] ug (rar, 


(8)  y,(r)= exp [ikr] — 3 i+! [4x/(21 + 1)]? Vo h- 


‘[( + 1) sin d* exp[—id#]+! sin 6; exp [—id,]], 


Lee] 


@) ys(r)= > i [7/2 + 1)} [Ud + 1} Vth” {exp [— 24] — exp [—2i6; J}, 


1=0 


(0) ys(r)= > Pret +1) [Ut + DIE Y; Th? {exp [ —240;*] —exp[— 269, ]}. 


1=0 


The wave function w,(r) will be obtained by using a variational procedure. 
Since the A-particle is known to be in the s-state, a trial function of the 
Hulthen form 


(11) u(r) = N(exp [— ar] — exp [— fr])/r 


will be employed, in which « is related to the binding energy of the hyper->He, 
and f is the variational parameter. Assuming a K-exchange interaction for 
the A-n force, and using the results of the Stanford experiments (1°) on the 


(19) R. HorsraptER: Rev. Mod. Phys., 28, 214 (1956). 
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shape and radius of the «-particle core, the potential seen by the A-particle 
is then 


(12) V(r) =— VY, exp [— (r/2)*], 


where R is equal to 1.31:1071* cm. 

Since the binding energy of the A-particle in *He, is known, we can obtain 
by the variational method (11) the value 6 = 3.4%, which corresponds to V, of 
about 57 MeV. The total volume integral for this particular well depth jis 
then found to be 720 MeV f*, a value in good agreement with the results of 
other authors, in particular, that of DALITZ (12). 


3. — Energy spectra and angular distributions. 


The energy spectrum of the pion may be obtained by integrating Eq. (6) 
over the element of solid angle 49, in the center of mass of the «-p system. 
With the use of Eq. (11) for u,, the evaluation of AT/dE_ is straightforward. 
As the wave function y used in the integrals does not take into account the 
repulsive Coulomb force between the «-particle and the proton, the spectrum 
so obtained still needs to be corrected for this interaction. In view of the 
fact that the Coulomb field accelerates the proton away from the «-particle, 
we must, therefore, represent its effect in modifying the spectrum of the pion 
by including the Coulomb correction factor F (1%), and by observing that in 
Eq. (5), k is to be replaced by [k? —(2u/#2) V(r.) ]*, in the presence of the Coulomb 
interaction, where V(r) is the Coulomb ‘potential between the «-particle and 
the proton, and w is the reduced mass. As a reasonable estimate, we shall 
take for 7, the radius at which the wave function w, has its maximum value; 
consequently, V(r) is found to be approximately equal to 1.7 MeV. However, 
it is to be emphasized that treating the Coulomb interaction in this way is 
not very accurate, especially when the relative energy in the «-p system is 
low. Thus, for energies of the pion exceeding about 31 MeV, the values of 
AT/AE_ so obtained may be quite inaccurate, and should be taken only as an 
indication of the order of magnitude. 

The energy spectrum of the pion is shown in Fig. 1. The curve a) is for 


(1) See, for example, R. G. Sacus: Nuclear Theory (Cambridge, Mass.), p. 33. 


(12) R. H. DALITZ: Proc. of the Sixth Annual Rochester Conference on High-Energy 
Nuclear Physics (New York, 1956), Sect. 5, p. 40. 


(*) For the meaning of the Coulomb correction factor F, see ref. (°), p. 680. 
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the case in which the final state interaction has been taken into consideration, 
and the curve (b) is for the case of no interactions () among the final par- 
ticles. Both curves have been plotted with p set equal to zero. 
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Fig. 1. - The energy spectrum of the pion: (a) with final state interaction and (b) without 
final state interaction. 


To arrive at the energy spectrum of the pion, the expression (5) for y, 
while valid only in the region outside of the range of the «-p potential, has 
been used throughout the entire domain of integration. The accuracy of our 
result is, however, not greatly impaired by this approximation since the con- 
tribution to the integrals comes mostly from regions of space larger than that 
for which the «-p potential is important. 

When the A-particle lies within the core nucleus, it is likely to decay by 
the non-mesonic mode, with the consequence that no pion appears among the 
final particles. To show the effect of this mode on the shape of the pion energy 
spectrum, we have evaluated J; of Eq. (7) by integrating from a radius R, 
(taken to be approximately equal to the intrinsic range of the A-« potential, 
i.e., 2-10-13 em) to infinity; the use of R, as a limit of integration comes from 
the consideration that when the A-particle is within the sphere of radius R,, 
non-mesonic decay mode will probably dominate. The energy spectrum of 
the pion calculated in this manner exhibits the same general features as those 
for the case in which the non-mesonic decay mode is neglected, and hence will 
not be given any further attention here. 


(4) By no interaction among the final particles, it is meant that we set the 
anomalous phase shifts 5, equal to zero. 
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We shall now determine the angular distribution of the final particles. 
Let 0 be the angle between the direction of emission of the proton in the 
center of mass of the «-p system and the direction opposite to that of emis- 
sion of the pion, and let us define 


ae) 


d27 ; 
DE; Se Sin 0 dû ; 


ar a 
=| ayaa, 8%, 


12 


then with the use of Eq. (6), the ratio Z,/1; can be found quite readily. As is 
clear from the definitions of I; and 1, expressed by Eq. (13), the ratio 1,/I; 
has the following meaning: for /,/I; greater than unity, the angle, 0, for a 
given pion energy, will probably be between 90 and 180 degrees; for J,/I, 
less than unity, 0 will likewise lie between 0 and 90 degrees. 


zi E, in MeV 
WA ee Der 6 To aie CE a 


Fig. 2. — The ratio of probabilities for backward to forward scattering of the proton 
in the center of mass of the «-p system: (a) with final state interaction and (b) without 
final state interaction. 


The ratio I,/I, is plotted against the energy of the pion in Fig. 2, again 
with the curve (a) for the case in which the final state interaction is consi- 


dered, and the curve (b) for the case of no interactions among the final par- 
ticles. 
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4. — Discussions and conclusions. 


The energy spectrum of the pion as shown by the curve (a) of Fig. 1 shows 
that for approximately 90% of the events, the pion is emitted with an energy 
in the range from about 21 to 30 MeV. This is in good agreement with Hill’s 
analysis of recent experimental observations, which indicates that in only one 
of the nine events studied, the pion has an energy less than 21 MeV. On the 
other hand, the curve (6) which is calculated by setting the anomalous phase 
shifts equal to zero, demands that the pion should be emitted with an energy 
not exceeding 21 MeV for about 30% of the events, a conclusion not in good 
accord with the results of the experimental findings. | 

The distinct peak near 30 MeV in the calculated energy spectrum of the 
pion indicates that there should be à number of decay events of 5He a recorded 
with the pion energy in the neighborhood of this value. However, only one 
event was reported. This discrepancy may therefore be used as a strong ar- 
gument against the validity of our theoretical calculations. It is to be noted, 
however, that this disagreement might conceivably come from the fact that 
it is difficult to observe low energy «-particles in photographic emulsions. 
For instance, with H_ equal to 30 MeV, the energy of the a-particle is 
[1.48 + 1.47 cos (180°— 6)] MeV. Since at this pion energy, /,/I; is equal to 
about 0.6, 0 is likely to be less than 90 degrees. Therefore, the energy of the 
a-particle may be quite low; for example, with 0 = 30°, the energy is 0.2 MeV, 
which will leave a track length less than a micron in the emulsion. In fact, 
some of the helium hyperfragment decays which are now classified as non- 
uniquely defined cases of *He,, because the recoil He tracks are too short for 
accurate identification, could be cases of relatively high energy pion emission. 
It is quite probable that the inclusion of these cases would bring the exper- 
imental distribution into better agreement with the calculated curve. 

According to the curve (a) of Fig. 2, the ratio /,/I; has an average value 
of about 1.35 in the range of pion energy from 21 to 28 MeV. On this basis, 
the expected number of decay events with 0 larger than 90° is then 5.2, which 
is only slightly more than one standard deviation from the observed number 
of 8. On the other hand, within the same energy range, the average value 
on I,/I; for the curve (b) is only about 0.1, which is in total disagreement with 
the experimental evidence. 

In view of the above discussions, we conclude that the final state inter- 
action must play an important role in determining the energy spectra and 
the angular distributions of the final particles in the decay of the hyper-5He. 
When this interaction is considered, the decay characteristics of this hyper- 
nucleus can be reasonably well understood. In addition, it is worth noting 
that the shape of the calculated energy spectrum of the pion may permit the 
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determination of the ratio p/s of the s- and p-channel amplitudes in free 
A-decay. However, this determination does not appear feasible at present, 
since the number of ‘He, decay events reported is still too few to justify such 


a calculation. 
* K * 


The author wishes to thank Professor R. D. Hu for suggesting the problem 
and for helpful discussions. He is indebted to Professor H. W. Wy tp for his 
continued encouragement and guidance. Valuable advice by Professor J. D. 
JACKSON is also gratefully acknowledged. 


RIASSUNTO (*) 


Si è studiato l’effetto dell’interazione dello stato finale sul decadimento dell’iper-* He. 
Si trova che l’interazione forte tra la particella x e il protone prodotto in questo deca- 
dimento ha una parte importante nella determinazione degli spettri d’energia e delle 
distribuzioni angolari delle particelle finali. Tenendo conto della forza «-p che si deter- 
mina in base agli spostamenti di fase sperimentali, un calcolo semplificato indica che 
il pione è emesso nell’intervallo d’energia da circa 21 a 30 MeV per circa il 90% degli 
eventi, e che lo scattering del protone avviene prevalentemente in direzione retro- 
grada nel sistema del centro di massa dell’x-p. Questi risultati sono in buon accordo 
con l’analisi di recenti osservazioni sperimentali sul decadimento dell’hiper-5He fatta 
da Hill. 


(*) Traduzione a cura della Redazione. 
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A Note on the Inelastic Contributions to the Elastic Scattering 
of High-Energy Electrons from Helium (*). 


P.. 0. DAyzy (*) and H. 8. VALK 
Department of Physics, University of Oregon - Eugene, Oreg. 


(ricevuto il 19 Agosto 1958) 


Summary. An estimate is made of the contribution of non-radiative 
inelastic processes to the elastic scattering of high-energy electrons from 4He. 
It is shown that at the presently available energies and scattering angles 
the dominant process is electrodisintegration, and that this and other 
effects may be expected to give a correction to the differential elastic 
scattering cross-section of several percent, the maximum values occurring 
at high-energies and backward scattering angles. 


1. — Introduction. 


In view of the continuing experiments (') on the scattering of high-energy 
électrons from ‘He, it is of interest to extend earlier estimates (?) of the in- 
elastic contributions to elastic scattering from light nuclei to include helium. 
For this purpose, it is convenient to employ the phenomenological treatment 


used in I. 


(*) A preliminary report of this work was given at the Ithaca, N.Y. meeting of 
the American Physical Society, H. S. VALK and P. O. Davey: Bull. Am. Phys. Soc., 


II, 3, 268 (1958). 
(+) Present address: Physics Department, University of Nevada, Las Vegas, Ne- 


vada. 
(1) G. R. BureLson and R. Horsrapter. We wish to thank Dr. HOFSTADTER 


for informing us of the progress of this work. 
(2) H. S. VaLK: Nuovo Cimento, 6, 173 (1957), referred to as I. 


51 - Il Nuovo Cimento. 
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It has been shown in I that the differential elastic scattering: cross-section, 
; including inelastic corrections, can be written in first Born approximation as 


i (1) do,/dQ2 = (do,/dQ2), (1 a C) ’ 
where do,/dQ represents the observed cross-section, (do,/d2), is the .un- 
. corrected cross-section, and C is the inelastic correction term. given by 


o? (k/4x)? | g(q) |? cos? (7/2) 


7) = G5,/d2 — 02k /4zx)? | g(q) | 608? (8/2) * 


Here; o,,.is the total non-radiative inelastic cross-section (which for incoherent 

processes is the sum of the individual cross-sections for these processes), k is 
the electron reduced wave number, # is the angle of scattering, and g(q) is 
an «inelastic form factor » (a function of the momentum transfer qf} taken 
to be of the form corresponding to an exponential type interaction, 


(3) g(q) = (1 + a#q?)?, 


where a is a parameter characteristic of the range of the particular inelastic 
process being considered. 

Since the helium nucleus possesses only one bound state, the principal 
: contributions to o,, will come from excitation to the continum : (electrodis- 
integration) and pion production. As will be shown later, however,.the former 


will make the dominant contribution for incident electron energies below 
1600 MeV. 


‘2. — Inelastic ccontributions. 


.In..order: to calculate the correction C in Eq. (1) it is necessary to deter- 
mine. the total inelastic cross-section, o,,, arising from both electrodisintegration 
. and pion production. 

The main contribution to the former will come from reactions of the type: 
= 4 = 3 
aN e+ *He > e-+*He +n, 
e-+ ‘He > e-+*H +p. 


As pointed out in a recent paper by Muto and SEBE (*), any other modes of 
disintegration may be expected to be of much lesser importance (~ 10%) 
and may be safely neglected in our calculations. 


The differential inelastic scattering cross-section for the reactions (A) has 


(*) T. Muro and T. SeBE: Progr. Theor. Phys. (Japan), 18, 621 (1957). 


» 


oil 


bel 


A NOTE ON THE INELASTIC CONTRIBUTIONS TO THE ELASTIC SCATTERING ETC. 791 


recently been derived (?) using gaussian wave functions containing empirical 
Spread parameters determined by comparison with the earlier Stanford exper- 
imental data (4). It should be pointed out that the values selected for these 
spread parameters are such as to give an rms radius for the «-particle of 
1.69 fermi rather than the 1.41 fermi appropriate to point nucleons (#5). In 
this manner the finite nucleon size effects are roughly included without intro- 
ducing any explicit form factor for the proton (5). 

Thus the contribution to o,, from (A) can be obtained by a direct numerical 
integration of Muto and Sebe’s expression over final electron energies and 
angles. Since, however, we are concerned with the ultrarelativistic case, we 
can employ their simplified Eq. (39) rather than the more exact Eq. (32). 
Unfortunately, all of these cross-sections diverge as the scattering angle ap- 
proaches zero, so that a logarithmic extrapolation has to be used in the angular 
integrations for Ÿ < 40°. Since a 
sizeable fraction of the electrodisin- 
tegration contribution, oî,, to 0, 
comes from this small angle region, 
the results obtained in such a manner 
can only be expected to provide, at 
best, rough estimates. 

By using the above procedure for 
computing of, at 400 and 600 MeV, 
we find the values 51 ub, and 34 ub, 
respectively. Itis encouraging to note 
that these numbers are approxima- 
tely twice those previously reported 
for ?H in I, corresponding to the fact 
that we are considering a system 
having essentially two modes of dis- 10% 
integration instead of one. At such 
energies where the binding effects are 


1079) E =150 MeV 


6= 60° 


10742 


T_T Tt 


d°o/4k; dQ in cmsr 


10°43 


quite small, one may well expect the 
cross-section for one of the processes 
(A) to be about the same as for 
electrodisintegration of the deuteron. 


E in MeV 
20 40 60 80 100 120 


10°45 


Fig. 1. — Differential inelastic cross-section 
at 150 MeV and 60°. 


(4) R. McArristER and R. Horsraprer: Phys. Rev., 102, 851 (1956) and R. Hor- 
STADTER: Rev. Mod. Phys., 28, 214 (1956); Ann. Rev. Nucl. Sci., 7, 231 (1957). 

(5) Actually a better choice of these parameters would be ye! = wr} = 4.3 fermi, 
leading to an rms radius of 1.61 fermi in agreement with the experimental data; how- 
ever, this difference will not affect our results significantly. 

(6) We would like to express our thanks to Dr. T. Muro for an interesting letter 


on this phase of his work with Dr. SERBE. 
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In addition to the preceding calculation, and as a supplement to Muto and 
Sebe’s paper, we have computed representative values of the inelastic. dif- 
ferential cross-section over a selected range of laboratory energies and angles 
so that their Eq. (39) can be more easily compared with experimental data 
when they become available. These results are shown in Figs. 1-5. All of 


E. = 400 MeV 
LO %g=120° 
E,=150 MeV 10°44 
10743 @=120° 3 i = 
sd [ e | 
a WAL, F E 1 | 
"E pa: = 
10744 = = 
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D 
= 104% < 
S FAS 
a i 
© F 
10745 [ 
Re 
L is 
E, in MeV 
1049 MR I ira 
20740 2) 00 2m 80 100 29 50 100 150 200 250 300 350 
Fig. 2. — Differential inelastic cross- Fig. 3. — Differential inelastic cross- 
section at 150 MeV and 120°. section at 400 MeV and 120°. 


them demonstrate the expected kinematical behavior with regard to energy 
and scattering angle: 1) The more the energy of the incident electron, the 
more probable it is that the electron will have to lose a greater amount of 
energy in order to be scattered into the same angle. 2) For a given incident 
energy, the greater the scattering angle, the more probable it is that the 
electron will lose a greater amount of energy. 

As in I, some estimate of the pion production contribution to &. can be 
obtained by considering the reactions (here we again neglect multiple pro- 
duction as a much smaller effect): 


SP seen bm cbt ay > a 
(B) | Bien Sie aaa 


oD eee per e-+n—>e-in-+7° 
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to arise virtually from the corresponding photo-production processes. The 
latter have studied extensively over the entire range of nuclei (**) and the 
results would indicate that for light nuclei: 


(5) ANZI i of (Z, N) PZ). 
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Ee = 600 MeV 
6=60° 
10743 

G 
5 104L 
ae LE 
E o 
[=] = 

10744 = 
3 i 
n°) sv 
3 è 
HS 10 ES 

ad 

1045. © 


E in MeV E, in MeV 
1) pay SERE ig 2 D 300 
200 300 400” 500 i 150 200 250 300 350 
Fig. 4. — Differential inelastic cross- Fig. 5. — Differential inelastic cross- 
section at 600 MeV and 60°. a section at 600 MeV and 120°. 


Since ‘He falls into this category we will assume that its photo-pion cross- 
sections are given by Eq. (5) with the additional supposition that of (n) ~ 
~ ot (p). We can then write in the Weizsäcker-Williams approximation (°) 
the same relations between the electro-pion cross-sections and make use of 
the values found in I. This procedure gives the results shown in Cols. 2-4 
of Table I. It may be expected that these values represent overestimates in 
so far as pion production is concerned because of the Pauli inhibition against 


J) W: Wiccrams, K. Crowe and R. FRIEDMAN: Phys. Rev., 105, 1840 (1957) and 
Ww. Imuor, E. KNAPP, H. EASTERDAY and V. PEREZ- MenDEZ: Phys. Rev., 108, 1040 
(1957). mad 

(2) B. Govorkov, V. Gotpanskis, 0. KARPUHIN, A. KUTSENKO and V. PAVLOV- 
SKAIA: Sov. Phys. Dokl., 2, 4 (1957). i 
' (®) This method is discussed in R. Datirz and D. YEeNNIE: Phys. Rev., 105, 
1598 (1957). i 
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producing charged pions; however, this is more than offset by our neglect of 
the other electrodisintegration processes. The increasing relative importance 
of pionic effects can be seen from Table I: At 400 MeV, the pion contribution 
is only 20% of the total, while at 600 MeV it has risen to 38%. 


TaBLE I. — Values of the total cross-sections. 
En Orn oT sa oF Fin 
400 MeV 51 wb 6 ub 4 ub 14 ub | 65 ub 
600 MeV 34 ub 9 ub 6 ub 21 ub | 55 ub | 


As a very rough check on the order of magnitude of the preceding results 
and the use of the logarithmic extrapolation procedure, we employed the latter 
to find the elastic differential scattering cross-section at zero angle along with 
the total elastic cross-section, o,: These results were then combined with 
those in Table I and found to satisfy the optical theorem inequality (1): 


(6) (do,/4.Q)y_ 9. > Ou + 0,,)2/1622 . 


3. — Calculation of the inelastic correction and discussion. 


The calculation of C according to Eq. (2) requires, in addition to o,,, some 
knowledge of the appropriate inelastic form factor g(q). This function has 
a Shape and range dependent upon the nature of the inelastic processes being 
considered. Since the only conditions placed on it mathematically are that 
g(0)=1 and g(q)<1 for |q|<za-!, there are many admissible functions 
that can be chosen, each corresponding to a different form of the additional 
term in the Hamiltonian. Of these possibilities, perhaps the simplest is the 
Yukawa type chosen in I; this function, however, corresponds to an interaction 
possessing a non-physical singularity at the origin, thus placing undue em- 
phasis on the region of high momentum transfer, and leading to unreasonably 
large values of the inelastic correction. The invalidity of such a model is 
verified by comparing the results of Eq. (1) (using a =1.61 fermi) with the 
experimental data at 400 MeV: Eq. (1) predicts too rapid a change in the total 
form factor for increased angles of scattering. Physically, one would expect 
the strength of the inelastic term to increase with decreasing radius r, until 
y =a, after which it would level off or increase to a finite maximum in the 


(!°) See, for example, L. I. Lapipus: Journ. Exp. Theor. Phys. (Sov. Phys.), 4, 
937 (1957). 
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manner: of a gaussian or exponential function. Here a would represent some 
characteristic distance of the order of-nuclear dimensions. 

For the present estimate we have chosen -g(q):to -be of-the form of-an 
exponential, Eq. (3), and taken a — 1.61 fermi, the rms radius of the «-par- 
ticle. The only remaining factor is that for the elastie cross-section (do,/dQ). .- 
According to the Stanford data, this. may.,be. taken as 


(7) (do,/49) = [e*/(k#c)]?[cos? (9/2) /sin* (8/2)]|F(q) À, 
where 
(8) F(q) = exp [— q?a?/6] 


is the form: factor and a has its previous meaning. 

Using Eqs. (2), (3), (7), and (8) and the values from Table I, we have. 
computed C as a function of scattering angle in the electron-‘He barycentric« 
system. The results are shown in Table II. 


TABLE II. — Values of the inelastic correction C:° 


aa 40° 60° 80° 100° | 
| 


400 MeV 4.6-10-9 7.02 107% 2.1 1055 
600 MeV | 5.0-10-9 3.5-10-5 7.110 


As in the earlier cases of the proton: and deuteron we see that the cor: - 
rection term increases with increasing momentum transfer. Especially note- - 
worthy is the sharp rise with increasing angle, characteristic of the corres- 
ponding decrease in the elastic cross-section. This, of course, is due to the 
rapid decrease of the form factor, Eq. (8); over the same range. In all cases 
it can be seem that the correction is not greater than 10%. ,It should be noted, 
however, that although the correction term is small at the present exper- 
imentally available energies and angles, as represented in Table II, it tends 
to become much larger at increased backward angles. We have not considered 
this region in the present paper since the uncertainties in the inelastic form 
factor at high momentum transfer and the limitations on the present: theory 
imposed by the use of first Born approximation do not warrant it. Such-large 
values of the correction: must be taken only as indicative of the inereasing 
importance of inelastic effects at high momentum transfer and no more. Indeed, . 
more knowledge of the inelastic form factor g(q),; for large q is required before 
further quantitative results can be given. It is clear, however, that if the - 
strength of the electron-nuelear:inelastie term levels off: appreciably for r.< @:1 
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g(q) will decrease more rapidly for increasing q than (1+ q?a?)*. In such 
a case the correction can be expected to remain of the order of magnitude of 
the values listed in Table II over the entire range of momentum transfer. 


The authors would like to thank Professors H. PRIMAKOFF and H. EASTER- 
DAY for enlightening discussions on various aspects of this work. In addition 
we would like to express our appreciation to Mr. A. JOHNSON for supplying 
us with his helium wave function calculations before publication. 


RIASSUNTO (*) 


Si valuta il contributo dei processi elastici non radiativi allo scattering elastico 
degli elettroni di alta energia derivanti da ‘He. Si dimostra che alle energie e sotto 
gli angoli di scattering attualmente disponibili il processo dominante è l’elettrodisin- 
tegrazione e che presumibilmente questo o altri effetti debbono apportare una corre- 
zione di alcuni percento alla sezione d’urto elastica differenziale, presentandosi i valori 
massimi alle alte energie e sotto angoli di scattering retrogradi. 


(*) Traduzione a cura della Redazione. 


Ling 
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On the Relation between Velocity and Mass of the Electron (*). 


S. RABOY and C. C. TRAIL 


Argonne National Laboratory, Lemont - Ill. 


(ricevuto il 30 agosto 1958) 


. Summary. — It is suggested that .the elastie scattering of electrons, from 
electrons be used to verify the Lorentz law of mass variation with velocity. 
The kinematics of such a process yield angular dependences which are 
sensitive to the model used. An ‘independent measurement of the velo- 
city of the incident particle is necessary to complete the measurement 
in a. logical fashion. It is shown that two other experiments reported in 
the literature are deficient. in .this respect. à 


1. — Introduction. 


In a recent article, FARAGO and JANOSSY (1) reviewed the experimental 
evidence supporting the Lorentz formula relating the mass of a particle to 
its velocity. They conclude that the Lorentz formula has less experimental 
support than is desired for such a fundamental relation. While the exper- 
iments indicate that mass does vary with velocity, a distinction cannot be made 
between the: Lorentz.formula and the results from other theories such as that 
of .ABRAHAM (?). _ Furthermore, the experiments actually determine the va- 
riation of e/m with velocity. In view of the universal use of a special functional 
form relating the mass to the velocity, it is desirable to design an experiment 
which will achieve that specific. measurement. 


(*) Work performed under the auspices of the U. S. Atomic Energy Commission. 

(1) P. S. FaraGé and JAnossy: Nuovo Cimento, 5, 1411 (1957). : 

(2) M. ABRAHAM: Ann. der Phys., 10,.105 (1902). We have translated this article 
into English and reprints are available upon request. aio . 
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In 1932, CHAMPION (*) reported a cloud chamber measurement which pur- 
ported to confirm the Lorentz law of mass variation quantitatively. His error 
assignments were criticized by FARAGO and JANossy, but were defended by 
CHAMPION (4). We would like to add some additional comments about the 
measurement and to comment on another paper which describes a similar 
measurement. We will attempt to demonstrate that it is impossible in prin- 
ciple for these measurements to provide unambiguous verification of the Lo- 
rentz formula. 


2. — Discussion. 


We suggest an experiment which will indeed measure the variation of the 
mass with velocity. Such an experiment involves scattering electrons from 
electrons, measuring the angle between the scattered electron and the recoil 
electron, and measuring the velocity of the incident electrons, perhaps with 
a time-of-flight technique. Fig. 1 is a sche- 
matic diagram of the proposed experiment. 

Application of the principles of conser- 
vation of momentum and energy yields the 
following collision equations. 


(1) Pi = P2 COS P + Ps COS (0 — gp), 
Fig. 1. — Schematic diagram of ; è 
the electron-electron scattering (2) Da ing = paia pla 


experiment. (3) PILE = LE. 


where p is the magnitude of the momentum and £ is the energy, H, is the rest 
energy g is the angle one of the recoil particles makes with the incident di- 
rection and 0 is the angle between the paths of the recoil particles. The sub- 
script 1 designates the incident particle, 2 denotes the scattered particle, and 
3 the recoiling particle. All quantities are measured in the laboratory co- 
ordinate system. Classically, 0 is always 90° since the particles have the same 
mass. 

To obtain a solution of these equations one specifies that momentum and 


energy are related by an equation which depends upon the model. Classically 
this relation is 


(4) p? 


I 

to 
= 
I 


(*) F. C. Campion: Proc. Roy. Soc., A 136, 630 (1932). 
(4) F. C. CHamPion: Nuovo Cimento, 7, 122 (1958). 
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The Lorentz model gives 
5 fa 202 
(5) E? = p'c? + mes, 


and in this theory 


(6) DAI I 
RA 

- m 2 

(7) E = af} ; 


where f is the ratio of the speed of the particle to the speed of light. The 
expressions for the energy and momentum based upon the Abraham model (5) 
may be written 


eee (eer, Lee 2 
È) Ex 8 ae | eas ae 23 | k 

au a [1 is os ee i 
(9) E=] et [prog 5 i ! 


/ 
LORENZO / 


Fig. 2 shows the relation between p and f for the 
latter two models. A similar curve is obtained for 
the energy relations. Fig. 3 shows how the energy 
varies with momentum. These curves illustrate the 
similarity between the models and, more importantly, 
show that those experiments which depend exclusively 
upon the relationship between energy and momentum  Fig.2.- Relation between 
are not very sensitive to the model used. Indeed in momentum and velocity 
the extreme relativistic case (~~ 1), E=pe for both for the Abraham and 
models and the solution to the collision equations (1), Lorentz models. 

(2) and (3) are identical for both models. 

For purposes of designing an experiment (‘), we present in Figs. 4, 5 and 6 
relationships between 0, g, and f,. In Figs. 4 and 5 the difference between 
9 and 90° is plotted against œ for different values of 6, for the two models. 
Fig. 6 gives the relationship between (90° — 6) and f, for p= 40°. These 


(5) In Abraham’s expression for the rigid sphere with uniform surface charge, we 
have identified the coefficient e/62ac? with the rest mass where e is the charge and a 
is the radius of the electron. This identification gives proper classical limits for the 
kinetic energy and momentum but the rest energy becomes (3/4) me?. However, in the 
collision equations these coefficients drop out. 

(6) M. B. Scorr, A. O, Hanson and E. M. Lyman: Phys. Rev., 84, 638 (1951). 
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quantities have been calculated with the digital computer GEORGE at Argonne 
National Laboratory. From Figs. 4 and 5 it is seen that the two theories give 
quite different results for the 

100 same B,, but the shapes of 
j the curves are similar and, 
for example, the curve for 

Ei B:= 0.99 in the Abraham 


ve ate | incides witl 
È ge e theory nearly coincides with 

Ge a A that for B, — 0.95 in the Lo- 

E 7 (CALCULATED POINTS) ILA 

AA rentz theory. This illustrates 


the need for an independent, 
È measurement of the velocity. 
[ j This type of experiment 
ae fa tie SI RATE has three virtues: 1) No 

assumption need be made 
concerning the relationship 
between charge and velo- 
city. All other experiments 
of which we are aware involve the measurement of the ratio of charge to 
mass and the assumption is made that the charge is independent of velo- 
city. 2) No law of force bet- rift 


È pr) 


Fig. 3. — Relation between momentum and energy 
for the Abraham and Lorentz models. 


ween the particles is assumed. ‘50 

Only conservation of momen-, | de 
tum and energy is required. : op a 

3) The coefficients.. involving E) SIOE 

the rest energy cancel, remov- : : ole - 


ing an ambiguity.in the Abra- 
ham model. 

An uncertainty in the ex- 
periment is introduced by the. . 
motion of the electrons in the 
atoms of the target. ForD 
and MULLIN -(7) have shown 
that in the. ease: of a beryllium 
target the spread in 0 is about 
3°, but this spread decreases 


0 c D IN DEGREES . = 
10 20 do 70 “50 60 70 80 


Fig. 4. — Variation of (90° — 6) with @ for the 
Abraham model. À 


as Z?. Our calculations (5) show that it is necessary to obtain the velocity to 


0.1% and the angles to 0.1° if one is to es a parigeulae uate ho OLO 


(7) Gs W. ForD A C.J. Minnie Phys Rév., a, 520 SES), ui o com 


munication. tH) 
(5) Calculations on the Abraham model vo been sede ft Di na 


0.99 and 
all permissible laboratory angles. These results are available upon request. 
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Such measurements seem just possible if a Van de Graaf generator is used as 
a source of electrons, and movable scintillation counters in a scattering cham- 
ber are used to measure 0 and g. As mentioned above, the incident bi must 
be determined with a time-of-flight measurement. It appears that an optimum 
6, is about 0.9. 


“ot LORENTZ 
| 
| 


B=099 


20 B=095 
| =090 
wi 
| ® IN DEGREES ‘ 
Oo Hama 20 so 60 70 80 
Fig. 5. — Variation of (90° — 0) with © Fig. 6. — Variation of (90° — 0) with f 
for the Lorentz model. for pg = 40°. 


We discuss now the results of two experiments which we will call the 
Champion experiment (*) and the Illinois experiment (°). 

Equations (1), (2) and (3) reduce, upon assumption of the Lorentz model, 
to the equation 


(yi — 1) sin peospo 
(y +1)? sin? y + 4 cos? y}? 


cos 0 =- 
{ 


where y = (1—f?) 74. 

In 1932 CHAMPION, who used radium Æ as the source of electrons, reported 
some results from electron-electron scattering in a cloud chamber. His measu- 
rements of g and 0 for various incident velocities are given in Table I. The 
subscript (obs) is used to indicate the observed angle 0. The incident velo- 
city f, in Table I was determined by measuring the curvature of the path of 
the particle in a magnetic field. This measurement gives momentum (actually 
Heo) and one must assume a model to obtain the velocity. Apparently CHAMPION 
assumes the Lorentz model and obtains a consistency check as indicated in 
Table I. We have calculated the corresponding quantities on the Abraham 
model and find that the angle @ calculated from this model agrees with Cham- 
pion’s experimental results. We do not agree with the footnote on p. 634 
in CHAMPION’s article. He reports 0,=58° for 6, = 0.91. We find 0, = 76.3° 
for this collision. 
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TABLE I. — Champion’s (3) measured angle between the recoil and the scattered particles 
compared with the angle calculated with the Lorentz model and with the Abraham model. 
Here f, and fy are the velocities of the incident electron (in units of the velocity 
of light) as calculated from the measured momentum from the Lorentz and Abraham 
model, respectively. The uncertainties in 8, and f, result fromthe uncertainties in the 
measurement of the momentum. Angles y and 6,,, are as defined in the text and their 
accompanying uncertainties result from the errors of measurement. Angles 0, and 64 
are calculated from the respective models with resulting uncertainties caused by errors 
in the corresponding f and 9. 


p+0.5 0 6, 0, 

Pa vA (°) (°) (°) (°) 
0.85 +0.01 0.90 +0.01 20.0 83.6 + 1.0 82.7 +1.0 82.9 +1.0 
0.83 +0.01 0.88 +0.01 26.6 81.2+1.0 81.7 +1.0 82.3 +1.0 
0.83 +0.01 0.88 +0.01 31.4 81.0 +0.1 81.0+1.0 81.8+1.0 
0.82 -+0.01 0.86 +0.01 22.0 84.1+1.0 83.2+1.0 83.9 +1.0 
0.85+0.01 0.90+0.01 wae 82.2+1.0 81.7+1.0 82.2+1.0 
0.83 +0.01 0.88+0.01 DIA 82.1+1.0 82.3 +1.0 83.1+1.0 
0.84+0.01 0.89+0.01 23.4 82.7+1.0 81.7+1.0 82.4+1.0 
0.90 +0.02 0.95 £0.02 24.5 79.6+1.0 77.4+1.5 78.2+1.5 
0.88 -+0.01 0.93 +0.01 35.4 76.8 +0.5 77.3 +1.0 78.3 +1.0 
0.85 +0.01 0.90+0.01 21.1 82.7 +1.0 81.8 +1.0 82.6 + 1.0 
0.91 +0.02 0.95 +0.02 36.9 75.2 +0.5 75.2 +1.5 76.3 +1.5 
0.93 +0.02 0.97 +0.02 29.6 72.5 +0.5 72.6 +1.5 74.1 +1.5 
0.85 +0.01 0.90 +0.01 21.8 82.4 +1.0 81.4+1.0 82.4+1.0 
0.82+0.01 0.86+0.01 36.9 80.6 +0.5 81.0+1.0 82.2+1.0 
0.87 +0.02 0.92 +0.02 24.7 83.6+1.0 79.6+2.0 80.6 +2.0 


These results lead us to conclude that the Champion experiment gives no 
evidence to support the Lorentz model. If, however, a direct measurement 
of the incident velocity (or of the velocity of one of the recoil particles) is made, 
then a check of the validity of the models can be made. This measurement, 
of velocity must be made independent of momentum-energy relations. 

We close with a brief discussion of the Illinois experiment. That experi- 
ment was similar to Champion’s except that a magnetic analyzer was used 
to measure p, as a function of y for a fixed p,. (The Lorentz energy of the 
incident particle was 15.7 MeV). However, to compare the data with theory 
it was necessary to obtain a relation between momentum and energy. We 
have seen above that the models are not sensitive to these relations and in fact 
become identical for 6, close to 1. We conclude that this paper offers no evi- 
dence for believing that the Lorentz model describes the electron more ac- 
curately than the Abraham model. It should be noted that the Illinois group 
does not claim a decisive experiment in this respect. 


» 
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* * * 


We wish to express our thanks to Dr. DAVID PHILLIPS of the Argonne 
Applied Mathematics Division for making the calculations. 


RIASSUNTO (*) 


Si suggerisce di usare per la verifica della legge di Lorentz sulla variazione della 
massa con la velocità lo scattering elastico degli elettroni su elettroni. La cinematica 
di tale processo fornisce dipendenze angolari soggette all'influenza del modello usato. 
È necessario eseguire indipendentemente la misura della velocità della particella inci- 
dente per completare la misura in modo logico. Si dimostra che a tale riguardo altri 
due esperimenti riportati in letteratura sono deficienti. 


(*) Traduzione a cura della Redazione. 


IL NUOVO CIMENTO VOLS ks Ne aD 1° Dicembre 1958 


Space-Charge Effects in Electron-Synchrotrons. 


C. BERNARDINI 


Istituto Nazionale di Fisica Nucleare 
Laboratori Nazionali del Sincrotrone di Frascati 


(ricevuto il 16 Settembre 1958) 


Summary. Space charge effects are considered (taking into account 
positive ions) with the aim of deciding whether they can set up serious 
intensity limitations for the beams of weak focusing electron synchrotrons, 
or not. The possibility of linear resonances and non linear beam splitting 
is shown. An estimate of the rate of production of positive ions in the 
doughnut is given. 


1. — Introduction. 


Space charge effects in electron-synchrotron are usually settled on by the 
statement that, if the injection energy is high enough, magnetic forces com- 
pensate nearly at all electrostatic forces (the well known 1 — f? factor (1)). 
This statement is, of course, valid in view of the till now obtained intensities 
for such machines, i.e. having in mind a figure of 101% for the order of mag- 
nitude of the number of circulating electrons per pulse. 

It is still assumed, even if not explicitly said, that the electron beam moves 
in vacuum but for the problem of scattering losses. Consideration of scat- 
tering losses gives an intensity independent rule for a troubleless pressure of 
the residual gas; now, we want to show that this rule does not eliminate the 
possibility of ion production, in the beam-occupied region, in such a quan- 
tity as to give rise to troubles with the focusing properties of the external 
magnetic field plus space charge forces. 


(1) See for instance: W. HARMAN, Fundamentals of Electronic Motion (New York, 
1953), pag. 291. . 
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A complete formulation of the problem of positive ions production and of 
their interaction with the electron beam is complicated by the number of 
concomitant circumstances; nevertheless, orders of magnitude and semiquali- 
tative results can be obtained by the use of simple models. 

Two main results will be shown in the following paragraphs: 


a) The possibility of reaching dangerous resonant values of the betatron 
frequencies. 


b) The possibility of beam-splitting into two or more parallel beams. 


It is hard to maintain that the strange facts observed in some electron- 
synchrotron laboratories (?) find an explanation in the results of this paper. 
However, we think that such results can give a guide for an easy experimental 
program with the aim of deciding if positive ion effects are to be accounted 
for (or rejected) in planning high intensity machines. 

Eventually, we anticipate that typical appearance values of the effects in 
weak focusing machines (which we are interested in) are 5-10-*mm Hg for 
the pressure, 5-10!° for the number of electrons per pulse. The main pertur- 
bation parameters depend on the product pressure X intensity. 


2. — Resonance effects. 


In this paragraph we want to study the effect of an uniformly distributed 
sea of positive ions, charge +e, on the betatron oscillation of the electrons 
in the beam. We think of the ion density as constant. 

Let us call N, the total number of circulating electrons, N, the number 
of ions in the beam occupied region and put 


eu 


NA 


Suppose for semplicity the beam is made up of four bunches, uniformly 
filled, with electrons and cylindrical in shape, the length, along the axis (pa- 
rallel to the direction of the bunch), being RA®, the circular section having 
a radius »,. The total radial force on an electron of the beam, at a distance r 


from the beam axis will be 


e? 
F, = à N, 


a 2 | 
9 RAGE" “sin 


(2) CERN symposium: 1, 67, 301 (1956). 


52 - Il Nuovo Cimento. 
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The 1— f? dependent part is due to electron-electron forces, the e part 
is due to electron-ion forces. 
We neglect the 1—? part, thus assuming 
Lis Pie 


the equations of betatron oscillations take the form 


e+ (1—n-+ Ôn)x = 0 (radial motion) 
"+ (n + dn) = 0 (vertical motion), 
where 
i) eR 
Da oe = meA@r’ 


and the other symbols need no explanation. 

The dependence on electron energy of this space charge parameter $n is 
not easy to guarantee: one could say that the quantity mA@r behaves 
roughly as B * (B being the main magnetic field in the machine) due to adia- 
batic dampings. But this energy dependence could be wrong because of a 
number of influencing circumstances (like scattering, radiation fluctuations or 
space charge effects themselves). We take a constant dn value, corresponding 
to injection values of the parameters; «injection values » means, of course, 
values reached soon after RF capture, when the electron beam density is nearly 
steady. 

For typical values as 


R = 360 cm, (met). = 2-5 MeV, 
FL CMS AD = 1 rad, 
one obtains | 
N, 
on = 0.le (aaa 7 


We can write down a resonance relation for the betatron frequencies, for 
a racetrack synchrotron with four straight sections of length L, in the form: 


ie piercer cao LE “ete 
p+ (1+ piatt va ton) =0, 


where p, q, r are integers (2 0) (°). 


() E. PERSICO: report T1 della Sez. Acc. dell’ INFN (1953). 
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In the case of L/xR = 0.106, n = 0.61 (as for the Frascati machine), there 
are three linear resonances near the unperturbed working point, namely 


pel SV) r=—-1 dn 0.2 
p=—2 g= 1 r=1 On ~ 0.3 
rh q=—1 r=) one 0.4 


The first resonance (dn = 0.2) reveals itself by the build up of vertical 
oscillations due to first harmonic median plane irregularities. The second one 
(dn = 0.3) is a coupling resonance excited by field or median plane second 
harmonic irregularities. The last resonance (dn — 0.4) is a field first harmonic 
excited radial motion (*). 

Non-linear resonances are not considered here because of the surely longer 
amplitude build-up time, which could allow for ion density rearrangements 
destroying the resonant value of dn. 


3. — The ion density. 


We want to estimate the quantity dn as a function of pressure, beam in- 
tensity, atomic number Z of the residual gas, etc. We neglect the secondary 
electron production since, as can be easily shown, secondary electron velocities 
are > than ion velocities, so that they rapidly leave the beam region. 

Let us call N the number of the atoms per cm’. At 300 °K of temperature: 


N = 3.2-101P (em), 


where P is the pressure in mm Hg. 
c; is the ionization cross-section: it is energy independent for relativistic 
electrons (4). We shall take as in reference (*) 


0,7 0.2Z-10-18 cm’. 


Put t, for the time an ion needs to escape from the beam-occupied region. 


(*) Resonant blow up of the beam, of the just considered kind, seems likely to be 
observable only in the vertical direction because of the need of having ions over the 
whole path of the electrons during the oscillation amplitude build up time. This has 
been kindly pointed out to me by doctor A. TOLLESTRUP. 

(4) M. J. Moravesix: Phys. Rev. 100, 1009 (1955). 
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It follows that 


t 
Ne cow [NE exp[— t/t] dt’ . 


We suppose N, is made up of two parts: 
N.= N,+N,exp[—t/r], 


where N, and N, are constant; t is a characteristic time of slow losses. 
When t>t,, T, e reaches a steady value 


Fix via 
© = CotyoN (1 rh oad ) 


t+t, Na 
the factor 1+ (t/(t-+t,,))(N,/N.) is likely to be about two, so that 
e 26ctno; N . 


The main problem is to calculate t,: The ions, since the time of their pro- 
duction, move along the lines of force of the external magnetic field, that is 
along the 2 axis. Thus, for an ion moving with velocity v, and produced at 
a random point of the beam section, 


8To 


ire | 
“  8x|v| 


What we need is the average «1/4,» over the possible |v| values, |v| being 


_ the modulus of the 2 component of the ion velocity (see Appendix A). One 


could take the average on thermal velocities (*), thus obtaining 


al 3 (2rkT\* 
ty 27 M È 


(*) Recoil energies of the ions do not change appreciably the termal energy distri- 
bution. The transferred momenta q are certainly less than 10 keV/c in a collision, so 
that the ion mean recoil energy is certainly less than 


ha 108 0,025 
SM + “471080 (E 


eV , 


which is less than thérmal energies by a factor 1/Z. 
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where M is the mass of the ions; but ion-ion and electron-ion electrostatic 
effects are certainly not negligible from the point of view of the ion distri- 
bution. There is the possibility of ion trapping by the main beam, which we 
shall examine in Appendix B. For the moment being, assume that the thermal 
distribution of velocities is valid; it follows that 


ty 37 V2 us (ro in em) 


assuming M ~ 2Z (mass of the H atom). 
The formula giving € is 


e= 1200 fr, Z°°P, 


where P is in mm Hg. The formula for dn is thus 


on =~ 120 


CZAP PN, fr, in em 
107 |P in mm Hg. 


0 
We see that for 
Ces 210,  P—6:10 4,  7,= lem, Md 10 


one obtains 


ne 0.1. 


4, — Beam splitting. 


So far we were concerned with a simplified model having cylindrical sym- 
metry. Beam splitting effects can be put in evidence by simply removing this 


assumption. 
We will treat in Appendix © the general case of ion distribution; for the 


sake of clearness we now consider a simple model in wich the ion distribution 
has elliptic asymmetry: 


Ne 1 + 26, Sin? p 
n+(r, p)rdrdg = Ask AD face, r dr dy Tr <T 


0 rst 


where n, stays for the ion density as a function of polar co-ordinates around 
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the electron beam axis, r and g. & is a constant measuring the ellipticity, 
satisfying e >—} because cf n, > 0. SN, 

The electric field components due to such a charge distribution are, at 
points far from the head or tail of a bunch 


Dee r\ 02 — 2 , © 1 I a eZ 
B, = o(1 (+) ex 40 gt; x 


2 
RAR To ve 1+64 Toft 


‘ Pon Ai e 1 r\ xz 
B,=0(1 ety (È -inz)* = }+ +40 (+ | =X, 


1 +e, \4 To Te 1+ @ Tibet 
where 
olio 0h z2=arteq, 
sen 
_ 2r:RAD' 


Limiting the analysis of possible motions to a simple case, i.e. to the motion 
on the median plane 2 = 0, we see that the equations of betatron oscillations 
for the electrons of the main beam 

are of the following kind: 


Pot. energy 


à | 
1) è + (0; + ew? In EE ip 


0 

The log term is a new feature for 

betatron oscillations; &,, m,, are two 

constants readily obtainable in terms 

of the parameters defining the fields 

» Ex, E, and the external magnetic 
field. 

The motion described by eq. (1) 

can be visualized by means of the 


graph (Fig. 1) of the potential energy 
of the force 


(0; + €,&i In =) æ, 


a 


for the case e > 0. (Were actually e, <0, the same conclusions hold true 


for the motions out of the median plane, round x — 0). 
Two equilibrium positions are found for the beam; the axis (7 = 0) is no 
more a stable equilibrium position. 
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The parameters defining the well-depth and well centers are quite sensi- 
tive to the ellipticity e,; in the case shown here, for instance, 


| 
|| = exp 


2 

re 2, 
2)? 

CAO 


where +, are the well-center co-ordinates (*). 


5. — Conclusion. 


This report summarizes a list of possibilities for space charge effects due 
to residual gas in the doughnut. 

We do not maintain these effects will be so strongly evident in a well 
working machine as can be conjectured from some of the reported formulas: 
these are, in fact, rough calculations containing too many unknowns. 

Nevertheless we are afraid that intensity limitations for weak focusing 
machines follow on this reasoning line; the point is whether they are so near 
the actual working points as here shown or not and this, we think, can be 
only decided by a suitable experimental program. 


xe * 


We would like to thank Professor G. SALVINI for stimulating help. 


APPENDIX A 


Simple Boltzmann equation for the ion density. 


Let us call n,(a, 2, v, t) the density of ions in a space of coordinates #, z, v 
(v being the 2 component of the velocity). Put n,(x, 2, t) for the electron density 
in the beam and f(v) for the maxwellian distribution function of the velocity v. 
Then 
on: On+ 


da v à Lo, Nen,(@, 2, t)f(v). 


Let us consider the case in which 
nel) 2 t) = y(2, 2) N.(t), 


(*) Once two beams are produced further splitting could be envisaged for each 
one of the two beams on the same basic mechanism. Prof. R. Witson pointed out to 


me this possibility. 


co 
= 
= 


C. BERNARDINI 


p(x, +) de de = 1. 


_ whole space 


e; Der A solution is readily found with the only condition that there are no ions 
(«© at the timet=0. 
tom Di n4(0,8,0, t) = co No) y [2, z— v(t'— t)] N.(t) dt . 


0 


This formula allows us to 
perform the correct average 
over the velocities. Namely, 
what we need is © 


T(z, x, t) = rover, z — vt) dv. 


Now, f(v) is a gaussian 
function of width <v*>; y, as a 
function of 2 — vt, is a strongly 
peaked function around the 
value 2 — vt = 0 of the argu- 
ment. The width of this peak 
: is of the order r,. When (Fig. 2a) 


$ dl , Co»? > ri ? 


‘ Fig. 2. 
SR we can substitute for y a 6 function: 
a pease px, 2)-= u(x)d(2 — vit). 


oi 
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Here uw takes into account the x dependence. Thus, for this case 


1 M \ Ma 1} 1 M \ 
I(x e u(x)| —— | exp |— -—| ~ — u(x)| ——] . 
Peas ae wee sean) Fo | eT al H (sa) 


The last approximation comes out from the fact. that we are mainly inter- 
ested in 2 values less than 7. 
On the other hand, if 


Ob = 19, 


and if 2 # 0, the p(z — vt) function, as a function of v, displaces its peak toward 
infinity as t > 0, so that it does not overlap appreciably with the velocity distri- 
bution, (Fig. 25). 

There is an intermediate region in which 


CITA 


where things are somewhat more complicated. This region is, however, of no 
interest for the final result. One can check all the qualitative considerations 
just made by calculating what happens for a gaussian shaped y. Eventually 
one obtains that: 


I(@,z,t) = 0 : CODE TE 


L M 


FL Fa) COLE: ME 


land 


QukT 


Using this result to calculate 


+00 
ny) =| n+d0 , 


—@o 


one obtains 


t 
i PI M \} atl. 
n> = eo. [NI z,t—t') dt’ = co, Nu(x) sat) N,(t —t rr ; 


ro/V<o® 


neglecting the time dependence of N,: 


L M \} Cw 1? 1 M \i 2kTt 
<n+> aed Len Mula) 5-5) N, In (2) = 9 04 TC er N, In Mr e 


0 
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The log term is of no importance. 


. 1 M \} ae 
~ x = ‘£ Il 
Ni EST <n» dx 5 CNT & a) N, n n) 


This calculation shows that the qualitative arguments of Sect. 3 can be 
followed with some confidence. 


APPENDIX B 


Ion trapping. 


The ion trapping mechanism can be described by means of matrix technique 
as follows: the ions in a given azimuthal position are periodically inside a 
bunch of electrons for a time RA@/c, outside of the bunches for a time 
(RJo)((x/2) — AD). 


When inside, the ions equations of motion are of the form 


Z+o?z2 = 0; 


when outside, however, 


where 


This is a time alternating focusing defocusing strueture, whose matrix 
trace is given by 1 


1 
3 Tr (transf. matrix) = 


= ns eli 3 = 
= cosh y,Ve cos y,V1 —e era == sim y, Ve sin y Vite. 


where we put for brevity 


ee CN.Re®? \t/(x/2) — AD (x/2) — AD 
vn) \2Mer Ad AD )=>( AD ) 


When y <1 the analysis of the trace is quite easy: 


®» 


3 Tr (transf. matrix) ~ 1 — $y(y, + pe) +d ef + mn), 
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and there follows that trapping is possible until e reaches the value 


2A 
PC tte, 8 


È) 


IT 


the assumption of this ¢ value as equilibrium value would lead to a pressure 
independent shift of betatron frequencies. But the just described trapping 
mechanism must be interpreted with some care: in fact we claim that the 
situation we shew in Sect. 3, based on thermal motion considerations for free 
ions, is more likely to be true for the following reasons. Let us suppose an 
ion oscillates along the z axis (i.e. it is trapped). If it leaves the beam axis 
with a velocity <v> (thermal average), the amplitude of its oscillations will 
be <v>/Q where £ is the characteristic frequency of the composite structure 
we deal with in this appendix. Trapping is effective only if 


<v) 


Q < T0 » 


because if the ion goes outside the electron beam boundary repulsive space 
charge forces prevail and make it reach the doughnut’s walls. Now 2 is given, 


in the y <1 approximation, by 2 ~ V87y(AD — (7/2)e)*(c/R) and 


Kay » CO» È RAD i 
Q ~ V8ay (AD — (xf2)e)t e °° £(N,/10*°)(A®D — (x/2)e) 4 


and this shows that trapping will be effective in fastening the ion cumulation 
but, as soon as © —0, the trapping mechanism becomes less efficient and 
the pure thermal motion considerations give a good approximation. 

We remark that because of this phenomenon the pressure dependence 
of e could be far from linear. 


APPENDIX C 


The ion distribution. General case. 


To show the possibility of appearance of log terms in the electric field 
formula, we develop here the potential due to an ion distribution of density 
n+(r, gy). The potential V(r, g) is given by 


27 (co) 


V(r,œ) = const 26e ag? drn.(r, g)In Vr? + 72 — 2rr cos (yp — Y) = 


0 0 


= const — 2CeK(r, ®) ’ 
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Expanding in Fourier series 
+0 
nA, p) = #2 qu(r) exp [ike] (g, = 9") 


and 
K(r; 9) =a In r Faw dr + x In rq,(r) dr — 


0 r 
T 


wt?’ exp[ik int = the 1 ar ee 
9 Poi LE + fra dr Lia qr(T) vi ‘ 


rl 
0 


Logarithmic terms can appear in the last term of the expansion, whenever 
g(r) — rl? (for not too large r; k #0) 


The case we studied in Sect. 4 is that of qs = const for r<r,. Attention 
must be given to the condition n}(r, g) > 0; in the model of Sect. 4 this condi- 
tion requires, for instance, ¢ > — }. 


RIASSUNTO 


Si studiano gli effetti di carica spaziale (tenendo conto della presenza di ioni posi- 
tivi) per stabilire se possono dar luogo a serie limitazioni per l’intensità del fascio di 
un elettrosinerotrone a focalizzazione debole. Si mostra che sono possibili effetti lineari 
di risonanza ed effetti non lineari di sdoppiamento del fascio. Si valuta inoltre il numero 
di ioni positivi prodotti nella camera a vuoto. 
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Lagrangian Formalism in Relativistic Hydrodynamics 
of Rotating Fluid Masses. 


F. HALBWACHS, P. HILLION and J.-P. VIGIER 


Institut Henri Poincaré - Paris 


(ricevuto il 17 Settembre 1958) 


Summary. — The general theory of relativistie moving fluid masses is 
studied in terms of a new set of variables originally proposed by Einstein 
and Kramers. The general equations of motion established directly by 
Bohm and one of us (J.-P. V.) then appear as conservation equations 
deduced from a scalar Lagrange function of these new variables. This 
procedure paves the way for a quantization of the internal motion of a 
model of extended particles in real space: model strongly suggested 
by the experimental results of Hofstadter. 


Introduction. 


In a recent paper the general theory of the relativistic. mechanical behaviour 
of rotating fluid masses has been treated by BOHM and one of us (1) and it was 
possible to show that the motion of various significant global physical quan- 
tities (such as the total linear momentum G,, the total angular momentum M,,, 
the center of mass, and the center of matter density) satisfy equations of 
motion which generalize simple relations postulated by WEYSSENHOFF (?). 
The problem we want to treat here consists in the deduction of these equa- 
tions as conservation equations resulting from a scalar Lagrangian depending 
on 2 new set of variables connected. with global properties of the droplet. 


Such a theory is interesting from two main points of view. 


(1) D. Boum and J.-P. VIGIER: Phys. Rev., 109, 1882 (1958). In this paper we have 
used, as far as possible, the same notations. Greek indices such as # vary from one 
+o four; latin indices such as k, from one to three. Indices repeated twice are sum- 
med over all possible values. We work in Minkovski’s space-time. 

(2) J. WEYSSENHOFF: Acta. Phys. Pol., 9, 8 (1947). 
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First it might pave the way to treat in a covariant way the problem of the 
classical and quantal behaviour of particles extended in space. This seems 
absolutely necessary in the present stage of physics for various experimental 
and theoretical reasons. 

As one knows, all recent experiments of high energy scattering between 
various types of elementary particles seem to indicate that such particles are 
not point-particles but have indeed areal extension in space (with a radius of 
the order of 0.6-10-! em). 

This results from many indirect evidences (which we shall not analyze here) 
and also from the remarkable proton-electron scattering experiments of HoF- 
STADTER and his collaborators (*). 

On the other hand many of the divergences. of quantum theory result 
precisely from the point-like character assumed arbitrarily for all particles in 
the present field theory. 

As a consequence it is not unreasonable to hope to eliminate them by 
starting a priori from a classical extended model. 

Secondly it provides a convenient basis to study internal stable excited 
quantum states of such droplets; thus generalizing the classical non relativistic 
treatment of POINCARE (4) and other authors. Indeed as we shall show in 
a subsequent paper resulting of common research in Bristol with Prof. Boum, 
such states are characterized by a set of quantum numbers which reproduce 
exactly a classification analogous to the classification proposed by TIOMNO (°) 
for elementary particles. This classification differs from the Nishijima Gell-Mann 
scheme by the essential fact that fermions and bosons are respectively iso- 
fermions and isobosons. 

However before following these lines of research it is necessary to give a 
Lagrange formulation of the equations of motion of such extended models in 
terms of convenient parameters. 

At first sight such a step raises many difficulties. 

If one drops the point-like character of particles and assumes them to be 
comparable to some field densities enclosed within a limited time-like tube 
(with radii of the order of 10-! cm) we are immediatly faced by the well 
known trouble that in relativity such a model depends evidently on an in- 
finite number of variables. This results from the fact that it is not possible 
to extend to the relativistic domain the non relativistic theory of rigid bodies, 
thus reducing that infinite number to a few characteristic variables (such as 
the Euler angles and the position of the center of mass). 


(*) R. HoFSTADTER: Rev. Mod. Phys., 28, 214 (1956). 
(*) H. Poincaré: Acta. Math., 7, 259 (1885). 
(9) J. Tiomno: Nucvo Cimento, 6, 69 (1957). 
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However, as we shall see, there exists a way of surmounting this trouble. 
If we consider the extended particle as a droplet enclosed within a time-like 
tube in space-time and define it in the most general way by an energy-mo- 
mentum density ¢,, and a current density j, inside the tube, we can introduce 
integral tensor quantities defined by integration over the total volume of the 
droplet. These quantities have been shown to satisfy relations which deter- 
mine their behaviour in time. In this way we describe a global sort of ave- 
rage motion of the body as a whole which leaves aside finer and irrelevant 
detailed motions within the droplet. Such an average hydrodynamical treatment 
evidently reduces the number of mechanical variables needed. 

In fact as we shali see, the general equations of motion need only a limited 
number of variables which correspond to average physical properties of the 
droplet considered as a whole (such as the total angular momentum M,,, the 
total linear-momentum G etc.). 

In principle it is then possible to determine a Lagrange function depending 
on such global parameters which gives these global equations of motion as 
necessary consequences of its Euler equations. 

This brief analysis determines the plan of this paper. In Sect. 1 we shall 
briefly recall the fundamental equations of motion of rotating fluid masses. 
In Sect. 2 we shall introduce the Einstein-Kramers « beingròssen » variables 
and introduce a general Lagrangian for the general case of rotating particles. 
In Sect. 3 finally we shall treat the particular case of liquid droplets satisfying 
Weyssenhoff’s conditions and discuss further possible extensions of the general 
theory. 


1. — According to our programm let us first give a brief summary of the 
relativistic liquid droplet theory. If one assumes that such a droplet is defined 
by a current density j, (equal to zero outside the time-like tube which encloses 
the droplet motion) and a symmetrical energy-momentum density tensor {,,; 
both satisfying the conservation relations 


(1) 
we can define the four-vector energy-momentum 


(7) 


1] 
G ar tua Vo , 


Z 


constant in time, where dv, represents the three-dimensional volume element 
associated with a space-like plane cross-section. 
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This four-vector allows us to determine a special Lorentz frame //, in which 
G® —0. In this frame we can define two physically important points: 
— the center of matter density Y, determined by the relation 


y: | jtan, | è | aide 


20 =0 


where the index 0 denotes the fact that the quantities are evaluated in 
the frame /7,. This center, which corresponds to the classical geometrical 
center of the particle, moves along a world-line within the tube characterized 
by a proper time 7 and a unitary tangent four-velocity ®, 


(u,v, =— €); 


— the center of gravity X, determined by the relation 


Xi | fe ar = [x ti A , 


nt.) | 


which moves in a straight line with a four-velocity parallel to G,. 
_ We can also define a total angular momentum of the droplet with respect 
to the center of matter density by the relation 


1 
Mu = a free == Y)tya = (x, — Yy)t a] dv. 


yy 
= 


Now the essential point on which is based this paper is that the global 
quantities G,, v, and M,, are (as a consequence of equations (1)) bound 


uv 


together by two sets of equations of motion, that is the four relations : 


(2) G, = G, = vy Oy Gi = 0 ’ 


d 
dt 
and the six relations: 

d 


(3) Mi = oe My = 00, My = Gv, G,0,, 


where A denotes dA/dr, the derivative along the world line followed by the 
center of matter density. 

Evidently these set of 10 equations do not determine completely the motion 
of the 6 M,,, 4@, and 3 v,, 80 that 3 more equations are needed which re- 
present various physical restrictions on the states of the droplet. One can 
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utilize for example the well known Weyssenhoff conditions 
(4) Mio = 0 


or various other sets such as those introduced by one of us (F.H.) (°). 

To simplify the Lagrange formulation, we can adopt the usual method of 
passing from the isolated particle point of view to an hydrodynamical one, 
that is we substitute to the isolated particle a fluid of such particles without 
interaction (*). More precisely, if we neglect the dimensions of the particles, 
we can consider such a fluid as a continuous distribution of variables mul- 
tiplied by a scalar density g which is assumed to be conservative. The center 
of matter density of each particle will follow a four-dimensional line of flow 
of that fluid, so we can multiply M,, and @, by © and we have the 11 equa- 


wu 


tions for the 11 quantities 0, g,—06,, m,= My: 


ui i ragni 


(5) o —Ÿ; 
(2) 9, — 0, 
(7) Li = Go, 9,0, » 


where the dot denotes derivatives of tensor densities along the line of flow (?) 
(for instance g,= à,(g,v;)). 
The condition of Weyssenhoff is then: 


(8) Veen). 


The next step is now to look for a Lagrange function. 

If we first concentrate on the 10 equations of motion (6) and (7) we can 
make the following remark: 

Suppose we dispose of a Lagrange function depending on parameters ql. 
If the Lagrangian contains only the first derivatives of the variables: 


Be page 


BELINFANTE (7) and ROSENFELD (*) have established the following theorem: 
the expression ¢,, = (0L/0(0,q’) ga Ò, £ satisfies the 4 conservation re- 


lations 


(9) RA D 


(6) F. HALBWACHS: Théorie relativiste des fluides à spin (Thèse de Paris, 1958) to” 
be published further. | 

(7) F. J. BELINFANTE: Physica, 7, 887 (1939). 

(8) A. H. ROSENFELD: Acad. Roy. Belgique, 18, 6 (1940). 


4 


53 — Il Nuovo Cimento. 
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as a consequence of the Euler equations on the g”s. Moreover if one denotes 

by Zi, the operator corresponding to an infinitesimal Lorentz transform on 
uy 

the variables 9g’, then the expression 


Weak 
fu = a Ziuv (0, 4) q b) 


(where the indices » and s take all possible values of 2) satisfies the 6 con- 
servation relations, 


(9') 9, e ait ty ro bu ’ 


also because of the Euler equation on the q’’s. 

We see then immediately that the 10 conservation equations (9) and (9') 
furnish exactly equations (6) and (7) provided that the #,, and the 7,,, deduced 
from the Lagrange function take the general forms: 


(10) L39010, 
(11) fox = 3m, (0) -0,(9)) - 


As a consequence, the general problem of deducing the equations of motion 
(6) and (7) of a liquid droplet from a Lagrange function depends now on our 
ability to discover a physically reasonable Lagrange funetion with corres- 
ponding 4,, and f,,, of the form. (10) and (11). 

Naturally, such a Lagrangian gives also Euler equations determining di- 
rectly the behaviour of the variables q*. From these relations we can deduce 
in principle 3 relations which complete the set (6) and (7) and play the same 
role as Weyssenhoft®s relation m,,v,= 0. The number of Euler equations is 
of course not restricted since a complete description of the motion of the 
droplet needs in principle an infinite number of variables q*. The more va- 
riables introduced, the more complete the description will be. However we 
must evidently try to introduce reasonable global variables in limited number 
which will define some sort of global average simplified behaviour of our 
droplet leaving aside details which evidently escape complete description. This 
treatment is somewhat analogous to the kinetic theory of gases or of plasmas, 
and the variables introduced play the role of collective variables corresponding 
to average regular motions inside the droplet. 

The deduction of such Lagrangians will be attempted in the next two 
sections. 


2. — Our first problem now is to introduce a convenient set of parameters 
g which will correspond to global properties of the droplet. 
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A first variable is evident, namely the velocity v, of the center of matter 
density which we shall write henceforward icb\, for reasons which will appear 
later (the index 4 is not a tensor index, but caracterises the vector). 

The Lagrange function must evidently contain a term generalizing the 
classical total energy of rotation of the droplet. Such a term should be the 
product of the angular momentum by a term corresponding to the angular 
velocity, so it can be tentatively written in the form T=} M,,o,,, where 
©, 18 the realtivistic generalization of a global angular velocity. 

Now it is very difficult a priori in hydrodynamics to define what is an 
angular velocity for the droplet considered as a whole with respect to a given 
point (such as the center of matter density) inside the droplet. 

The simplest procedure is to go to the rest frame 2, of the center of matter 
density. In this frame we postulate that there exists a three-dimensional set 
of axes D", whose instantaneous rotation describes the instantaneous global ro- 
tation of the droplet, so that it is bound to its global distribution of matter. 
(We can take for example the principal axes of inertia with the aid of the bare 
matter density jj). The projections of these axes on fixed axes constitute a ro- 
tation matrix 07 which naturally satisfies the relation : 


OO 0. 
The total instantaneous relativistic motion is then evidently obtained by 


adding to this space-rotation a Lorentz transform with the velocity v along 
the four velocity bd, v and b, being related through: 


po == aw, by? = of = (1 — 02/0?) À. 
(16, 


In other words, the three space-like vectors b, and the four-velocity b’? 
form an orthogonal set of unitary vectors moving with the body, which sa- 
tisfies the relations expressing the Lorentz transform (°): 


(A) tu ay 
bi ——1%,, 
sik 
by = 10;0;; 
LA 
be UA Cv; 
A DR PA 


and furthermore the following set of equivalent equations expressing their 


(9) T. TAKABAYASI: Nuovo Cimento, 7, 118 (1958). 
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orthogonality and unitary character: 


&) Enti (n) (©) pi) 
(12) bY == e Esa Oy bz bg 
‘ HE AI 
(13) BED = 6, 
(È) p(n) __ SEN 
(14) db, bi = 0°", 


(where £,,,; or e is the completely antisymmetrical symbol of Levi-Civita) 
eas . . D . 
The instaneous angular velocity associated with this moving frame 2, 
is just (°) 


(15) Oy = + (bi pe — be pe) 


and we assume that this is the term which we must introduce in the kinetic 
energy T and which defines quantities ¢,, canonically conjugate to the m,, 
Physically these assumptions could be justified by saying that during a very 
short interval of proper time the droplet rotates in the rest frame like a rigid 
body. 

We can define directly such quantities in special cases. Suppose for instance 
that the current density j, inside the droplet varies very slowly. Then if we 
denote by «, the unitary four-vector parallel to j, it is well known that one 
can write: 


O, = }(9 Mc 0, U,) (1°) i 


But this is not true in general since the droplet can contain subvortices and 
many complex motions so that one cannot suppose that « ,, varies slowly over 
the droplet as a whole. 

However we can still define such a quantity even in more complex cases. 
Suppose for instance that the current motion inside the droplet can be split 
into a regular slowly changing part, and an irregular part. If this latter part 
is sufficiently complex, we can neglect it in average and define w,, With the 
help of the regular part. 

Independently of these physical considerations, we can say that the com- 
ponents 6 constitute a set of variables attached to the global kinematical 
properties of the droplet. These variables where first introduced in physical 
theories by Einstein and Kramers and so we shall call them henceforward 
EK variables. 

A next step is to express m,, in terms of suitable tensor variables. This 
can be easily made in relation with the preceeding EK variable We AS was 


(°) L. BRILLOUIN: Les tenseurs en mécanique et en élasticité (Paris, 1938). 
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shown by Boum and one of us (J.-P. V.) (*) one can use 2°° to split m,, in a 
sum of two bivectors, namely: 


(16) My, = Ey ap8. bs + io(bt,— DL), 


nov v n 


where the vectors s, and 4, are respectively: 


n pi 
(17) Si So Eos Mag , 
a 
(18) ty = 0 Miu bi . 


It is then clear that s, and 4,, which lie both in the space of 2, define a set 
of six independent variables which determine m,,. 

Before we go any further in the Lagrangian, let us make a physical remark. 
It is not astonishing that m,, and «,, are in general different skew tensors. 
As we have seen, it is a characteristic property of liquids that conjugate 
variables are not parallel in general. Thus g, and v, are not colinear, except 
in very particular cases, such as the spinless fluid. Physically, it can be 
understood in the following way: the angular velocity moves over the droplet 
and changes its shape as it rotates, so that the internal tensions change in 
general, modifying the angular momentum in a complex way. The special 
case where the two tensors are parallel, which is the general case for rigid 
classical bodies, cannot be used in relativistic theories El 

To complete the Lagrangian, we must evidently add to T the term cor- 
responding to the rest mass energy of the matter: 097(,c?. 

The sum of these two terms corresponds evidently to the total energy. 
Furthermore we must add a certain number of terms with Lagrange multi- 
plyiers imposing restrictions on our variables. The first term imposing the 
conservation of the current (5) can be written: 


icob 7,8 . 


And then we must have a term imposing unitarity and orthogonality to the 
b BK variables, that is 1,,(b5b; — 6,,) (4, is evidently symmetrical). Fi- 


pri 
nally, we get the following Lagrangian 


(19) L= 0M0 + Ole pap Sad? + i (Opty — BE, )] OY DD + 
+ 100d 0,8 + Ay (BBY — dyn) » 


(*) In a sense this is analoguous to the well known fact that in the neutron case 
the charge repartition is such that there is no total charge but there remains an ano- 
malous magnetic moment. 
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in which we frequently replace for the sake of simplicity the bracket by m 
and ich, by © 

Batore we den the equations of motion, it is necessary to remember that 
the droplet is describable by a set of 10 dynamical yariables (G,, M,) and 
6 kinematical ones describing the angular velocity. These 10 variables plus 
Des satisfy the 10 conservation equations, so that these equations leave 
3 degrees of freedom open. In fact there exist 3 more degrees, because of 
the br, which do not appear explicitly in these equations. Now any given 
Lagrangian determines all variables through the corresponding Euler equa- 
tions, that is, in our general case, 16 variables. As we have shown, Weys- 
senhoff’s case is very special in the sense, that 3 relations more between 6? 
and m,, were introduced, and must be derived as Euler equations from the 
Lagrangian. 

As a consequence the 13 equations (6), (7), (8) can be integrated separately, 
but the other 3 variables b' are also determined according to the 3 remaining 
Euler equations which evidently depend on the explicit form of the Lagrangian 
chosen. Of course, many Lagrangians give the correct set of the 13 relations; 
so that Weyssenhoff’s motion for the center of matter density can be asso- 
ciated with many modes of internal motions of the liquid droplet, which cor- 
respond to various possible distributions of the internal Poincaré stresses which 
hold its elements together. 

Let us now write the equations of motion. First the Euler equations 
relative to the variables S and 2, provide immediately the condition of con- 
servation for 0, and the condition (13) for the set of EK variables. Then we 
can use some of the Euler equations to eliminate the redundant variables, 


namely the Lagrange coefficients 0,8 and 4,,. The equation relative to the 
variable 6 gives: 


Me? a IO 
that is: 


Sh QS = — eme — Em, F202. 
Then the equations relative to the space-like EK variables eC” give: 


0,(4m vu © vb” sa im b+ 27 bp” 


(2) Ye 4 


that is: 
(21) | 2,07 + rid — 42d = 0, 


while the equation relative to b gives in the same way: 


(km, 0,0) Li im, Oe + 27 ges +tofe 8,58 8,9 bé 4 it Ab boo — bb®)) 


uv u 
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or else: 
99 9 7 (4) . (OU ey (4) 3 (E) £) 
(22) m0; + m,, 5 44,30, + + 08,4589; 8,05 ay 


wa 


— git (bb — BP DE) + iom,,b© 8,5 — 2ico 9,8 = 0. 


Ss 
7 a 4h 


Multiplying (21) by bi (with contraction with respect to r) and (22) by be 
and adding, we get, taking into account relation (13): 


2m,,0 DE + m,,— 44, + 0 


1(&) (£) p.(4) 
wa 2 uv Eva 03 8,05 ë, + 


—oùt,0 0 (DLE — bP DY) + ibm, be de — 2icb, 00,8 = 0. 
Taking the antisymmetrical part, we eliminate 1, because of its sim- 
metrical character. Multiplying the resulting equation by b®, we get, after 
a short calculation in which we use equations (18), (13) and (20): 
(23) 2ico 9,8 = — 2970 bi — 2io(i, + 1,05 bY) + 


+ som, gb 0,09 + 06,80 8,05 (Sy — be Ds) - 
An important consequence of this relation is that: 
(24) | PIO 


Now we can easily show that the Lagrangian (19) gives the correct results, 
for the tensors 4, and, f,,, take the forms (10) and (11). 


VA 


OL OL 
a Y fa ne fe) (5) 
buy (0,8) 6,8 ae (0, BS”) [21 


because of (24). An immediate calculation furnishes: 


ine 9,10 08 Ar im, bi 9,67) 


Mm x 
so that we have a linear momentum density : 
(25) I Al 0,8 a 3M,, bY a,b; > 


One sees immediately that the form (25) is the generalization of the clas- 
sical relativistic expression for the momentum. S is just the relativistic J akobi 
function, while the second term represents the contribution of the internal 


rotation. 
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If we replace 0,8 by its expression (23) we get: 


; ; Q E), we (4), (4) 
: CY ay DLE 0 
Gu = icoMTtb, e (ts + tb) di, ) Die EyanpO, D; ( a wr ) 


That is, for the so called « transversal momentum »: 
4 i E “E 
(26) p,= 07m — fa = ; t, +00 bl + à Evinp 0, 8308 (Ow — D DE )| . 


Of course the conservation equation ¢,t,,= 0 provides as usual the hydro- 
dynamical equation: g,= 0. 

Let us show in the same way that the Rosenfeld-Belinfante tensor f,,, 
takes the correct form (11). 

We have 


0 1 
GS == 5 May icb DE 
IG 5 


The operator of an infinitesimal relativistic rotation for a vector is, as is 


. well known: 


3% = 4(6,,0 


uv au Bv 


— 6,,6,,) 


a” Bu 


so that we have immediately: 


eran ey (4) BE) ee CE) ET EME (E)\ HE) 
funn = +i0m,,b, bY (5,45 ay O09 ay) D5 7 £0, (m,,,b; m,,b, LE 


Or because of relation (13): 
ge = 3M,,0, b) 


which is just the correct form. 

As we noticed above, we can obtain from our Lagrangian 6 relations too, 
concerning the 6 kinematical quantities whose evolution is not determined 
by the 10 dynamical equations. Namely: 

Derivation with respect to t, gives: 


ied) bY (0, bb — 0,0 bY) = 0, 
that is 
(4) {1 (E) BE NOT RE) __ 
b (oP bP — bi bY) 220 


or in other words: 
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Derivation with respect to s, gives on the other hand: 


eh à più peo) (a) 
ich; 0,07 OF E4490, = 0. 


That is 
(4) RÉ RE __ 
Ewap0, by bg == 0. 


Or, in other words: 


Exp O30 y — 0. 


Therefore the tensor ©,, is zero. 

So we see that the above Lagrangian formulation is not at all equivalent 
to the bare set of the 10 dynamical equations expressing the general motion 
of the relativistic droplet. It represents a very particular case of that motion, 
namely the one where the matter of the droplet, considered in a purely kine- 
matical sense, has no instantaneous rotation. In this point of view, such a 
motion generalizes in a natural way the classical motion of the newtonian 
point-particle. 

It might seem strange at first sight that one should have an angular mo- 
mentum density without kinematical global internal rotation. But à close 
examination shows that m,, and w,, are very different kinds of physical quan- 
tities. m,, depends on the internal stresses of the droplet, so that if these 
stresses vary with time in such a way they would not disturb the global distri- 
bution of matter, one could have angular momentum without angular velocity (*). 

Taking account of the above supplementary relations in the expression (26) 
of the transversal momentum, we get simply: p,= (o/¢)t,: In other words, 
in the case considered, the non classical transversal momentum does not result 
from any internal rotation, but only from the variations of the vector ¢,, 
which represents, as is well known (') the distance between the center of 
matter density and the so-called «center of mass ». 


3. — According to our program, we shall particularize the general motion 
defined by the dynamical equations (6) and (7) with the aid of some supple- 
mentary conditions. The most interesting way is to state some relations between 
the dynamical angular momentum m,, and the set of kinematical vectors b\°’. 
The simplest binding relation is naturally Weyssenhoff’s one, which assumes 
that m,, is orthogonal to the four-velocity. 

As is well known, such a relation added to the 10 preceeding ones, fur- 
nishes a set of 13 equations for 13 variables, these equations can be completely 
integrated and it seems that the motion is determined in an unique way. But 


(*) This case corresponds to a complete decoupling between M,, and ,,, which 
is evidently not the case in general. 
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we note that no statement has been made about the 3 other EK variables, so 
they are entirely arbitrary. 

Now to treat this problem by Lagrange’s method, we must add to the 
Lagrangian (19) a term 2,t,, which will immediately furnish by derivation 
with respect to the Lagrange multiplyer J, the condition needed. 

But then we get a supplementary term in the Euler equation relative to #,, 
which becomes: 

(ER l, = 0 
and then does no longer represent a restriction on ©,,, but gives the expres- 
sion for the new variable Z,: So we are left with the restrictive condition 


(27) Eva Da Vy = Ox 


If we express w,, with the aid of the EK variables we have: 
Ego = 0 
and we see, multiplying by 5° and taking the dual, that we have: 
(1) DA 2 (A 2 
(28) Ud 


In other words, the derivatives of the three space-like EK vectors are 
parallel to the four-velocity. Or else, noting that with regard to the ‘proper 
frame 27, of the center of matter density, ©,, has only time components, we 
see that the kinematical motion in the frame Z, is purely accelerative without 
any rotation. 

Furthermore, multiplying (27) by d°° and taking the dual, we obtain: 


Bk) JG) LR UE) 
bb — bb” = 0 
that is 


(29) UO = 0 


the set (28), (29) beeing equivalent to (27). 
If we take into account the latter relations, and remove the terms con- 
taining t,, which vanishes, the expression (26) becomes 


lo : 
Pu = oe) ic Enta80 4 8,08 , 
that is 
iL 
# Pu So, e Muy ? 


which is a well known consequence of Weyssenhoff’s laws of motion. 
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Finally, we can impose a further restriction by the assumption that the 
dynamical spin vector s, is parallel to the kinematical vector db. We must 
then add to our Lagrangian a further Lagrange term, 


(3) 
7] uv & uv ab Sa bp 


so that the derivation with respect to the Lagrange multiplyer 7,, (which is 
evidently antisymmetrical) gives the condition of parallelism: 


8 b® — 0. 


E uva B Sap 


Then the condition (27) becomes 


Ew oa COM E + NapË OB | D — 0. 
This introduces no longer any restriction to the motion, but simply pro- 
vides the expression of the new variable 7 
The total Lagrangian is then finally: 


tuv° 


2 = ome + € 5 lEmap8a bp + i (bet, — By? t,)] D2 UDO + 
+ icob® 9,8 + 140? bY — by) + lit, + QuvEprapSadp è 


We note that, with the latter assumption, the angular momentum density 
takes the form: 
mM = ohyé bY DE, 


(2 uva a 


where h, designes the magnitude of the spin s,. 
Then the rotation energy becomes: 


T= E Olin Egy, be 09 bY DE. 


jaa B 
The terms corresponding to €=3 and £— 4 vanish 
T = oh bo (bP OY + bP bP), 


Ox p 
which becomes, taking into account the relation (12) 
T= ohb = icohbS bi? dd. 


Then, removing the redundant variables s, and t, and the corresponding 
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Lagrange conditions, we get for the Lagrangian the most simplified form 


e = 0Me + icgh, bb 0,0% + ice b\? 0,8 + À,,(b% D — du) » 
which was already proposed by two of us (F.H. and J.-P. V. (*11)) and which 
appears to correspond to the particular case where the angular momentum 
and angular velocity are related in a very strict way. It is easy to show that 
in this case the expression (30) is also valid. 


4. — Conclusion. 


In conclusion we wish to indicate a few possible extensions of the theory 
developped in the paper. 

Evidently this Lagrange formulation opens the way to wide developments 
of the relativistic hydrodynamics. Introducing more variables into the La- 
grange function, it is possible in principle to describe more and more properties. 
of the isolated liquid droplet. 

Another possible extension is to introduce coupling terms between these 
droplets in order to develop the theory of relativistic hydrodynamics for fluids 
presenting internal angular momentum density (spinning fluids). The Lagrange 
formalism is very convenient to treat such theories, since it is only necessary 
to introduce in the Lagrange function suitable coupling terms. Some steps 
in this direction have already been made by one of us (F.H.) (°). 

Another possibility is to study the isolated droplet in special physical con- 
ditions. In a further paper, for instance, we propose to study the relativistic 
hydrodynamics of rotating fluid droplets with negligible rest mass, moving 
with velocities approaching the velocity of light. 

The Lagrange formalism also offers the possibility to quantize the pre- 
ceeding theory in the conventional way: one deduces the hamiltonian from 
the lagrangian and replaces the Poisson brackets by the usual commutators. 
In this way one discovers a new class of relativistic wave equations corres- 
ponding to all possible states of relativistic spinning droplets. Some of these 
equations could eventually correspond to new types of particles which have 
been shown to exist in nature. 

Finally, and this seems to be the most promising line of research, the pre- 
ceeding model paves the way to quantize internal states of motion of relativistic 
fluid droplet. This results from common work made in Bristol by Professor 
Boum and one of us (J.-P. V.). 


(1) B.C. UnaL et J.-P. VIGIER: Comptes-Rendus Acad. Sciences, 245, 1787, 1890 (1957). 
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If we define internal stable quantum states by the physical property that 
the angular velocity ,, and angular momentum m,, come back periodically 
in the same relative positions inside the droplet, one can show that such states 
determine a new set of internal quantum numbers which characterize dif- 
ferent types of periodic motions of these two quantities. As we shall see (!) 
these quantum numbers are very similar to the numbers introduced, by NISHIJIMA, 
GELL-MANN and TIoMNO to classify elementary particles. 


* * * 


We wish to express our gratitude to Professors DE BROGLIE, BoHM and 
TAKABAYASI. This paper is a natural development of common researches 
carried out in Bristol with Professor BOHM; and professor TABAKAYASI’s work 
and advices have greatly contributed to the clarification of the basic con- 
ceptions used therein. 


(2) In a further paper. 


BLASS UN AO! (*) 


Si studia la teoria generale delle masse fluide in moto relativistico in termini di 
un gruppo di variabili originariamente proposte da Einstein e Kramers. Le equazioni 
generali del moto impostate direttamente da Bohm e uno di noi (J.-P. V.) appaiono 
allora come equazioni conservative dedotte da un lagrangiano scalare di queste nuove 
variabili. Questo procedimento apre la via ad una quantizzazione del moto interno 
di un modello di particelle estese in uno spazio reale; modello fortemente suggerito 
dai risultati sperimentali di Hofstadter. 


(*) Traduzione a cura della Redazione. 
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Hyperfragment Decay Energy and Angular Distribution - I. 


J. SZYMANSKI 


Institute of Theoretical Physics - Warsaw University, Warsaw 


(ricevuto il 19 Settembre 1958) 


Summary. — Impulse approximation is applied to the hyperfragment 
decay. x-mesons energy distribution is obtained for the *He, case neg- 
lecting x-« interaction. Results are compared with experiment and the 
validity of impulse approximation is discussed. 


1. — Introduction. 


Decays of bound A, both mesonic and non-mesonic, can serve as a source 
of a large quantity of data on hyperfragments, and, indirectly, on the inter- 
actions of A with nucleons. From the analysis of the decay we can draw con- 
clusions concerning the role of singlet or triplet potential as suggested by 
DALITZ (1) and CHow-Stroxov (?). On the other hand, it seems that much 
valuable information can be obtained by analysing the angular and energy 
distribution in the case of the three-body hfr decay, e.g. ‘H,>t+p+x. If 
the strong interactions of the decay products can be successfully taken into 
account, it may be expected that with somewhat better statistics than those 
at present available it will be possible to check the wave functions obtained 
by calculating the binding energy (#1) and indirectly the A-N potential. In 
this paper we deal merely with the mesonic decays, particular consideration 


. H. DaLITZ: The nuclear interaction of the A-hyperon (preprint). 

C. Crow and M. J. Stroxov: Nucl. Phys., 6, 10 (1958) 

ZIELINSKI: Nuovo Cimento, 3, 1479 (1956). 

SCHNEPS, W. F. Fry and M. S. Swami: Phys. Rev., 106, 1062 (1957). 
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being given to the interaction of the final decay products; that this influence 
is important is indicated by the angular distributions as regards cos 0 = pk 
in the apparent rest system of the A in which k=—q; here k, q, p is the 
momentum of the z-meson, proton, and the recoil, respectively for the decay 
hi—R-+p-+7. By à we denote the versor with the direction of vector a. 
For this decay we have a distinct anisotropy and asymmetry independent of 
the hyperfragment. If the interaction is neglected and it is assumed that A 
has a spin 3 (which now seems rather certain) and is bound with an orbital 
momentum /= 0 we obtain e.g. for ‘He, an isotropic distribution. If A were 
to be bound with a higher orbital momentum, which is not very probable, 
then we can obtain an anisotropy, but it must be symmetric with respect to 
the plane cos 0 = 0, while experiment (5) shows a distinct asymmetry. Eight- 
een cases have been found in which cos 6 is negative and only one case in 
which it is positive. 

In principle, our calculations, which involve the decay products, coincide 
with the impulse approximation. The modification introduced is that we give 
the method by means of which the virtual protons may be taken into account, 
these protons becoming real ones only after scattering. The ratios of the cross 
sections for the scattering of protons and pions on light nuclei allows one to 
assume that the main interaction will be that of p with the core, and there- 
fore in the present calculations too we limit ourselves to consideration of only 
this interaction. In addition to neglecting the interaction we make use of the 
following other assumptions: a) in light hyperfragments Z < 2 the greater 
part of the time is spent outside the core (57); b) if the decay occurs outside. 
the core then the operator for the A+p-+7 which already takes into ac- 
count the basic part of the strong interactions has the form 


(1) H, = [yo + Bok) pp, AT - 


where g and B are constants that are independent of the energy. For a free 
A the only possible shape of the transition operator is 


(2) Hi = of pols + vs) apr AT . 


(5) R. Levi-SeTTI, W. SLATER and V. L. TELEGDI: Reported on the VII Roches'er 


Conferen-e, 1957 pag. vill-9. A world survey of experimental data on hypernuclei preli- 


minary report preprint. 
(6) N. DaLLaporra and F. FERRARI: Nuovo Cimento, 5, 111 (1957). 
(7) H. WiLteLmson and P. ZIELINSKI: Nucl. Phys., 6, 219 (1958). 
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This follows from the requirement of relativistic invariance (see e.g. STAPP (8), 
J. WERLE (*)). The non-relativistic approximation of (2) is (1). 

The adoption of the form (1) for the operator of the vertex of a weak 
interaction corresponds to an artificial separation of the decay process into 
two stages. In the first stage we have only p interacting with 7, in the second 
we have p interacting with the core and x with the core. This separation seems 
justified if we note that the range of the x-p potential is of the order of the 
Compton wavelength of the nucleon. We shall now discuss assumption (a) 
is somewhat more detail. For this purpose we expand the wave function of 
the hyperfragment in terms of the eigenfunction of the core 


(3) ie) == oy Xo(&) for) + YS WnXald) LAr) = Pi+P, 


n=1 


here £ denotes the internal variables of the core r=r,—r,, r\andr, de- 
notes the corresponding co-ordinates of the center of mass of the core and 
the A particle. The functions f,(r) are defined by the formula 


(4) fn (r) = W [Voir 28) a8 


N is the normalizing coefficient and «, the coefficient of fractional parentage. It 
is seen at once that the function f,(r), n #0 can be essentially different from 
zero only in the region of the core because the range of the A-N forces is 
small, and therefore a deformation of the core may be expected only when A 
is outside it. For nuclei with Z= 2 the probability of mesonic decay inside 
the core is decreased by the Pauli principle. Besides, inside the core there 
is a competitive non-mesonic decay because of the direct proximity of the 
nucleons. These the facts which allow us to suppose that even if there is some 


core deformation the chief contribution to the mesonic decay is made only. 


by the first part of the function (3) = %yo(&) fo(r), since f,(r) is essentially 
different from zero outside the core. 

Our work has been divided into two parts. In the first part we discuss the 
computational method and we give the results obtained for the energy distri- 
bution of the z-mesons from the decay *He,+>«+p-+z7. In a forthcoming 
paper, which will be published after completion of the numerical calculations, 
we shall discuss the angular distributions for the decay of light hyperfragments. 


« 


(5) H. P. StAPP: Phys. Rev., 103, 425 (1956). 
(°) J. WERLE: Phenomenological theory of the S-matrix (to be published), 
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2. — Method of calculation. 


The Schrédinger equation for the decay of a hyperfragment to an ac- 
curacy of g/g is the coupling constant of weak interactions and has the form 


(2s HH, 


here 
A-1 
H=H+K+K+WM+W W=d Vs 
i=1 
e wh A-1 A 
H,=Y Kn+3YVin=K+V Wie Vs 
n=1 nam i=1 


where A is the atomic number of the hyperfragment K,, K_ are the kinetic 
energies of the n-th nucleon and the x meson respectively, Vin. Vig are/thé 
potential energies of the i-th nucleon with respect to the j-th and the x meson 
with respect to the i-th nucleon. 


In the general case (4) the solution has the form (see e.g. CHEW-GOLDBERGER) 


1 


(5) Se re pena 


HP, 
where ®, satisfies the equation (E — H)D, = 0. 

In our case ®,= 0, which means that initially we have only hyperfrag- 
ments. We make the transformation 


iz di 
— H, — Hy. 
E—H +in a — w 
where 
(6) H=E HE, -K,+", W=W+ Wa, 
L pi 1 RC | 
Si a = H, — 
FF NA ww) a+] A 


1 


1 1 ni 
4 vl A, tem =e! + T-|H,. 
w |a a a a 


We denote by T the scattering operator defined e.g. in the papers of LIPPMAN- 
SCHWINGER (1°) or CHEW-GOLDBERGER (11?) 


T = 0 + w 
a—w 


) B. Lippman and J. SCHWINGER: Phys. Rev., 79, 481 (1950). 
(11) G. F. CHew and M. L. GoLDBERGER: Phys. Rev., 87, 778 (1952). 
) G. F. CHEw and G. C. Wick: Phys. Rev., 85, 636 (1952). 


54 - Il Nuovo Cimento. 
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The probability of a transition on the basis of (6) is given by the formula 
il 
(6) W = 2a0(E — E,)|M |? Me GI(1 +27) H, | Yo} - 


For the sake of discussion, in our further calculations we shall limit our- 
selves to the decay *He,—a+p+z. We denote the moments of the x par- 
ticle, proton and x meson in final state by pqk and q'p'k' in the intermediate 
state; we take Y in the form 


= | 
(7) VASE | GA Hea) Ian) - 


We omit here the indices describing the state of the core, remembering that 
all the while we have to do with the basic state of 4He. X, is the center of 
mass momentum of the hyperfragment, p, is the À momentum with respect 
to the core, f(p,) is the Fourier transform of f,(7), » is the projection of the 
hyperfragment coinciding for *He, with the A spin. By (7) and (6) we obtain 


1 © 
ab) + in 


(8) M — [<pale (PAK oy ee +3 |[cpake|T|p'q'k’> BE 


*<p'q'k' | A, | pykv>-f (py) = dp’ d?k' dq’ dp, . 


A. 
(270)?? 
Introducing the variables 


M, 


M, 
DoT bee da t=p+qK=p+q+k, M=M,+M,, 


neglecting the interaction between x and «, and taking H, in the form (1) 
we obtain 


7 


(9) <pqk"| A, pk) = g(2x) 6(q + k— py) 6(k— ky) U (1 + Bok)u 


(10) <pqku|T|p'q'k'u'> = 2x) d(k—k')d(e—t)Ku|T|lu'), 


(11) maar 


pi 
—— 
4 


ED) TETE) 
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where 
MM. M, 
D ET TA M re PRE CSS RS 
1 1 M, M. 
E) =—P, EU) =— 12 = Spied Pre DE 
D a a 


In order to integrate over 1 we put 1/ (EQ — E(U’) +in) in the form 


(12) — ind(E(l) — E(l')), + 


El) — BV’) 


Making use of the relations between the phase shifts and the diagonal in 
the energy part of the operator 7 (see ¢.g. LIPPMANN-SCHWINGER) we obtain 


i dl’ 
E E ALE ED O) — BE) ul + fol) flo) 


vp Pe ñ 
= D sin bn exp [1040] CT MM — jou) CT MM — plu") Yu (À - 


UU JnM 


Yow» (l) (EQ) — EV) u(1 + Bok)u,f(l — v)odE(l)}4Q,. 


We shall now show that if the decay takes place outside the core then the 
part of the matrix element containing the integration in the principal value 
can be expressed by means of (13) by a known function D,(1V). To show this 
we transform into the space configuration that part of expression (11) which 
contains the operator 7. We make the transformation only with respect to 
l and I’, that is, the center of mass of the « particle and the proton; we con- 
tinue to describe the meson in the momentum representation. After trans- 
formation (11) has the form 


(14) Zara Iu re ><r|L [ra Gr") + @,(r'r")]: 
i ‘exp [ivr'] uy(14+ Bok)u,fo(r") , 
sin l|r'— r"| 
LE; | ray | 
cos l|r’— r’ 


Gr") = <r'u'| — id (E(B) |r'u'> = iA 


G,(r'r") = SR erp |r'u'> + 


rr" 


| 
(15) | 
| 


l' denotes the operator whose eigenvalues were previously denoted I' 
(i= i(2/0r)). We note that if r or r’ is greater than the range of the potential 
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W, then T,,(rr')= <ru|]T]r'w) vanishes; this is seen at once from the formal 

expression for 7. Therefore if by À, we denote the radius of the x particle 

obtained from the scattering of nucleons on x, then we can assume that 

Ty = 9 ifrorr>R. Ifr<r" we can expand (G(r'r’) and G.(r'r") into 

a series of Legendre polynomials, respectively cos y = r'r', in the form 
G(r'r") =iA > (2n+1)jn(lr"")jn(lr") P.(co8 y) ; 


(16) i 
G.(r',r")=-AY(2n+1)jn(r"'")n,(tr") P,(cos y) . 


Inserting (16) into (14), expanding exp |[ivr"] into spherical harmonics and 
integrating over dQ,,, we obtain, if the decay occurs inside the core 
Ry 
(17) BY I inclu|rp)<ru|T|r'W!Yjn)Yam®)Yn(È)uy(1 + Pok)u, dir der’ 
we 
0 


. UOTE r'2 ar" fior} : 
Ri Ro 


From (17) we find the ratio of matrix elements g,, and g, containing the 
respective parts of the Green function G, and G, proportional to P,(cos y); 
we obtain 


[in (or), (Ir) fo(ryr? ar 
(18) aa ee it si D, (Ww) 
"Sin (ir) in(or) folr)r? dr 


Ro 


By virtue of assumption (a) we neglect the decays taking place in the region 
r"< R,. By (18) and (13) we obtain 


(19) M,, = lu, a + fak)u,jl — v) + i] > (1 + D, (tv) sin 6,, exp[iô,,]- 


Inu 


Gis(J MM — up)eny (J MM — p'p') Tei Yrar-w(0) (1 + fok)u, 


{A — v)5(E(l) — EU) dB’) da i 


Integrating over l’ we have 


(20) M,,=9 {u,(1 + Bok) u, f(l—v) +4 > (1+D, (Iv)) b, (Lo) sin è, exp[iòd,,]: 


InM 


Ong (TMM — pp) On (TMM — pp) Vie 10 Yo x 0} 
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and 


(21) b, (lv) = 2a perni Ji l—v)Y,(cos d)dcosd. cos.ìd=T%. 


3. — Energetic distribution of the x-mesons. 


Since we are not interested in the absolute value of the probability of 
transition and only in the relative energy distribution, we shall omit the energy- 
independent multiplicative constant. The probability of transition as a function 
of energy is given by the formula 


(22) (0) Y My 2Ô(E — Ey) 010,02,d2, dE, . 


uv 


In the center of mass system we have 


2 a hats ay eels Me 
ere 1= [2,0 faite HT LM: 


Q being the heat of the reaction. Integrating over dQ, 2,dH we obtain 
(23) W(o) ~ 2F(0)(1 + 2k?) + = ‘a + Re D,,(lw)) bi(lv) sin? ô + 
+3Im D(lv)bi(lv) sin 23,2 (2J + 1) + pr keG Jn) + 


+ Y|1+ D, (lv) 2b; (Ww) sin? 6,, È i ta 5 PER ( Im]. 


Jn 


Here 
F(v, 1) =| ire —v|?dQ,dQ, , 
ue ; 9 2 
jesi (2,(L000) On4(7 MM — up) 

4) Gna E I 
Ca MM — u'u')Cy(1u +», 4, 7) Cyu'+ v,u',v)- 
‘Chad; M Col La +4) Onn(LMN a iar v) ’ 

as | 9 

= — — 1)4+#+1 Onn(L, 0, 0, 0) Cn (JM, M — n, u): 
(25) ROm= FS (yet pi Cans 0, 0, 0) Oni FHL, | 


On(TH, M — py 1) CIM, M—p', 10) CIM, M — pu", pw") Oye! + oy’, »): 
Oy (1u"+ vv) CAEN, M — y", pu" + v) Onn(LNM —p'p' +7). 
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On the basis of data obtained from the scattering of protons on helium 
we see (JUVELAND JENTCHKE (1*)) that only the phase shifts Oxo) ya Og: are 
different from zero in the energy interval of interest to us. 

Hence only the first three G(Jn) and R(Jn) will be of interest to us. After 
carrying out the summation we obtain the following values: 


3 2.17 6.81 
| a= 5, aan="—", an, 
(30) ‘ 
3 2.78 8.01 
| nao = 2. R(3 1) = 4x ? RGD = 


Furthermore, as f,(r) we take the function of Dallaporta and Ferrari 
f(r) = N exp[— ar], a = 561102, 


which gives a good binding energy. For such a form of it it turns out that 
D, (ol), v= lare of the order of 4, and consequently, in the first approximation 
we neglect them. This considerably simplifies the calculations. For b,{(vl) 
and F(vl) we obtain 


82? il al 
È 2 
Eee È =B) (A+ | 
4x 
by (lv) = A? B?° 
122 2 A — B 2AB 
5 Lie 
be) = NE METE am 


A=l +0 La, 
B= 2». 


And for W(w) we find, after making the above approximations, 


(32) W(w) ~ Flo) + 3bî(0) sin? dy + 3d%(v)[ sin? d,, + 2 sin? d,] + 


+ Btk*[F (lo) + 3bi(t0) sin? 64, + bi(Iv)[2.37 sin? 6,, + 7.2 sin? d,,]. 


The results are given in Fig. 1 where the curves for the parts corresponding 
to scalar coupling and pseudoscalar coupling for weak interactions are given 


“ 


(18) A. C. JUVELAND and W. Z. JENTHKE: Zeits. f. Phys., 144, 521 (1956). | 
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separately. In both cases we also give the energy distribution that is obtained 


if the interactions of the decay products are neglected. All distributions are 
normalized to 1. 


Fig. 


Number of events 


In the experimental distribution the cross hatched area 


4 


30Energy in MeV 25 0 
La P.V. Coupling without interaction. ----- P.V. Coupling with interaction. 
--- S Coupling without interaction. —— S Coupling with interaction. The theoretical 


distribution is given for the pion energy in the interval (23-35) MeV. For w smaller 
than 26 MeV, if the D,,(vl) (proton virtual states) is taken into account the theoretical 
distribution can be considerably modulated. 


represents cases that are not uniquely defined, while the remaining area repre- 
sents the uniquely defined cases. The experimental data are taken from LEVI- 


SETTI. 


We shall discuss the results in the next part of this work along with 


the discussion of the angular distribution. 


RIASSUNTO (*) 


Si applica l’approssimazione impulsiva al caso del decadimento di un iperframmento. 


- : À he 
Si ottiene la distribuzione d’energia dei mesoni x nel caso del 5He, trascurando l’inte- 
razione n-x. Si confrontano i risultati con l’esperienza e si discute la validità dell’appros- 


simazione impulsiva. 


: 3805 


(*) Traduzione a cura della Redazione. 


IL NUOVO CIMENTO ViGie wx No 1° Dicembre 1958 


Penetration of a Magnetic Field into a One Dimensional Plasma. 


A. RON 
Technion, Israel Institute of Technology - Haifa 


(ricevuto il 22 Settembre 1958) 


Summary. — In connection with the proposed compression and heating 
of a plasma through the application of an external magnetic field arises 
the problem of the penetration of such a field into a plasma. It is here 
investigated how an external magnetic field penetrates into a one dimen- 
sional bounded plasma which behaves as a compressible hydromagnetic 
fluid whose conductivity is proportional to density. The resulting mathe- 
matical problem is complicated, in particular as the movement of the 
boundary is not known beforehand. Transformation to a system of co-or- 
dinates moving with the fluid simplifies the equations. In particular, 
in the case of a slowly varying magnetic field a single equation obtains 
with the magnetic field as only independent variable. 


1. — Introduction. 


Interest in the problem of the penetration of a magnetic field into a plasma 
arose recently in connection with astrophysics and high temperature plasma 
physics. The penetration of a magnetic field into a plasma must be studied 
in order to determine its pressure and movement and to gain some insight 
into the mechanism of its confinement. This investigation is of special interest 
in connection with recent proposals (1) to confine and heat plasma with the 
aid of an external magnetic field. 

In previous works (?*) the problem was introduced under the idealization 
of the plasma as an incompressible fluid of constant conductivity. LOUGH- 
HEAD (*) was more realistic by considering the plasma as a compressible fluid 


(1) G. Scumrpr: to be published. « 

(?) E. N. PARKER and M. Kroox: Astrophys. Journ., 124, 216 (1956). 
($) S. CHANDRASEKHAR: Astrophys. Journ., 124, 232 (1956). 

(*) 


R. E. LOUGHHEAD: Phys. Rev., 99, 1678 (1955). 
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but introduces the simplification of infinite conductivity. We shall consider 
a compressible plasma of variable conductivity, but shall restrict ourselves. 
to a simple geometry. 


2. — Formulation of the problem. 


We consider a one dimensional plasma confined by an external force to 
a given interval along the z-axis (— à < x < a) and unbounded in the other 
two dimensions. At a given moment (t= 0) an external, homogeneous and 
2 directed magnetic field is applied to the plasma, and the external force is 
removed. The magnetic field, B,(t), outside the plasma is known, while its 
penetration into the plasma is to be studied. An ideal source of the magnetic 
field is assumed. 


21. Fundamental equations. — The plasma is considered to be a magneto- 
hydrodynamic compressible fluid, to which the Maxwellian and Eulerian equa- 
tions can be applied. Furthermore, we assume that the plasma is a homo- 
geneous, isotropic, neutral and conducting medium. The fields we are going 
to deal with vary so slowly that the displacement current can be neglected 
and Maxwell’s equations take the form 


oB 
(1) Vallata: 
1 À 
(2) —VxB=j, 
U 
(3) V:B=0, 
(4) j=c(E+tvxB), 


E and B represent the electric and magnetic fields in the plasma, y the per- 
meability of vacuum (M.K.S. system), j the current density, o the conductivity 
of the plasma and v the velocity of the fluid. 

Neglecting any external forces except those caused by the magnetic field, 
we have for Euler’s equations: 


(5) ef + Ve] — Vp + jxB, 
00 
(6) pier, 


where 9 and p are the density and pressure of the fluid. 
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Eliminating E and j we have from (1)-(4) 
0B i 
(7) —B + vxloxB)=Vx[2vxB), 
and from (5) 


Ow | 


| 1 
(8) Oia, + (Vo DE LE HR DE be 


The equation of state and the law giving the relation of the conductivity to 
the other quantities in question complete this set of equations. The pressure 
is assumed to be polytropic 


(9) p= C0’, 


with y= 5/3 and C, a constant. If we use for the conductivity Spitzer’s re- 
lation (5) 


(10) o6=0,7', 
where C, is a constant, then, with the equation of state, 
(11) PT; 
0 
where C, is another constant; with (9) we get 
(12) a= CQ, 
c is a constant which depends on C,, C, and 0,. 


2°2. One-dimensional equations. — In our one dimensional case, as the only 
space variable is x and the external magnetic field is 2 directed, we have 


(13) B= {0, 0, Bia, i}, 
(14) D = {o(a, 1), 0,0)}, 
(15) 0 Sole) 


(°) L. Srirzer Jr.: Physics of Fully Ionized Gases (New York, 1956). 
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This reduces equations (7), (8) and (6) to 


oB 0B ov 0/1 0B 
16 — Mes = e 
qe A on da 5s Al 
à = 
(17) d o 1 op gt: 
ot dx o dx puo da 
eo 00 dv 
18 we + - 
Ge) Non Corot 


The initial values of B and v are zero and those of @ and p are p;, p; res- 
pectively, while on the moving boundary B(x,t) is known to be B,(t). 

Equations (9), (12), (16), (17) and (18) with the initial and boundary con- 
ditions given, form a complete set which governs the problem. Unfortunately, 
however, these equations are non-linear and one of them, (16), is not homo- 
geneous in the first derivative. Moreover because of the compressibility of 
the plasma the boundary is moving in an unknown way, which can be found 
with the complete solution only. This gives the problem as formulated in x 
and t as independent variables a very formidable appearance indeed. 


3. — Transformation to new coordinates. 


31. The equations of transformation. — Let us introduce new independent 
variables £ and 7 for x and t 


(19) a = x(&, T) , 


(20) BS te, Gh. 

The functions æ(£, t) and {(£, t) are to be chosen so as to give a stationary 
boundary in (é, 7) and, if possible, simplified equations. 

Accordingly let £ and 7 be the space and time variables of a co-ordinate 


system of an observer moving with the fluid. Then 


(21) tte, The, 
(22) x = “(&,t) =È +] 06, t) dt , 


È is interpreted as the initial value of « 


(23) E = 2(£, 0) 
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and we obtain a stationary boundary 
È, = + a. 


The equations of transformation of the derivative are 


(24) © O 0,8 
Ot dv e, OE” 
(25) lo) L e 
si oa ag 0€” 
and the non-linear expression 
d © 


(26) —_ = 


that appears in (16), (17) and (18) is simplified to 0/07. 
In the Appendix (24)-(26) are derived in detail and it is formally shown 
that (21), (22) are indeed advantageous transformations. 


32. The transformed equations. — Applying (24) and (25) to (16), (17) 
and (18) we have 


© Ox © 1 cB 
27 = i 
oh AU 3 0g Cas ce) 
A av 1 oP BoB 
Ot o(0a/0E) \ GE u ay 
(29) mat 4 = 0. 


The equation of continuity (29) can be directly integrated, rendering 


(30) nie 


F(é) is defined by initial conditions as 


on 
F(é) = o(é, Ue: (5,0), 
F(é) = 
since 
o(Ë, 0) = Li 
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and from (22) 


Ox 
ES LE, ) 

Thus 
- Ox 
(31) org 
which simplifies (28): 

ot L'/0p « B 06 
(32 = = 
or 

ov 0 PB? 
(33 —= : 
Using (12) to eliminate o from (27), we get 
To Cx lo) 1 0B 
(34 5 = - 
7 a =) ai a 
and with (31) 
© da eB 
(35) at (5 si oi 
where k is 
di 

36 A 
(36) ma 


An important conclusion can immediately be drawn from (35) without 
solving the problem explicitly. Integrating (35) over € from §=0 to the 
arbitrary plane é = €' 


er 


PT 4 È 9 /eB 
aa sb coli 
or ; 

af. AB 2B 
(38) [Bae = */ (é") AU 


0 


The left-hand side integral is the magnetic flux, ®(é’) between the planes 
æ— 0 and v=2(é'). Because of the symmetry of B(é), (0B/0Ë)(0) = 0, thus 


0 0B 
(39) a OE = RE), 
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or, with the help of (24) and (31), 


dna de Des epee 
oe = a ee) o(E) da ©), 
and. 
À gery pO (gr 
(40) n PE) =i (E) a (E95 


For incompressible media it can be shown on similar lines that 


£ Oa) =" (a), 
where k’ is a constant. 

Comparison of the last two equations shows that the flux through our 
particularly chosen surface changes more slowly in the compressible fluid the 
greater the density. This conclusion is important for the estimation of mag- 
netization and demagnetization times of a compressible plasma. 


4. — The magnetic field equation in the quasi-stationary case. 


In the case of slowly varying magnetic fields we are led to considerations 
which resemble those appearing in connection with piston engines. The quasi- 
stationary case is of importance in plasma problems, just as the problem of 
the piston engine is solved under the assumption of a train of quasi-stationary 
States. 

In the quasi-stationary case 


ov 6 

SEL 
giving for (33) 

© B? 

(41) Ble +p) =% 
integrating, 

B: 
42 ato 2 
(42) p+ Da Gz). 


Outside the plasma p=0 and B= B,(t); we thus have 


(43) p(é, 2) + aa T) Pa Bot) i 
lb 24 
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If we write for (9) 


Pi y 
PE y o, 
and substitute o from (31), we have 
(44) SIT, 
P (0w[06} 


We can now eliminate 04/0£ from (35), (43) and (44) give 


(45 2up; 1/y 
de Sa — BYE, 1) 


and (35) changes to 


0 


(46) a 


B 0:B K — k 
CE Bn ped (up 

The only dependent variable in (46) is the magnetic field B(£, 7). In 
principle this equation can be solved with the given initial and boundary con- 
ditions; in practice numerical methods would have to be applied. 
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APPENDIX 
Formal derivation of (22), (24) and (25). 


In (19) we demand 


(47) E= x(£, 0), 
thus 

(48) a(t) = &+ f(é, 7), 
where 


(49) f(é, 0) = 0. 
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Now, f(&7) can be chosen so as to linearize 
(50) 


which appears in (16)-(18). 
(50) is trasformed by 


a. dt 0 dé 0 
{51) Di Ot Oc Ob OE? 
(52) CIO, PRC 


to the new independent variables, if we express the coefficients t, = 07/0t, 
E, = 0E[etete. as functions of € and 7. 


Now 
Ô(E, 7) fe t)\ = 
(53 = - 
(95) ca, t) \dEn) ? 
i.e. 
on fell ive a! 1 | ka 
Ta Us | le À | Let, — dle |— te =e 
Using (21) we get 
(55) a, 
(56) tye OG 
yielding for (54) 
1 œ 
(57) ER e aa ayes | 
TOUTES 0 1 
(51) and (52) now reduce to 
(58) lo) z CE te. 
di OE ay OF” 
(59) Re 
Ov a, 0Ë 
Substitution into (50) gives 
(60) hee ce hi 
dé Or Nay Vag) ee” 
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it is therefore advantageous to fix x, as 


(61) L7 0 9 
and thus 

i d 0 
62 pete Mee 
ou de, «of 


which linearizes (50). 
f(€, t) can now be defined: from (48) and (61) 


(63) HE ME ris 
OT 


thus 


T 


fe, =[v(g, Dar + 918), 


0 


which is reduced by (49) to 


T 


(64) fé, T) =a v(§, T) dt, 


0 


giving (48) in the form 


e(E, 7) = Ë + fo ad 


RIASSUNTO (*) 


85: 


> 
> 


In conseguenza della proposta di comprimere e riscaldare un plasma per mezzo 


dell’applicazione di un campo magnetico esterno sorge il problema della penetrazione 


di tale campo nel plasma. Si esamina nel presente lavoro come un campo magnetico 


esterno penetri in un plasma unidimensionale limitato che si comporti come un fluido 


idromagnetico compressibile la cui conduttività sia proporzionale alla densità. Il pro- 


blema matematico che ne risulta appare complesso, particolarmente in quanto il moto 
del limite non è noto in precedenza. Il passaggio a un sistema di coordinate in moto 


col fluido semplifica le equazioni. In particolare, nel caso di un campo magnetico len- 


tamente variabile si ottiene una equazione unica in cui figura il campo magnetico quale 


sola variabile indipendente. 


(*) Traduzione a cura della Redazione. 


55 - Il Nuovo Cimento. 
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IL NUOVO CIMENTO VOLANO 1° Dicembre 1958 


Covariant Description of Spin. 


C. B. VAN WYK (*) 
Department of Physics, Columbia University - New York, N.Y. 


(ricevuto il 22 Settembre 1958) 


Summary. — The most general, finite, proper Lorentz transformation 
is generated in terms of an infinitesimal transformation. This generator 
is interpreted as a relativistic spin operator which is capable of describing 
in a relativistically covariant way the spin and electromagnetic properties 
of particles with spin } and 1 and finite mass. Extensions of the method 
to particles with vanishing rest mass are indicated. 


1. — Introduction. 


According to a well established procedure finite Lorentz transformations, 
including space rotations, can be generated from suitable infinitesimal trans- 
formations. This procedure has the familiar advantage of linking these trans- 
formations of space and time directly with the spin angular momentum operator. 

However, it has become customary to consider only those types of Lorentz 
transformations that leave invariant either two space vectors (four-vectors 
with vanishing time components) or one space vector and one time vector 
(four-vector with vanishing space components). The infinitesimal transfor- 
mation generating such a finite transformation corresponds therefore only to 
a special type of relativistic spin operator that usually takes a non-covariant 
form. 

In this paper the most general finite Lorentz transformation containing 
six independent parameters will be constructed in terms of a generator or spin 
angular momentum operator which is used as basis for the relativistically co- 
variant description of the spin of particles with spin 4 and 1. This Spin oper- 
ator contains a six-vector which, in the case of a particle with electromagnetic 
properties, is identified with its magnetic and electric moments. 


(*) Now at the University, Bloemfontein, South Africa. 
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The covariant treatment of the spin of particles with vanishing rest. mass 
presents a more complicated problem and two approaches to this problem are 
indicated. 


2. — Discussion of an antisymmetric tensor. 
The four-vectors m, and n, satisfying 


(1) MyM, = N,N, = 1, M,N, = 0, 


where 
Ma = VMN, » ha = 1 (m, and n, are real) 
define the antisymmetric tensor 
(2) Lu=— UMN, — MyN,) = Moltun)ag Mp 
and the dual tensor 
LO, = FE uvoodioo = — ÎEuvooMoNo 5 


where €,, is the antisymmetric Levi-Civita tensor. These tensors have the 
invariants 


(3) Lili 09, Digg Ligy = = 2 


It follows that 


(4) Liyy(my, + ny) = + (My +in,), 
and 

(5) Lek = 0, Lido 0, 
where 

(6) Kn, = by Ny = Oy in = ln = 0, 


Without loss of generality g, and k, can be chosen to satisfy 
(7) kudu = 9, GG, =], TAI 


Then the dual tensor can be written 


(8) Di = Luk Qky ; 
and corresponding to (4) and (5) 

(9) Likyt dy) = + (fut); 
(10) Li,m= 0, Lo, = 0. 
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The vectors 


(11) ut = (m,+in,)2*, ot = (ky + du)2 


” 


are null-vectors and satisfy 


Ears osato Lg Ee 
(12) Hu = = uso, = 0 
Define the real three-vectors 
(13) a=m)/n, b= mn, — nM , 


where where x denotes a three-vector as well as the space part of a four-vector 
%y. Then (3) becomes 


(14) ab=0, @—b'=1 or (a+ib}=1. 


The antisymmetric tensor L,, defines a skew-symmetric 4x4 matrix 
L= (L,,) such that (?) 


0 — ig 1s bi 
143 0 ae > a; bs 
L => x I 
7 * 143 da; 0 bs 
di —- b, — ds 0 
(15) oS in RE 


where, in terms of the matrices defined in (2), 
K; Ca ibys ’ S; = tim ’ 


(j, k, mis a cyclic permutation of 1, 2,3) are the generators of Lorentz trans- 
formations and space rotations respectively and have the commutation pro- 
perties 


[tuys tool ra U(tnoOve + tao 0a Cho Ore cri tyo0 0) . 
The matrices S are the relativistic analogues of s of (1). Similarly the dual 
tensor defines the matrix 
L? = i(S:b— K:a). 


These two matrices satisfy 


(16) LP 0 = Pig ie ye 


(1) This part of the discussion follows the pattern of C. B. van Wyk: Nuovo 
Cimento, 6, 522 (1957). 
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Since they commute they have the common eigenfunctions (11) which can 
be written 


Ux 
(17) : 
Cee 2-#(N, + ON,, + ia) exp[ + fl, 


Fr 24(N, + taN,, F b) exp[ + ia) 


where N, are the unit vectors 
(18) IN; = (bb; N, = (a/\b)/ab , N; = aja. 


and x and # are undetermined parameters resulting from the lack of uni- 
queness with which m,, n,, k,, q, are determined by given L. In spite of 
this lack of uniqueness the structure and eigenvalue properties of L and L? 
are most clearly exhibited by (2) and (8). The representation (15) of L has less 
structural simplicity but a and b have direct physical significance as will 
appear below. 

From (16) follows that 


exp [iL0]=— 1 + iL sin 0 + L? (cos 6 — 1) 


and similarly for exp[ZL?p]. The most general, proper, single-valued Lorentz 
transformation containing six independent parameters can now be written 


(19) æ'— exp[iL0 + L’olx, With 4, = ict, 


where 0 and g are real « angular » variables. 

The effect of the transformation (19) is easily understood in terms of (4), 
(5), (9) and (10). The parameter 0 is the «angle of rotation » in the plane 
of m, and n, while g is the hyperbolic angle measuring the « velocity » for the 
Lorentz transformation in the plane of %, and q,. If b—0 then (19) des- 
cribes a rotation about the space direction a as well as a Lorentz transfor- 
mation relating to a velocity with direction a. 


3. — Application to particles with finite mass. 


31. Particles with spin unity. — In interpreting L, the generator of the trans- 
formation (19), as a relativistic spin angular momentum operator, one is guided 
by the limit for which b= 0. The operator L then reduces for all practical 
purposes to the 3x3 operator s-a of (!) which describes the spin of a spin 
unity particle in its rest system. 

This operator has the eigenvalue zero belonging to the eigenvector a which 
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can be called the non-relativistic spin vector. Generalizing from this obser- 
vation it will now be assumed that a relativistic spin vector can be defined 
such that it is an eigenvector of L belonging to a zero eigenvalue. 

Secondly, it is observed that the energy-momentum vector in the rest 
system has a vanishing space part and can in a formal way be regarded as a 
fourth eigenvector of s-a belonging to a zero eigenvalue. Therefore the second 
assumption is that the relativistic energy-momentum vector of a spin unity 
particle is another eigenvector of L belonging to a zero eigenvalue. 

Since in the non-relativistic limit the spin and energy-momentum vectors 
can be regarded as orthogonal, it will be assumed that this is also the case in 
the relativistic description. Hence the spin pseudovector s, and the energy- 
momentum p, are vectors of the same type as k, and q, and satisfy 


(20) Daf Sto La SuPu= 0. 


From (8) follows that 


(21) Li, = Pu8v— Pr8y 

whence 

(22) &, = Lp, 

and. 

(23) a = SPo— PS, b=p/\s={(P/\a)/p- 


This vector s is the space part of the spin vector s, and must not be confused 
with the matrices s of (1). Note that (22) expresses the spin vector in terms 
of the energy-momentum and the spin operator. 

Within the above relativistic framework it is now possible to study the 
effect of Lorentz transformations. From (23) follows that in the rest system 
of the particle 


(24) i= sea b 0", 


where a° is the spin vector of the particle in this system. Since a and b re- 
present a six-vector it follows from (24) that when the particle has momen- 
tum p, they become 


(25) a=a:+pa,, b=-p/a=p/a, 


Our 0 0 0 0 
a= 4,71 4,, a;|P , a; lp. 


pe 
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Similarly, the spin being a pseudovector is given by 
9 PAP 0 0 0 
(26) S — Poa, # as, So = pa, . 


The consistency of (25) and (26) with (20) and (23) is evident. From (25) 
follows that b vanishes in the rest system while in any other system it depends 
on the spin direction in the rest system and the momentum. | 

In the case of a particle with a magnetic moment the six-vector a, b is 
proportional to the magnetic and electric moment six-vector of the particle. 
This statement follows from (24) and (23). 

The functions s, and m,+ în, are eigenfunctions of L belonging to the 
eigenvalue 0, +1 and can therefore be regarded as the three relativistic spin 
states of the particle. | 


82. Particles with spin 4. — The above procedure for spin unity particles 
can also be applied to spin & particles (?). In this case the Lorentz group is 
represented by the Dirac matrices — i(y,y, — yry,)/4 where y;= — iBa;, ya=f. 
The most general double-valued Lorentz transformation is given by 


y' = exp [1760/2 + T?p/2]y , 


T = (1/2)yuyyLu = — WY My Yo Ny y 
T? = (1/2)yp voli, = Yukuyr ds, TT? =T°T=y,, 
and m,, ny, ky, q are defined as in (1), (6) and (7). 


If 7 is to be the covariant spin operator of the particle it must commute 
with the projection operators 


A, = (1/2)(1 4 ya) 5 


which distinguish between negative and positive energy states. Here p, is 
the energy-momentum vector of the particle normalized as above. Now Hh: 
and A only commute if 


(28) Lospy, = 
If L,, is written as a matrix L and L? defined as above then (28) implies that 


esso LIAS TP Ir 


(2) The relativistic description of electron spin has been approached from another 
angle by H. P. STAPP: Phys. Rev., 103, 425 (1956); Boucrrar and MicHEL, Nuclear 


Physics, 5, 416 (1958). 
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Define s, as in (22) then s, and p, are again eigenvectors of L belonging to 
zero eigenvalues. From 


Lye = ae Ye 
and (27) follows that 


(29) Vu(My ane My) Pe 0. 


Hence (29) defines the four linearly independent eigenstates of 7 in terms of 


m,+in,, the eigenstates of L belonging to the eigenvalues + 1. 
If 7 commutes with A, so does 
(30) Py Pu = WoPu Su = B(C'S — 07384) 


with s, defined by (22). If this pseudovector s, is interpreted as the spin vector 
of the spin 3 particle, as the space part of (30) suggests, then the present si- 
tuation is exactly the same as in (20)-(26) for spin unity particles. The spin 
operator can be expressed as 


T=—o:a-ia:b 


and if the spin } particle has a magnetic moment then a and b are again pro- 
portional to the magnetic and electric moments respectively. 


4. — The case of zero mass particles. 


In this case no rest system exists which can indicate a unique way of 
fitting relativistic spin into the above framework. The vanishing rest mass, 
however, requires that the energy-momentum vector must satisfy p,p, = 0 
which suggests that p, and the spin vector s, be identified with the eigenstates 
vi of L and L”. (17) can be normalized and rewritten as 


Pu = (Po/a)(Nz + ON,, ia), Su = (a/2p9)(N; — DN,, — ia) 


whence 

LiyP» = Pies ie Sy = — Su 
and 
94) Pu Pu = 0 = Susy, Su Pp = 1. 


With this physical interpretation the relations (21) and (23) between the spin 
and momentum four-vectors are still valid for zero mass particles. However, 
(22) is no longer valid. 


où 
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The spin + and 1 particles for which (31) would be valid would have some 
kind of inner structure described by the six-vector a, b occurring in the spin 
operators L or T which are still 4x4 matrices. At present none of the known 
particles appear to be of this type. 


5. — The photon and the neutrino. 


Apart from the description of the spin of zero mass particles presented in 
Sect. 4, there is another description which is suggested by the finite mass case 
when the mass +0. An equivalent procedure is to allow P>D > © in 
which case (26) becomes simply 


(32) «AP 
where À is an undetermined constant and 


(33) Su Su = PuPp = SuPy =. 


Now the number of eigenstates of the spin operator has been reduced to three 
and L”, as defined by (21), no longer serves a useful purpose. Consequently 
the spin operator for the spin unity particle becomes the 3x3 matrix s- P 
of (1) with the eigenstates p (the momentum) and m + in (the spin states), 
where, by (1) and (6), 


m'p=n"p=mn—=0, m'm=nn—l. 


For the spin ? particle the spin operator becomes the matrix o-p — 
=—1Us-m)(o-n) with the eigenstates y, given by 


c'(m+in)y, = 0. 


In the absence of parity conservation for such a particle, o can be regarded 
as à 2x2 matrix. 

The most general 2x2 matrix with determinant unity can be written in 
the form exp [io-(a + ib)(0 + ig)] with a and b satisfying (14). Since these 
matrices serve as representation for the Lorentz group, one might expect the 
Pauli matrix o-(a+7b) rather than o-p to be the spin operator for a spin à 
particle that satisfies a two component theory. However this is not the case 
since the existence of a six vector implies the presence of an inner structure 
which, according to the above discussion, is only possible in the four com- 


ponent theory. 
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Où 
en 
ho 


Since both the proton and the neutrino have zero mass and each has two 
spin states which appear to describe their inner structure completely, it must 


be concluded that their spins are described by (32) and (33) rather than by (31). 


This work was completed with the assistance of the Carnegie Corporation 
of New York and the South African C.S.I.R. while the author was enjoying 
the hospitality of Columbia University. 


Note added in proof. 


The discussion in Sect. 4 does not take into account invariance under inversion 
of axes. Such invariance requires s,p,= 0 which, together with 


SyMy = Sn = Pully = Pully = PP 0.5 


yields 


8, = APu and 88h 0 - 
Hence in this case the spin of a zero mass particle is described completely in terms 
of the energy-momentum vector and two linearly independent spin states as outlined 


in Sect. 5. 


RLASS.U NTO (4) 


Si genera la più generale trasformazione propria di Lorentz finita in termini di 
una trasformazione infinitesimale. Tale generatore si interpreta come operatore rela- 
tivistico di spin capace di descrivere in modo relativisticamente covariante lo spin 
e le proprietà elettromagnetiche delle particelle di spin 3 e 1 e massa finita. Si indica 
l'estensione del metodo a particelle di massa a riposo evanescente. 


(*) Traduzione a cura della Redazione. 
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On Sommerfeld’s Approximation in High Energy Photoelectric Effect 
and one Quantum Annihilation of Positrons in the K-shell. 


H. BANERJEE 


Department of Physics, Indian Institute of Technology, Kharagpur - India 


(ricevuto il 26 Settembre 1958) 


Summary. — The angular distribution and the total cross-section for high 
energy photoelectric effect and one quantum annihilation of positrons in the 
K-shell have been calculated with the use of the Sommerfeld-Maue approx- 
imation for the wave function of the scattering state of the fermion. The 
error in this procedure is of relative order (a/e), where a= Z/137 and e is the 
energy of the scattering state in units of mc?, and is, therefore, negligible 
in the extreme relativistic case. Results have been obtained correctly up to 
terms of relative order a and methods for calculating the contributions 
from higher terms have also been indicated. The expression for the angular 
distribution of photoelectrons obtained in this paper does not show any 
marked qualitative difference from Sauter’s formula. One may conclude 
from our results that the ratio of the forward to the maximum intensities 
of the photoelectrons diminishes as the primary energy is increased. The 
absorption coefficient for photoelectric effect has been presented in a form 
which is similar to Hall’s results. Correction terms in the formulae of 
Bhabha and Hulme for the positron annihilation process have been found 
to be the same as those in Sauter’s formula. 


1. — Introduction and general discussion. 


Although the photoelectric phenomenon is one of the basic physical pro- 
cesses of quantum electrodynamics, the only theoretical result for the angular 
distribution of photoelectrons that is available to experimental workers for 
ready reference is that of SAUTER (1). SOMMERFELD (?) has also obtained the 


(1) F. SAUTER: Ann. Phys., 11, 454 (1931). 
(2) A. SOMMERFELD: Wellenmechanik (New York, 1939), pp. 482-494. 
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same result following a different method of approximation. These results have 
been obtained only after neglecting terms of relative order a (= Z/137) which, 
however, is not too small for elements that are generally used as converters 
in experiments on photoelectric effect. The inadequacy of theoretical results 
which has prevailed since 1935, naturally makes one enquire whether with 
any one of the available methods of approximation an improved result in a 
closed form for the angular distribution can be obtained and thereby the 
discrepancies between Sauter’s result and those recently reported by exper- 
imental workers (*) can be explained, at least partially. 

The Born approximation method, which is commonly used in high energy 
scattering processes, is, however, not suitable for the photoelectric pheno- 
menon. Although it gives the correct Z-dependence, the first Born approx- 
imation, which consists in representing the state of the scattered electron by 
a plane wave, introduces an error even in the approximation up to which Sauter’s 
result is valid. This will be apparent from the calculations in this paper 
(Sect. 3). The reason for the failure of the first Born approximation is that 
important regions for the integrals occurring in the matrix element are at 
distances of the order of #/mc, the Compton wavelength of the electron, from 
the nucleus. Consequently the wavefunction of the electron, which is ap- 
preciably distorted by the nuclear Coulomb field within this region, cannot 
be approximated by plane waves without vitiating the results. In fact, in order 
to obtain Sauter’s results by this method, one has to take into account the 
contributions from the next higher approximation. So that, for any impro- 
vement on Sauter’s results the third Born approximation in the nuclear Coulomb 
field has to be considered, which makes the calculation prohibitive. 

If wave functions separated in polar co-ordinates are employed for the final 
state, as was done by SAUTER (1), one is confronted with the difficult task of 
summing an infinite series whose different terms correspond to the angular 
and magnetic quantum numbers J and m respectively of the final state, in 
addition to integration difficulties. In the derivation of the absorption coef- 
ficient in the extreme relativistic case, HALL (*) succeded in summing this 
series only after neglecting terms of O(1/l). He has, however, expressed the 
absorption coefficient in the final steps of his calculation in terms of an integral 
which is not integrable in a closed form. For the angular distribution the situation 
further worsens and the only available result is that of SAUTER (!). A better 
approximation of the angular distribution by this method seems to be very 
difficult, even if some suitable simplification is found possible. 

So far as the photoelectric effect is concerned, Sommerfeld’s (?) method of 


(3) A. HEDGRAM and S. HULTBERG: Phys. Rev., 94, 498 (1954); S. HULTBERG: 
Ark. f. Fys., 9, 245 (1955). 


(4) H. HALL: Revs. Mod. Phys., 8, 358 (1936). 
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treatment has received the least attention. His method essentially consists 
in representing the state of the scattered electron by a wave function (5) 
known as the Sommerfeld-Maue wave function, of the form 


F) 


(1) Y= x, (1 


a 


)ur, 


where N, is the normalization constant, (p, e) the momentum and energy of 
the electron, F is a confluent hypergeometric function, v is the normalized 
Dirac spinor amplitude for an electron with momentum p and a= (— yy) 
is the Dirac operator. With the above approximation for the final state and 
neglecting relativistic effects for the bound state, which amounted to Pauli 
approximation, SOMMERFELD was able to obtain Sauter’s result. Sommerfeld’s 
calculation in its original form could not be extended to obtain a better result ; 
firstly, because the order of accuracy of the wave function (1) was not fully 
realized, and secondly, because of the Pauli approximation for the bound 
state one is forced to neglect terms of relative order a in the matrix element 

It is only recently (*7) that much attention has been focussed on the Sommer- 
feld-Maue wave function. Separating the wave function (1) in polar co-ordinates, 
BETHE and MAXIMON (°) have shown that the term of the resulting series cor- 
responding to the angular momentum / coincides with the exact continuum Dirac 
wave function in spherical co-ordinates if in the latter terms of relative order 
(a?/l) are neglected. As will be shown in Appendix C, in high energy photo- 
effect, the main contribution to the integrals in the matrix element comes 
from terms for which 1!~e«. Further, it will be apparent from the calculations 
in Sect. 3 that an error of O(a?) in the wave functions will lead to an error 
of O(a) in the matrix element. Consequently the error committed in the matrix 
element by writing (1) for the wave function of the scattered electron will be 
of relative O(a/e). The Pauli approximation of the bound state, in which terms 
of relative O(a?) are neglected, is, however, not necessitated by the above 
approximation for the scattering state. In fact, in the high energy case one 
may represent the final state of the electron by a Sommerfeld-Maue wave func- 
tion and still treat the bound state exactly without giving rise to any incon- 
sistency. This is exactly the argument on which Hall’s calculations for the 
photoelectric effect in the relativistic region have been based. 

In the present work an attempt has been made towards treating the wave 
function of the K-electron exactly, while the state of the photoelectron has 


(5) A. SOMMERFELD and A. W. MAUE: Ann. Phys., 22, 629 (1935). 
(6) H. A. Berne and L. C. Maxtmon: Phys. Rev., 93, 768 (1954). 

) H. OrsEn, L. C. Maximon and H. WERGELAND: Phys. Rev., 106, 27 (1957); 
Horton and E. PHiBBs: Phys. Rev., 96, 1066 (1954). 
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d-Maue type of wave function. The error 
h a procedure, is of relative O(a/e). In- 


tegrals occurring in the matrix element do not appear to be integrable in a 
closed form. They have, however, been reduced to an elegant form from which 
successive approximations in powers of a may be possible. For our purpose, 
contributions to the matrix element only from terms of relative O(a) have 
It has been found that the correction term in the angular 
expression for 


been approximated by a Sommerfel 
committed, which is inherent in suc 


been obtained. 
distribution leaves Sauter’s result qualitatively unaltered. The 
the absorption coefficient has been rendered in such a form which clearly brings 
out the similarity with Halls result. The calculations for the photoelectric 
effect have been extended to the case of one quantum annihilation of a positron 
in the K-shell in Sect. 5. For the latter process correction terms of relative 
O(a) in the only available results of BHABHA and HULME (5) have been obtained. 
In the high energy limit the error in our results is of relative O(a’). Methods 
have, however, been indicated how our results could be further improved by 
finding higher order terms. In the light of this paper it will also be evident 
how easily the polarization effects of the photoelectrons, which will form the 
subject matter of a subsequent paper may be included in our consideration. 

Although the present work has been carried out with Coulomb wave fune- 
tions both for the initial and the final states of the electron, deviations due 
to the screening effects of internal and external charge distributions will not 
materially affect our results. Within distances from the nucleus of the order 
of the Compton wave length of the electron (see Appendix C), where the be- 
haviour of our wave functions is important, the potential function V(r) may 
be approximated with sufficient accuracy by 


according to the Hartree self-consistent field. For heavy atoms the « inner 
screeing constant » is negligible, while the « outer screening constant », the 
only effect of which is to change the energy of the photoelectron by Yj, may 
be neglected in the high energy case. 3 


2. — Matrix element. 
Introducing m as the unit of mass, #/me as the unit of length, and me? 


that of energy, one can write the matrix element for the photoelectric effect 


(8) H. J. BnaBna and H. R. HuLmr: Proc. Roy. Soc., A146, 723 (1934). 
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e PISA 
H,=<==|['P;(e-YyPexp[ik-r]dr 
1 = yey | Male WM exp Like rar, 
where e is the polarization vector of the incident photon, and (k, k) are its 
momentum and energy. The energy of the photoelectron is given by Hinstein’s 
equation which may be written as 


(3) e=k+V1-aak+1, 


if terms of O(a?/e) are neglected. Y,, the wave function of the K-electron, is 
given by 


(4) ee 4 Ni Pr | ye exp [— ar] W , 
where 

———. D, Nei 
EROE dai zar SE (pd ea 
(4a) EVENT: LAN, = (a) re =) 


The normalized spinor w represents the polarization state of the K-electron. 
VY, is the Sommerfeld-Maue wave function as in (1). The normalization con- 
stant N, and the confluent hypergeometric function F are to be so chosen that 
asymptotically (5) 


| 


exp [— ipr] 
pe fel 


nan a ju exp [ip:r] + 


This condition may be satisfied if one chooses 


(5a) F = F{—ia,; 1; —ipr—ip-r}, 
and 

ei [Ti F 
(50) N = Cola exp È M1 + dia), 


where a, = ae/p. Substituting (1) and (4) in the expression (2) for the matrix 
element we get 


(6) Ay, = Cu [(e-y)Lo + (e*y)(Y: I) + (a: L)(e da I,|w , 


(9) J. M. Jaucx and F. RonrricH: The Theory of Photons and Electrons (Cambridge 
Mass., 1955), pp. 318-326. 
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where 


NME 


= 
Gi 


+ 7 


In the above we have introduced the following abbreviations, 


| 


(8a) Le fexp [iqr] F*r-° exp (—ar]dt ===. Le 0). 


ca 
ia 
(8b) l= sj fe [ig:r] pei si -§ exp[— ar]dt = — DES, V,J(— 9), 


(8c) = 3 fe [iqg:r](V, F*) r ° exp [— ar|dt = 


2e 
er *\ 8 ex Le Bag ies 
= — 5 Jespria- ri] Lv,F* exp[— arjar = — = (VJ(—9)) ; 


dE 


ia ; 43 os ; 
(8d) Je so | liq:r] (a: V,F*)(e-y) (y-7) r-exp[—ar]dt, 
where q= k—p is the momentum transferred to the nucleus and 
E dt 
(9) J(— 6) =|exp [ig:r] F*r ? exp [— ar] aa 


Contributions to the matrix element from I, have been shown to be of rela- 
tive O(a/e) in Appendix B and can, therefore, be neglected. 


3. — Evaluation of integrals. 


A convenient contour integral representation for the confluent hyper-geo- 
metric function F* in the integral J(— 6) in (7) is 


(2+,1+) 


; È 7 é 1 
(0) F*= Flia,; 1; ipr + ip-r) = sat dt ts-1(t—1)-* exp [i(pr + pri]. 


The contour encircles the point 0 and 1 in the t-plane once in the anticlock- 
wise sense. At the point where the contour crosses the real axis to the right 
of 1 both argt and arg ({—1) are zero. We substitute the above represen- 
tation of F* in the expression for J(— 6). Uniform convergence in ¢ of the 
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space integral requires that the contour should not contain any singularity 
fo the space integral, which condition may be easily satisfied. With this under- 
standing one can change the order of integration, perform the space integration 
first, and obtain in an elementary way 

(0+,1+) 


(11) HT) dat 1)-i4- 


{(a = ipt + ig)? — (a — tpt — ip 
CT ES à : 


(12) x = Vep + 2ip:q +4. 


The possible singularities of the integrand in (11) are the branch points of x 
and the zero of {(a —ipt)?+ x}: One can show easily that the integrand is 
regular at the branch points of y. The only singu- 

larity is thus at t= o, where 


(13) o =— (a? + q°)/2(q°p — ipa). 


Let us now extend the contour to a large circle 
T of radius 9 excluding the branch point at t=o 
by a small circle y of radius 7 as shown in Fig. 1. 
We have then 
(14) J(—6) = dit It It In: 


Fig. 1. — The contour in 


where the respective contours are denoted by the the pese 


subseripts and have been shown in Fie’ 1. ‘The 

lines L,, L, have been chosen to extend from o to the left so that they may 
not contain the branch points of y. It can be shown in an elementary way 
that as o > co and 7 — 0 both J, and J, vanish. We fix the branch of x 
in such a way that y —i(a — ipt) at t= 0. We further define 


: we. La (a — ipt — y) 
(15a) J = ra ofart CNE (a pe + yt?’ 
and 
| ee pt FAN 
i ei eee) at (s(t Loe — — si: 
(150) J ( ) ( ) PAU x ipt)? + Th tas 


The integrand of J" is regular at t= 0. Therefore, 
Jit J, = 9. 


56 - Il Nuovo Cimento. 


3831 


870 H. BANEJREE 


We choose arg (t— a) to be x on L, and — x on L,. After considerable sim- 
plification we then have. from (14) and (15) 


(16) JC 6) 0 [ea = ia) her 


(17) 0=I1-d)(2exp[-t2i)), 


co 


- (exp [276] — 1) Ree eae = tie ee 
as) x= CRESS [aga — aged — 98 BOE 
and 

An(a? + g?)ia- 
ai 5 {(a — ip) + (p + ql} 
In the above 
| ia po + ia 
== i = 25 — = 
eta p(o—1) 


A ee ’ po + p'q 
= AE) — . 
| Rive 2Ex +1 with œ eae er Ty 


It may be observed that the integral in the expression for K in (18) be- 
haves as 2/0 for small 6. The presence of the factor (exp [2710] — 1) however, 
ensures that the value of K is finite. The expression for K in (18) may be 
further simplified if one remembers that in our formulation we can always 
neglect terms of relative O(a/e). For this purpose, we observed that in high 
energy photoeffect the angle of emission 6(= <x È, p) for most of the photo- 
electrons is of O(1/e). Therefore, q is of 0(1) and v, as defined in (20), is of 
order 1/e. The integral K exists if v= 0, so that, if it be expanded in powers 
of v in the neighbourhood of v= 0, only positive powers of v will appear in 
the expansion and, that too, in terms of O(a’). We may, therefore, put v — 0 


in the expression for K, the error committed thereby being of O(a/e). We 
thus have 


Le) 


_ (exp[2xiô] — 1) i we 
(21) K n | dé(1 — ué)in: I 


0 


It can be easily shown that 


G2) K = 1 + O(at). 
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Neglecting terms of O(a?) we can therefore write from (16) 


4r(a? SL qa 
23 I(—8) & OF, = 073 LE 
( ) ( ) {(a — ip)? + (p + q)?}i™ 


This is exactly the value obtained by SOMMERFELD (2). Substituting this 
approximate value of J(— 6) in eqs. (8) we obtain after some simplifications 
and neglecting terms of relative O(a/e) and O(a?) 


; C'Ja(1— ia) f di i \ 

(24a) ih = (e 2 1) La —B ne 0) a , 
il'Ja(1— ia) | 1 E) 
(245) E = Ja — 19)€=1) led — A c08 0) take ; 
6 O' J ,a(1 — ia) 

(240) ley 


where f is the electron velocity and 0 =< (k, p) is the angle of scattering. The 
correction terms for I, and I, in (24a) and (24b) due to the approximation 
of J(-d) by J are of relative O(a) and O(a?) respectively. The integral I, 
is already of O(1/e) relative to I and I, so that its correction terms will 
contribute only terms of relative O(a/e) in the matrix element, which we shall 
neglect. It is thus evident, up to terms of relative O(a), only the expression 
for I, in (24a) will be affected. 

In eqs. (24) terms of relative O(a?) have been intentionally retained in order 
to preserve symmetry. It may be pointed out here that even in the approx- 
imation up to which eqs. (24) are valid, the Born approximation would not 
give the same results as ours. In fact, in the Born approximation the integral 
I, and the second terms within the braces in (24a) and (245) do not arise. 

The problem of obtaining a suitable expression for the integral. J(— à) 
in (16) thus reduces to evaluating K in a closed form, which, however, does 
not appear to be possible. Derivation of a closed expression for K correctly 
up to terms even of O(a?) will prove to be a difficult task. But so far as our 
matrix element is concerned we need find only the derivatives of K with 
respect to the parameters w and # in terms of which it has been defined. It 
will be presently shown that these derivatives, though not exactly integrable, 
admit of suitable approximations. 

From (20) and (21) we have, 


[ce] 


€ x 270 — i 6 — + R 1-ò 
(25a) 5 = (exp a ä 1) (1 4 iu) fasti dci né}? ( wi LIO i 


0 


wo 
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= 1-6 doh 
aK _ (exp[2xiô] —1) tr DE OS È È A 
+= SE fase Sti is BTS oe 


0 


(25d) 


The integrals in (25) are now convergent as 6 + 0 and the factor (exp [2710] — 1) 
ensures that oK/déw and 0£/0x are both of O(a?). We may, therefore, put 
6 = 0 in the integrands if terms of relative order a* are neglected in the matrix 
element and obtain 


foe) 


0K Gé" TE 


(26a) as Dai — ia, JE — yé)in-? = + O(a‘), 
and 

é IE 
(260) = =; [asa Te agp EE) + O(a4). 


0 


The values of 0K/0u and cK/cx in the first approximation, as derived in Ap- 
pendix B, are 


DEMO Tri a? TiPq 
27 sh ee sl ~F(- hax aston 
FF ou ii + Vato) 4 2p:q 
and 

x Ginn Owed, ie 2 
ee) de 21-% 8 \p'q 
Therefore 
(28) ok =F (se). 
da 4\p-q q 


The expression for J,, correct up to terms of relative O(a), may now be easily 
calculated from (8a), (16) and (28) 


; bar's, ts ET Wl eS FANO] 
‘(Ca ee ir e a 


Substituting the values of L, I,, and I, from (24b), (24c), and (29) in the 
matrix element H,, in (6) one now gets the result correctly up to terms of 
relative O(a). For a better approximation the values of 2X/0v and OK | 0x, 
and hence of I, correctly up to terms of relative O(a?) have to be calculated 
from eqs. (26). For the evaluation of the correction terms of I, the values 
of 0K/èdu and 0K /cx as given in eqs. (27) are sufficient. 
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4. — Angular distribution and absorption coefficient. 


Summing over the final and averaging over the initial spin states of the 
electron one now gets from (6) 


(30) | His|? = 


role 


l@|2((e—1) | Lol? + (e + DL —iB|? + 
4 4(e + 1) Im (e-I)(e-I?) + 4(p-e) Im J, (ef) + 21m {p-(L — iL)}]. 


In the above the symbol Im signifies that we are to take the imaginary part 
of the expression following. Substituting the expressions for ioe le onda 
derived in the preceding section and neglecting terms of relative O(a/e) and 
O(a?) we get after some simplifications 


igh te OR 
2(1 — B cos 0)? 


(31) | Hy» |? ae |M} 


TELE pmi et (e | mafe gina” 
ET sa —B 0080) 2(e — 1) cos! Gr — 5 — B cos 0} os 
where 

Li . , 27°a?(1 + a?) ia è a 
Go) para OO Pi {— ora + 20 ty 2: 


and g is defined by p'e— P sin 0 cos ¢@. 

We observe that the result given in (31) does not differ qualitatively from 
Sauter’s formula for the angular distribution. It is apparent from eq. (31) that 
up to our order of approximation there will be no emission of photoelectrons 
in the forward direction (0 = 0). The experimental value 0.29 of HULTBERG 
and HEDGRAM (°), the exact theoretical value 0.27 and the first approximation 
value 0.21 obtained by SAUTER and WUSTER (1°) for the ratio of the forward 
and maximum intensities of the photoelectrons ejected from Pb-converters 
by 1.3 MeV y-rays show that most of the contributions to the forward distri- 
bution of photoelectrons come from terms of relative order a/e in the matrix 
element. We are, therefore, led to conclude that the ratio of forward and 
maximum intensities of photoelectrons should diminish as the primary energy 


increases. 
The absorption coefficient of the K-shell is given by (°) 


(33) tp 2x (a) | | Hy, |?d2 , 


= 


(19) F. SAUTER and H. O. WustER: Zeits. f. Phys., 141, 83 (1955). 
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where N, is the number of atoms per cm*. Substitution of the expression for 
|H,, ? from (31) then leads to 


exp [— xa + 2a tg! (a/1 — 6)](1 + ((83 — 27*)/6) è) 


(34) vote N aa lin 


3 Anais», {| fra 2(e-2) sere, GE) 
(e —1)2 | 15 3e? Sechi 
approximately, where « is the fine structure constant. Although our result 
is correct only up to terms of relative O(a), we have retained higher order terms 
which appear naturally in our calculation so as to illustrate the close similarity 
of our result with the expression derived by HALL (4) for the absorption coef- 
ficient in the extreme relativistic case. It may, however, be pointed out that 
the a-dependent term within the braces in (34) viz., (— 4774/15) differs from 
the corresponding term 0.05 tg-! (.63/a) in the expression given at the end 
of Section 21 of reference (*). The difference may be attributed to the pro- 
bable inaccuracy of the approximation of a very complicated double integral, 
which Hall was forced to make in the final steps of his calculation. In fact, 
from the way this approximation has been made by HALL it appears that the 
error in Hall’s result is of relative order a? if not greater. 

In view of the close agreement of the results of the exact numerical calculation 
of HULME et al. (11) with the formula of Hall, the validity of our correction 
term may be questioned. One may, however, argue that in the particular for- 
mula which is given in eq. (17) of reference (4) and to which HULME et al. have 
referred, the term 0.05a tg-! (.63/a) does not occur. Secondly, although the 
formula of Hall was modified by HULME et al. so that it might coincide with 
Sauter’s formula in the limit a — 0, it is, however, clear that in any such 
modification contributions from terms of relative order a/e (which is not small 
in the range of energy considered by HULME et al.) can not be included. It 
is quite likely that these contributions would partially cancel the deviation 
of Hall’s formula from the correct one. For further justification of our re- 
sult (34) it may be pointed out that the diminishing effect of the term 
(— 474/15) would be sufficiently neutralized by higher order correction terms. 
This is apparent from the factor (1 + ((33 —27°)/6) 6) in (34) and also from 
the fact that in the case of large Z (viz., Z = 84), where these correction terms 
are significant, numerical results obtained from the formula (17) of reference (+) 
are less than the corresponding results of HULME et al. (11), whereas for low Z 
the contrary is true. 


(1) H. R. HuLME, J. McDouGaALL, R. A. BUCKINGHAM and R. H. Fow er: Proc. 
Roy. Soc., A149, 131 (1935). 
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5. — Cne quantum annihilation of positrons. 


The process of simultaneous annihilation of a positron and a K-electron 
with the emission of one quantum of radiation may be considered as the tran- 
sition of the bound electron to a state of energy — £ and momentum — p, 
where (p, €) are the momentum and energy of the incident positron. One can 
thus directly obtain the matrix element for this process if in the matrix ele- 
ment (2) for the photoelectric effect (k, k) and (p, e) are replaced by (— k, —k) 
and (—p, — e) respectively, according to the well-known substitution law. 
Calculating exactly by the method followed in the preceding sections one gets 


rada a8Z5 grein? Do. 
8a2e2(e + 1) |2(1— Boos 0)? 


4 cos? p el ee mafe MAN 
= # 0080) 2(e + 1) cos°pj— alt feos Ot 2 


’ 


sya 2 
pee) 


of the angular distribution of the polarized photon and 


€ Ana  2(e +2) Es +2 5) 
5 in | 5. : ae FASE LE colli 
A ia pt 38 Bet nl; TB)’ 
with 
cp [— ma + 2a tg (a/1 — d)) 33 —27 
(37) ri na + 2a tg (a/l IDÉES — 276 s), 
a2 6 


for the annihilation cross-section due to the two electrons in the K -shell. The 
undefined symbols used in the above exactly correspond to those employed 
in the case of the photoelectric effect. 

From eqs. (35) and (36) it is apparent that the expressions for both the 
angular distribution and the total cross-section can be obtained directly from 
the corresponding formulae for the photoelectric effect by applying the sub- 
stitution law except for a factor (e+1)?/B°e? in o,. The origin of this factor 
may be traced to the fact that while obtaining the total cross-section on 
positron annihilation, one has to integrate the expression for the angular dis- 
tribution over the momenta of the emitted radiation and divide by the velo- 
city f of the incident positron, whereas in the photoelectric effect the inte- 
gration is over the momenta of the photoelectron and the velocity of the 
incident quantum is unity. One may further observe that the correction terms 
in the formulae of Bhabha and Hulme (*) are the same as those in Sauter’s 
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formulae. Our results thus resolve the doubt (1?) about the applicability of 
the correction terms for the photoelectric effect to the positron annihilation 
process without any modification. 


* OK *% 


The author wishes to thank Prof. S. Gupra and Dr. G. BANDYOPADHYAY 
for many helpful discussions and encouragement throughout the present 
investigation. I am also grateful to Prof. J. M. JAUCH of Iowa State 
University, U.S.A., for his kind interest in this work. 


APP iN Dux, 


The eqs. (27) will now be derived from (26) employing a proper limiting 
procedure. The difficulty that arises in our approximation is due to the sin- 
gularity of the integrands in (26) at £ — 1/u which, according to eq. (21), 
approaches = 1 as a — 0. 

The derivation of (27b) presents no difficulty as the integrand behaves 
properly at £ = 1. One readily obtains setting a = 0 in the integrand in (26b) 


(A.1) = =-{—=-- 


From (26a) one gets 


Se (co) 


j iL SY ta,—1 
+ Jar + (1 w) fast — w&)iu-2 “a 


0 0 


It should be observed that the integrals in (A.2) vanish identically if one 
directly sets # — 1 and a = 0 in the integrands. But, as will be shown pre- 
sently, the first integral in (A.2) has a non-vanishing value in the limit a — 0. 
We write it as the sum of three integrals 


i ca LEM ee Eth an Be 
(1 — w&)ia-1 n (1 — ug)iani 
i leg <f+f+ [Jase 
0 0 ig T4 


(12) J. C. JAEGER and H. R. HULME: Proc. Camb. Phil. Soc., 32, 158 (1936). 
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where 7 is a small positive quantity, but otherwise quite arbitrary. We take 
n> \1—w| and a. The singularity £ = 1/u of the integrand lies outside the 
range of integration for the first and the third 
integral. One can, therefore, set u—1 and a=0 
in these integrands and ultimately find that they 


cancel each other. Such a procedure is not per- o) ere 1 
missible for the second integral because the sin- fi / à 
gularity shifts from 1/u to 1 in the é-plane 

as a->0 and thus lies within the range of in- Fig. 2. — Contour in the 
tegration. From (21a) it is apparent that the £-plane. 


point 1/w lies below the real axis in the É-plane. 

We can, therefore, replace the path (1 — n) 

to (1 +7) by a semi circle of radius 7 about £ — 1 above the’ real axis 
(Fig. 2). Further, as R is analytic in the neighbourhood of È = 1, one can - 
replace R by V2(1 — x)(1 + 6’), where |0'|, which is small, vanishes as 7 — 0, 
and this factor may be brought outside the sign of integration. So that 


Lt Jas ee Be 


0 


al i 
= Lt fi exp [10] d0(1 — u — un exp [10])!! = 
V2A4:="x)(1-4- 0") 70 È Lee] 


ati 


~ /9(1— w)(1 +0) 


We can now make 7 — 0 and finally get 


n+1 
(A — we) iar i 
L d aversi 
ee 2] - R /2(1 — 2) 
1-7 


In a similar manner one can show that the second integral in (A.2) tends 
to zero with a. Neglecting terms of relative order a we, therefore, get 


OK | a xi 
(A.5) Sediari 


APPENDIX B 
The term I;, neglected in the matrix element, can be written from eqn. (6d) 
in the form, 


ap 
(A.6) I, = — 2e(2 — 6) (a: V,)(e-y)(y: Var ’ 
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where J, may be obtained from J(— 06) by replacing lo) by (1 +6). One, 
therefore, gets from (16) and (21) neglecting terms of relative O(a?/e) 


I'(— 6)(exp [216] — 1 


ree ) (2 exp[— 4 wi]) °K, , 


(A.7) J, = 


where 
pidge a tre on M 
K, = (24 — ia) [aga ue 
0 


It follows from (20) that w and PR and, therefore, K, are of zeroth order 
in p. So that (V,K,) can at most have terms of O(1/e). Thus, if we neglect 
terms of O(a/e) in (V,X,) and, therefore, of O(a) in K, we get from Appendix A 


ELA TS I 
à p V2 x) 2q 
so that 


(V,K,) = 0. 


Therefore, (V,A,) is of O(a/e). Further, from (24c) it is apparent that (V,J,) 
is of O(a/e). Thus, from (A.6) and (A.7) it follows that the contribution of J, 
to the matrix element is of relative order a/e and can, therefore, be neglected. 


APPENDIX C 


i The following calculations are intended to show that the maximum con- 
tribution to the integral J(— 6) arises from: (a) regions within distances from 
the nucleus of the order of Compton wave-length (i.e., r & 1) of the electron, 
and (b) final state wave functions which correspond to angular momenta la e. 


pen’ a) Separating the confluent hypergeometric function in the integrand 
in (7) in polar coordinates (1%) one obtains, 


(A.8) exp[— ipr] F(ia,;1; ipr + ipr) = mon > (21+1)(—i)L(r)P,(cos 0), 


where 


Lt) = [MU + 1 — tay) (21 + 1)1] (2pr)* exp Lipr] FU + 1 — ia, ; 21 +2; — 2ipr), 


(°) W. GORDON: Zeits. f. Phys., 48, 180 (1928). 
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and 9 = x(p,r). Substituting (A.8) in (7) one obtains after performing the 
angular integrations 


oct I(l+1—ia,) 


49 JD Tr a? CUTS 


(pk)'P,(cos 0): 


œ 


-|drr2*2-° exp[— ar —i(k— p)r] F(1+1—ia; ; 21+2;—2ipr) F(1+1; 214-2; 2ikr). 


Maximum contributions to the integrals thus come from the region 


1 
A.10 Fe 5 = = = wil n 
pet) "7 Ta Lik —p)| 


in the high energy limit. For distance larger than this the exponential factor 
falls off rapidly. 


b) In order to find the relative order of contributions to J(— 6) from 
final state wave-functions corresponding to different J, it is enough to obtain 
an approximate expansion of the from 


(A.11) J(— 6) = Y M()) P.(cos 0), 


and consider the behaviour of M(l). Neglecting terms of relative order a we 
obtain (11) from (19) 


2x [1 En “ail 
WI "Di (; — cos 0) SPL, (3) Pcos 0), 
U 
where Q, is a Legendre function of the second kind of order J. For small values 
of J and as B >1, the behaviour of Q, is approximately given by (aay 


\ 


i var 
(A.13) (3) = ju a) SIAT), 


in terms of the well-nown Ÿ-function and the Euler constant y. From 
(A.11) and (A.13) we therefore conclude that for small values of 1, M (1) increases 


(14) A. ERDELYI: Higher Transcendental Functions, Vol. 1 (New York, 1953). 
pp. 168. 163. 
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as (21 + 1) In 2e. If, however, Lis of O(e), so that (I/p) = @ is finite, we get (7°) 


1 x \i 
(A.14) ol) (ent, 


asymptotically. Hence in this region M, increases as (2lp)t exp[— 1/p], which 
is maximum for / = p. 


(5) H. Jerrreys and B. S. JEFFREYS: Methods of Mathematical Physics (Cam- 
bridge, 1956), p. 656. 


RIASSUNTO (*) 


Si sono calcolate la distribuzione angolare e la sezione d’urto totale per l’effetto 
fotoelettrico alle alte energie e l’annichilamento monoquantico dei positroni nello strato A 
ricorrendo all’approssimazione di Sommerfeld-Maue per la funzione d’onda dello stato 
di scattering del fermione. L’errore commesso in questo procedimento è dell’ordine a/e, 
dove a= Z/137 ed e è l'energia dello stato di scattering in unità me? ed è pertanto 
trascurabile nei casi relativistici estremi. Si sono ottenuti risultati corretti fino a termini 
di ordine relativo a e si sono anche indicati metodi per calcolare i contributi di termini 
più elevati. L'espressione per la distribuzione angolare dei fotoelettroni ottenuta nel 
presente lavoro non mostra differenze qualitative marcate dalla formula di Sauter. 
Si può concludere dai nostri risultati che il rapporto dell’intensità equiversa a quelle 
massime dei fotoelettroni diminuisce col crescere dell’energia primaria. Si presenta 
il coefficiente d’assorbimento per l’effetto fotoelettrico in forma simile ai risultati di 
Hall. Si è trovato che i termini di correzione per le formule di Bhabha e Hulme per 
il processo di annichilamento del positrone sono uguali a quelli della formula di Sauter. 


(*) Traduzione a cura della Redazione, 


3842 


IL NUOVO CIMENTO VOLTA EN. 5 1° Dicembre 1958 


High Z Nuclei in Cosmic Radiation. 


V. Bisi, R. Caster, C. M. GARELLI and L. TALLONE 


Istituto di Fisica del Universita - Torino 
Istituto Nazionale di Fisica Nucleare - Sezione di Torino 


(ricevuto il 3 Ottobre 1958) 


Summary. — Adding the results obtained by measuring the charges of 
117 heavy nuclei (Z > 10) to those obtained in a previous paper, we deter- 
mined the values of the ratio N(Z >15)/N(Z > 10) at flight altitude and 
at the top of the atmosphere, and of the ratio N(Z—26)/N(Z > 10) at 
fight altitude, at the top of the atmosphere and at the C.R. source. 
It results that the chemical composition of the C.R. source differs from | 
the average composition of the universe, being much richer in heavy 
nuclei (Z >15) and in particular in Iron. The absorption mean free path 
in air was also determined studying the zenith angle dependence of C.R. flux. 


1. — Introduction. 


The knowledge of the detailed features of the Cosmic Ray charge spectrum, 
is of great interest in the problem of the Cosmic Ray origin. 

In particular from the relative abundance of nuclei of charge Z >10 
(H-nuclei) information can be derived on the chemical composition of the 
C.R. source, since this part of the spectrum is qualitatively the less affected 
in the diffusion through space and, therefore, may reflect its original com- 
position. 

Unfortunately this part of the spectrum is the less well known owing both 
to lack of statistics and to experimental difticulty in obtaining a good charge 
resolution. In a previous paper (1) which will be referred to as paper I, the 
Turin group used for charge analysis of the heavy primary spectrum a photo- 


(1) R. CESTER, A. DeBENEDETTI, C. M. GARELLI, B. Quassiati, L. TALLONE and 
M. Vigone: Nuovo Cimento, 7, 371 (1958). 
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metric device (2) that allowed a good charge discrimination up to the heaviest 
C.R. nuclei (Z = 26). 

The present research was carried out, under experimental conditions com- 
parable with those of (I), in order to improve statistics for H-nuclei. 

The interaction mean free path in air for these particles was also obtained 
studying the flux dependence from the zenith angle. 


2. — Exposure details. 


The emulsion plates used in this research were part of a large G-5 stack 
flown with vertical geometry over Texas (A = 41° N) in January 1955, at an 
atmospheric depth of 13.5 g/em?. The flight curve is given in Fig. 1. To the 

amount of residual at- 


Qa 


mosphere a total of 
30; 3 g/em? must be added 
80 in order to take into 
El account the packing 
SE Rane material (1 g/em?) and 
x | 1098FPM a layer 5mm thick of 
ey i emulsion over the scan- 
Zl ned area. 
a 30 
220 
a 
& 0 


Fig. 1. — Time-altitude 
curve for Texas flight. 


00 1 
Day time 


3. — Detection of particles. 


A total area of 32.6 cm? of emulsion was scanned at 5 mm from the top 
edge of the stack. In order to improve statistics for large zenith angle tracks 
an additional area of 13.6 cm? (6.8 em? from each side) was scanned at 5 mm 
from the lateral edges of the stack emulsions. 

The selection criteria were the following: 


1) range/plate > 1.5 mm. 
2) number of ò-rays/100 um> 5.5. This value had been found pre- 


viously to correspond in our emulsion to a Z= 8 track (identified by 
its splitting products). 


2) zenith angle 0 <70°. 


(?) M. Artom, V. Bisi and C. GENTILE: Ricerca Scientifica (1958), p. 287. 
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The problem of selecting the tracks to be followed through the stack, out 
of the observed ones, was made difficult by the fact that the range/plate of 
every track was changing near the upper edge of the stack due to emulsion 
distortions. We decided to trace through the stack only those tracks that 
had a range per plate > 3 mm in at least three plates. This criterium ensures 
the measurability of the photometric absorption of the track. 

In order to avoid systematic errors, we added to these the tracks that 
came to an interaction before the end of the distorted region but that would 
have satisfied the measurability criterium if they did not interact. The 
charge of these tracks was estimated by d-ray counting. 

The main disadvantage due to this selection criterium consists in the inde- 
finiteness of the solid angle cone of acceptance and, as a consequence, in the 
impossibility of estimating, in the present work, the H-nuclei absolute ver- 
tical flux. 


4. — Charge identification. 


Results of photometric absorption measurements on 68 H-nuclei tracks, 
entering the stack with zenith angle 0 < 35° are reported on a histogram in 
Fig. 2a. The zenith angle limitation ensures us that the measured tracks are 
all due to H-nuclei of relativistic energy. 

The measurement criteria were the same as described in paper (I). The 
photometric device was adjusted to give the maximum resolution in the region 


to one track 


ee corresponding 


4 
Photometric absorption (%) 


Fig. 2a. — Histogram showing the resolution in photometric measurements. Each 

track is represented by a triangle of constant area with a base equal to twice the 

standard deviation. The dashed areas correspond to the calibration tracks. Only 
tracks incident with a zenith angle 0<3°5 are reported. 
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of charge Z>10. In this way we obtained well resolved peaks in the photo- 


metric absorption histogram (Fig. 24, b). 


Area corresponding 
to one track 


35°<ÿ<70° 


30 35 40 45 50 55 60 65 


Photometric absorption °%lo 
Fig. 2b. — Histogram showing the resolution in photometric measurements for tracks 
incident with a zenith angle 35°< 0< 70°. The broken line areas correspond to events 
which were determined to be not relativistic from the analysis of the interaction stars. 


In order to establish a correct correspondence between photometric ab- 
sorption and charge values, events in which charge could be unambiguously 
determined from the splitting products had to be carefully analyzed. Unfor- 


rea 
I 11 
E 
18+ [ i 
L o-S 
14+ 3 
4 r= 
IS (a) 7° 
10} fo cha 
ae 5 
ts L 
8 
ISS i 
2-19 
Lr DRE 2 TE TY Te [ey Tr 1 
70 90 no BO MUSCLES + 
3 5 7. 9 nN 2 
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(b) a 
Area corresponding 
to one track 
Z=5 Area corresponding 
to one track 


Fig. 3.— a) Total gap-length as a function Fig. 4. — a) d-ray density per 100 um as 
of Z. b) Histogram showing the resolution a function of Z. b) Histogram showing the 
in gap-length measurements. The dashed . resolution in $-ray counting. The dashed 
areas correspond to the calibration tracks. areas correspond to the calibration tracks. 
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TABLE I. — Calibration events. 
Incomi Presi 
ei Target | Splitting Charge measurements 
A | = Nr 
Li | n xp Total gap length 
Be n | 2a Total gap length 
| | 
| 
B n | Li+a Total gap length and 5-rays counting 
| n | 5p Total gap length 
c n 2a+2p Total gap length 
n | 6p » 
_—_—_—$— | 
nucleus | Be+2p | Total gap length 
(N,=4, N,=0) | 
(0) nucleus | B+3p | Total gap length and 3-ray counting 
À (N3= 3, N = 0) | 
Ne p (grey) | B+x+3p | Ne: Phot. absorption; B: total gap length 


and $-ray counting. 


age n 0+2x Mg: 5-ray counting and Phot. absorption ; 
0: total gap length and 5-ray counting. 
nucleus (*) N+a+3p | Mg: 6-ray counting and Phot. absorption ; 
| N: total gap length and 3-ray counting. 
nucleus (*) | Li+3x+3p Mg: Phot. absorption ; Li: total gap length. 


Al nucleus | Na+ 2p Al: Phot. absorption; Na: Phot. ab- 
(N,=18, N,=9) | sorption and total gap length. 

A n | P+a+p A and P: Phot. absorption. | 

Ca p (grey) | P+a+3p Ca and P: Phot. absorption. 


(*) In these events energy determination allowed us to definitely discriminate - splitting 
protons from the other relativistic prongs of the star. 


57 - Il Nuovo Cimento. 
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tunately such events become less frequent for higher charge values of primary 


nuclei. 
All calibration events are reported in Table I. In most of them a heavy 


secondary particle is emitted and the accurate determination of its charge is 
essential in the interpretation of the event. Three methods have been used 
in the charge determination of heavy secondaries: 
1) Total gap length for tracks of charge 3 < Z<11 (Fig. 3). 
2) S-ray counting for tracks of charge 5<Z<12 (Fig. 4). 
3) Photometric absorption measurements for tracks of charge 9< Z. 
The photometric absorption values could then be plotted as a function 
of Z (Fig. 5) and were found to lie on a regular curve which can be represented 
by an analytic expression of the type: 


| Am 
A, = a(1 — exp [— 0(Z — ¢)]) , 990 0/0 


where: a=66,67; b — 0,073; c=1.45. 
The large peak at the upper limit 45+ 


of the spectrum belongs to this Am=a(t-e 9-0) 
curve, corresponding to a charge 40F 145 
value Z = 26. 35 62007 
In Fig. 2b are reported on a his- | =66.67 
togram 33 tracks entering the stack 30; 
with a zenith angle 35°< 0 < 70°. 10 12 14 16 18 20 22 24 26 
Up to a charge value Z = 14 we Fig. 5. — Variation of the photometrie 
calculated that at 41° geomagnetic absorption with the charge. 


latitude all measured tracks are 

relativistic. For charge values Z > 14 a percentage of the tracks incoming 
at the largest angles have a value of B < 1. Since most of the tracks in this 
charge region interact in the stack, an evaluation of the energy could be de- 
rived from the study of the interaction stars and as a consequence an indi- 
cation for the charge could be obtained. 


5. — Interaction mean free path in emulsion. 


The interaction mean free paths in emulsion for all H-nuclei and for nuclei 
of charge Z>15 were determined tracing through the stack the 117 tracks 
recorded. The values obtained are: 

A (emulsion) = (42.9 + 4) g/em? for H-nuclei, 
À (emulsion) = (28.0 + 5) g/em? for nuclei of charge Z>15. 
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This set of values was used in the extrapolation of the charge spectrum 
to the top of the stack. 

The result for H-nuclei is in very good agreement with the one reported 
in paper I. 


6. — H-nuclei absorption in the atmosphere. 


The analysis of the H-nuclei flux dependence from the zenith angle allows 
us to determine the absorption mean free path in air. Such method, of deter- 
mining the absorption mean free path is justified only if the cosmic ray flux 
incident on the top of the atmosphere, averaged over all azimuth angles, does 
not depend from the zenith angle. A recent experiment of DANIELSON and 
FREIER (*) seems to prove that this is the case at geomagnetic latitude 
41° North, for zenith angles 6 < 65°. 

In Fig. 6 are plotted in arbitrary units the flux values as a function of 
the thickness of air traversed, a) for H-nuclei, b) for nuclei of charge Z> 15. 


2 215.3230 g/cm? 


A=117:30 g/cm? 


@ (arbitrary units) 


Atmospheric depth gfem® 
2 


14 16 18 20 2 24 6 28 30 


Vig, 62 — Dependence of the C.R. flux from the thickness of air traversed: curve a): 
for H-nuclei (Z >10); curve b) for nuclei of Z> 15. 


(3) R. E. DANIELSON and P. S. FREIER: Phys. Rev. 109, 151 (1958). 
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Every point was obtained as a result of a weighted mean of the T-stack (pre- 
sent work) and Y-stack (paper I) data. Statistical errors on each point are 
indicated. 

The absorption curves are well represented by exponential laws. The 
absorption mean free path values, corresponding to the best fit exponential, 
are the following: 


À' (air) = (15.3 + 3.0) g/em? for H-nuclei, 
À (air) = (11.7 + 3.0) g/em* for nuclei of charge Z > 15. 


These results are in agreement, within the limits of statistical error, with 
the experimental values of FREIER et al. (4) and are not inconsistent with the 
values calculated following the method proposed by BRADT and PETERS (°), 
if the value for the fragmentation probability P,,, is small (= 0.1). It has 
been pointed out (*) that for atmospheric depths < 15 g/em?, the value of 
the H-nuclei absorption mean free path greatly increases, due to the fact that 
the main effect ot H nuclei interactions at this stage results in a change of 
composition within the group itself, more than in an overall variation of the 
H-nuclei flux. Im this case the H-flux dependence on atmospheric depth 
could not be represented by an exponential law. 

The knowledge of a larger amount of data obtained at smaller atmospheric 
depths will clarify this point. 


7. — Discussion of results. 


The histogram in Fig. 7 gives the charge distribution of H-nuclei, entering 
the stack with a zenity angle 0 < 35°, at flight altitude, as obtained by adding 
to the results of the present work (T-stack) the previous results of paper I 
(Y-stack and B-stack). This could be done since latitude (78) and atmospheric 
depth of the three flights are comparable. The total number of tracks reported 
in the histogram of Fig. 7 is 153. 


The good statistics and charge resolution guarantees the reliability of the 
results. 


(4) P. S. Freier, G. W. ANDERSON, J. E. NAUGLE and E. P. Ney: Phys. Rev., 
84, 322 (1951). 
() H. L. Brapr and B. Perers: Phys. Rev., 77, 54 (1950). 


() M. V. K. Appa Rao, S. Biswas, R. R. DANIEL, K. A. NEELAKANTAN and 
B. Peters: Phys. Rev., 110, 751 (1958). 


(7) C. J. Wappineron: Nuovo Cimento, 6, 748 (1957). 


(5) J. A. Simpson, K. B. FENTON, J. Katzman and D. C. Rose: Phys. Rev., 
102, 1648 (1956). 
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In Table II are reported the ratios of nuclei of charge Z>15 to H-nuclei 
and of Iron to H-nuclei as obtained a) at flight altitude, b) at the top of the 


| 


le) [1 [] (el Le 
10 n 12 13 14 15 16 17 18 9 20 21 22 23 24 25 26 Z 


Fig. 7. — Relative abundance of H-nuclei at flight altitude, as cbtained adding the 
results of experiment I (B and Y flights), and of the present experiment (T flight). 


40 


30 


20 


10 


atmosphere. The extrapolation to the top of the atmosphere was made under 
the following assumptions: 
sorbed according to exponential laws; 
2) the absorption mean free path values are the ones obtained in Sect. 6 
from a zenith angle dependence analysis. For Iron we used, as an 


1) H-nuclei, nuclei with charge 10<Z<15 and Iron nuclei are ab- 


approximated value, Ape (air) = 11.7. 


Taste II. 
N(Z > 15) N(Z= 26) 
| N(Z>10) N(Z>10) 
Flight altitude 0.29 + 0.03 0.08 + 0.02 
| Top. of atmosphere 0.65 + 0.07 0.11 + 0.02 
| C.R. source — 0.16 + 0.03 
Suess and UREY | 0.10 0.05 
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As pointed out in Sect. 6, this method of extrapolation through the upper 
layers of the atmosphere may be inaccurate, if the assumption of an expo- 
nential depth dependence of H-nuclei will be proved to be incorrect. 

An estimate of the ratio Iron to H-nuclei at the C.R. source was obtained 
by assuming an exponential attenuation for both the Iron and H-nuclei com- 
ponents. 

The cross-sections used in extrapolation are the ones quoted by RAJOPADHYE 
and WADDINGTON (*). The composition of interstellar gas was assumed to 
be of 90% H and 10% He, with an average density of 0.1 atoms/cm? (1). 
The mean path traversed by cosmic rays was taken to be R= 1.85-10?5 cm. 
This value was obtained using a set of one-dimensional diffusion equations 
with the following values of the parameters; for fragmentation probabilities in 
interstellar matter a weighted average of the Chicago (11) and Bristol (*) data: 


Pay = 0.29 Po US Pi = 0.07 
Pyy = 0.36 P,,, = 0,34 

Pin, = 0.27 

Prsoretroa = 95 


for flux values at the top of the atmosphere, those obtained in paper I 


D; = 1.67 particles/m? srs, 
Do = 5.52 particles/m? sr 8 , 


Dio = 2.82 particles/m? srs . 


The results of this calculation are reported in Table II. 

A comparison with the SUESS and UREY data (7) (last row of Table II) 
definitely indicates that the composition of the C.R. source differs from the 
average chemical composition of the universe, being much richer in Iron and 
15 < Z < 26 elements. 

From an analysis of the histogram of Fig. 7, giving the charge distribution 
at flight altitude, it appears that nuclei of even charge are more abundant 
than nuclei of odd charge. Interactions occuring during diffusion through 


the interstellar and atmospheric matter must have reduced the strength of 
this effect. 


(*) V. Y. RAJOPADHYE and C. J. Wapplineron: Phil. Mag., 3,9 (1958). 

(9) V. L. GINZBURG: Progress in elementary particles and C.R. Physics (1958), p. 340. 
(1) H. Kosuipa, G. ScHuLTz and M. ScHEIN: Nuovo Cimento, 9, 1 (1958). 

(12) H. E. Surss and M. C. Urey: Rev. Mod. Phys., 28, 53 (1956). 
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RIASSUNTO 


Riunendo i risultati ottenuti misurando le cariche di 117 nuclei pesanti di carica 
Z>10 a quelli ottenuti in un precedente lavoro, si è determinato il valore del rapporto 
N(Z >15)/N(Z > 10) all’altezza di volo e alla sommità dell’atmosfera, e del rapporto 
N(Z=26)/N(Z>10) all'altezza del volo, alla sommità dell'atmosfera e alla sorgente 
della Radiazione Cosmica. I risultati indicano che la composizione chimica della sor- 
gente della R.C. differisce dalla composizione media dell’universo, essendo più ricca 
di nuclei pesanti (Z > 15) ed in particolare di Fe. Si è anche determinato il cammino 
medio di assorbimento nell'aria analizzando la dipendenza del flusso dei primari pesanti 


dall'angolo zenitale. 
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IL NUOVO CIMENTO VOL AX SAN 6 1° Dicembre 1958 


On the Connection between the S-Matrix 
and a Class of Non-Local Interactions - I (*). 


Ku. CHADAN 


Ecole Normale Supérieure - Université de Paris - Paris. 


(ricevuto il 7 Ottobre 1958) 


Summary. — The problem of finding a phenomenological separable 
interaction from the S-matrix has been resolved by M. Gourdin and 
A. Martin who gave exact and explicit formulae relating the interaction to 
the phase-shift and bound state binding energy (1). We generalize their 
method to the case where the interaction is the superposition of a non- 
local separable interaction and a local one, both of them being central (3). 
The local part of the interaction is supposed to be known. We show how 
the separable part can be determined by the S-matrix. Except for a 
few modifications, the method and results are analogous to those of (1:7). 
The main result is that, given the local part of the interaction, the phase- 
shift and the bound states binding energies, it is then possible to determine 
exactly the separable part of the interaction if the phase-shift has a good 
behaviour (4). 


1. — Introduction. 


The mesonic theory of nuclear forces suggests that the interaction between 
two nucleons is local at great distances but becomes non-local and singular 
if the two nucleons come close together. For a phenomenological treatment 
of nuclear forces, it is then natural to superpose a non-local interaction on a 


(*) Supported in part by the United States Air Force through the European Office 
Air Research and Development Command. 

() M. GOURDIN and A. MARTIN: Nuovo Cimento, 6, 757 (1957). 

(*) M. GourpIn and A. MARTIN: Nuovo Cimento, 8, 699 (1958). 

(*) The case of tensor forces will be considered in a next paper. See also O.R.A.S 5 
245, 1597 (1957). 

(4) A preliminary report of this work has already been published in CR AS, 
246, 1513 (1958). 
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local one. We can suppose that the local part of the interaction is determined. 
so as to agree with experimental data at low energies. 

To determine a non-local phenomenological interaction from the S-matrix, 
it is always necessary to make some assumption about it. An arbitrary non- 
local interaction is too general to be uniquely determined by the S-matrix. 
There are a great variety of non-local interactions which are equivalent as. 
far as their agreement with experimental data is concerned. 

In absence of any detailed theoretical information about the nuclear forces. 
at small distances (2), it is plausible to assume that, near the origin, the inter- 
action is separable. This assumption has, as we shall see later, the advantage 
of simplicity. Then, there is some interest in connecting the separable part 
of the interaction to the S-matrix, the local part being known from low 
energy experimental data. 

We take an interaction of the form (?): 


Dr, r°) = V(r) der) + D emntr)utr') Pal cos 6) ; 
(1.1) 0 


r= |r], en À A =), ét EL. 


We have then restricted ourselves to purely central forces. 
The Schrôdinger equation therefore reads: 


2 


h 
2M, 


(1.1) AY(r) + EV (r) =| U(r, r')¥(r')dr'= 


= VPI) + > er] u (r) u(r’) Pi (cos 6)¥(r') dr’ . 


M,=reduced mass = M/2. In what follows, we put #—2M,=—1. 

It is well known that the /-th term of the expansion in the right-hand side 
of (1.2) acts only in the state of angular momentum I. ‘e, is therefore the sign 
of the separable interaction acting in this state. 

In the present paper, we consider only the case of zero angular-momentum 
(S-state). The extension of the method to higher angular momenta is straight- 
forward... 

The similar problem in which one tries to find a local potential from its 
S-matrix has extensively been treated, by several authors: LEVINSON, BARG- 
MANN, JOST and KoHN, GEL'FAND and LEVITAN, NEWTON and FULTON (°°). 
Because of the presence of a local interaction in our problem, we make exten- 
give use of the results or 9). 


(5) R. Jost and W. Konn: Dans. Mat. Fys. Med., 27, 9 (1953). 
(6) R. NEWTON and T. Futon: Phys. Rev. 107, 1103 (1957). 
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2. — S-state Schrodinger equation. 


21. Existence of the solution. — The Schrodinger equation in the case of 
zero angular momentum is: 


(cel 


(2.1) y"+ By = EU Ur) Ter) ple) dr! + Vir) pr), 


Un) = Vis ruo(r) , BS ys 


where y is the reduced radial wave-function. 

In order that this equation have a solution V(r) and U(r) must satisfy 
some conditions near the origin and near infinity. As we shall see later, it is 
essential that the «local» Schrédinger equation: 


(2.2) gy" + Ep = Vo 


offers good solutions. V(r) is then assumed to satisfy (*): 


œ 


(2.3) fr Vir) |dr<.00, 


(ce) 


(2.4) fr [F{(r) de es. 


0 


On the other hand, we suppose that U(r) satisfies the following conditions (7): 


(2.5) U(r) me a< à, 
(2.6) U(r) ~r? pae 


2°2. The «local » Schròdinger equation. — Here, we recall some properties (5) 
of solutions of (2.2). Consider first the case of E > 0 (scattering). (2.2) has 
many interesting solutions, each of them having either a given asymptotic 
form or a given behaviour near the origin. The most interesting are: 

(a) Î(+k,r): limexp[+dkr]f(+k,r)=1, OR A 
(b) P(E, r)= ok, 7): g(E,0)=0,  p'(E,0) =1; 
p(#, re a sin (kr + 0,), 


0, being the «local» phase-shift. 


(7) These conditions are obtained by considerations quite analogous to those of 
the begining of (!). Recall that U(r) = ru(r). 
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In what follows, we write indistinctly q(E,r) or g(k?, r) when H=k?>0, 
The solution g(k?, 7) can be obviously written as: 


1 
(2.7) p(k, r) = a IC kr) — f(— k) f(k, r)], 
f(k) =f(k, 0),  f—h)=f(). 


It is easily seen from (2.7) that: 


(8) pth, e) > SPL ed) sin (hr + 09); 
2.9) exp [2id] = MMM - 


If there are bound states, their binding energies are given by the zero of f(k) 
in the lower half of the complex k-plane. Under the conditions (2.3;4), 
f(k) has a finite number of simple zeros along the negative imaginary axis. 
Call these zeros —iy:, 1=1, 2,...,. There are therefore n bound states with 


the binding energies —y;=— £,. The corresponding eigenfunctions are gi- 
ven by: 
Rs È 
(2.10) pir) = 9(— Br) = afl) fiv”), 
U 


“a 


g(r) is not normalized to unity. 
The function f(k) enables us to define the spectral function o(£) by: 


0(0)= 0, 
do 
(2.11) de A k_ RO 
dE x |f(k)f a 
where 


(E fie ar| - 


The spectral function has the interesting property that for every square 
integrable function F(r) we get the completeness relation: 


co 


(2.12) exon — |do(£) fr (r) g(2, 1) ar : 
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This means that: 
+0 
(2.13) | o(B, r)g(E, r')do(E) = d(r—1"). 


— oO 


Moreover, the wave-functions g(k?, r) and y,(r) obey the usual orthogonality 
relations: 


- 2 __Ô(E—E) _ a]|fòk_k) 
[ove ,)g(k'?, r) dr = do/dE = Dh? , 
(2.14) | (pr) dr = cr * Ow , 
| (r)g(E, r) dr = 0 E=k?>0. 


These formulae permit us to define a kind of integral transform which is very 
similar to a Fourier sine transform. For any square integrable function F(r) 
we define 


F(k) = front, r) dr, 
(2.15) 


Tee i F(t) pilr) dr . 


The inversion formulae read: 


+0 


2.16) Fl) = fdo(M) FB) 9B, 1) = 5 arr + | dg(H) FH) gH) 


— © 


2°3. Solution of (2.1) with use of integral transforms. — We write: 


(2.17) vir) =| pl te" ae, 
(2.18) Ute) = [UE MAP). 
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The equation (2.1) then reduces to: 


+o +0 


(2.19) fauve, t= ex [tete U(E')g(E',r), 
(2.20) N =| U(t) p(t) dt LME 


Scattering, H=k. The wave-function is represented by p(k, 7) and. 
its transforms by y(k, p) and p(k). (2.19) gives us: 


9 e On U: 
(2.21) pilk)=e ke LE, N(k), 


(2.22) (k® — p?) plk, p) = eN(k) U(p) . 

The solution p(k, 7) which is of interest must have the following properties: 
lk, 0) = 0, 
p(k, r) ~ Asin (kr +0), 


to 


ÿ being the total phase-shift. Moreover, p(k, r) must coincide with œ(k?,r) 
when U(r)=0. The solution of (2.22) satisfying all above conditions is: 


Therefore: 
Pi ei pil”) on "| Ue") Ss 
(2.24) p(k, r) = > ec,U, N ea + £N (| ae (p?) P ger i + p(k?, 7) . 


0 


By replacing (2.23) and (2.21) in (2.20) we find: 


U(k 
(2.25) N(k) = = = 
SU tes st 
ef pe 000) 2°p4E 


al 


0 


When r goes to infinity, the first terms in (2.24) go to zero and therefore the 
behaviour of p(k, 7) is given by the remaining terms. It is easily seen that 
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one has (8): 


(2.26)  y(k,r) ~ exp T= 50] f(k) [sin (kr + 6,) + tg à, cos (kr + 6,)], 


T—o k 
F(k 
(2.27) È tg du = x a 
F(p) ec; Ui 
1 Pela ed Le LB 
2ek? U?(k) 
BE = 
1) = Hk) D 


The right-hand side of (2.26) is of the form A sin (kr+6), as expected, with 
6 given by: 
(2.28) 6=6,+ 6, 
6,(k) is the additional phase-shift due to eU(r)U(r'). The superscript v indi- 
cates that a local potential V(r) is present. 

If V(r)=0, then g(k?, r) = (sin kr)/k and f(k)=1. (2.27) reduces to the 
classical form : 


(2.29) 233 tg du = di 
k  F(p) 
1 + ap "a dp 


2'4. Bound states. — If there exists a bound state with binding energy +X 
(X<0), (2.19) leads to: 


Lee] 


Derw:(E. + X)qi(r) fam voy = Appia 


(ce) 


ENI D ¢, Up. (r) + facts U(P)g(p?, DI . 


0 


(8) To do this calculation, we substitute p(p?, Tr) by (2.7), 


1 b 1 1 1 
p°—k? di e 


and then we get (2.25) by calculating the integral along an appropriate path in the 
complex p-plane. 
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Suppose first that all U, 40 and U(p) #0. We deduce that: 


NU 

2.30' PRAIA 

( ) Yi E, X ? 
—eNU(p) 

2.30" ites: See 

( ) y(p) or 


By replacing these expressions into (2.20), we get the eigenvalue equation : 


eo 


(2.31) press nes P| agin” Send 


p? 


More generally, we consider the equation (2.31) without any restriction on 
the sign of X. If X=? >0, @(X) is nothing but the denominator of (2.27) 
multiplied by —e. In 
order to determine the 
number of negative roots 
of (2.31), we make use of 
a graphical method. 


a) e=+1. A(X) is 


represented in Fig. 1. One ye 
sees immediately that if 
ce, Ur 
00) = 2, — n=2, €=*1, D(O)<O, V:n:2 
eni) 
— facto Re ni Fig. 1. 
0 : 


is positive, the eigenvalue equation has n — 1 negative roots. Otherwise, there 
are n eigenvalues. One sees also that the absolute value of all these eigen- 
values are smaller than the corresponding Hs. This is to be expected, the 
separable interaction being repulsive. 


b) e=—1. This case is very similar to the case a). The conclusions is: 
6, Ui; U?(p) 
if Hie ooo lit 10) 
L 


0 


there are n eigenvalues ; if @(0) < 0, there are n-+1 eigenvalues. 
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Hereafter, we designate the negative roots of (2.31) by — IE ae) 
and the corresponding eigenfunctions by w,(7). y,, Means Jeske) esto) der. 


Consider now the case in which one of the U,’s, say U,, is zero. a 
P(X) is continuous at X=— €; and one of the {— /"7} desappears (see Fig. 2). 
As it can be directly 
seen from (2.1), the 
eigenvalue which desap- 
pears is replaced by 
— E;, the eigenfunction 
being @;(r). This means 
that if U(r) is orthogonal 
to ,(r), the superpo- 
sition of eU(r)U(r') to 
V(r) does not change 
the wave-function ¢;(7) 
n=2,€=-1,U,=0, @(0)<0 and its binding energy. 


yin+1=3 The remaining eigenva- 

The state -E, is degenerate lues are given as above, 

Fig. 2. i.e. by (2.31). From these 
later, one, say —J%, 


‘can be equal to — E,. In this case, we have a degeneracy at — E;. If it is 
so, the second wave-function, y,(r), can be chosen orthogonal to ;(r). 
This is simply done by taking Y= 0. 

The case where many U, are zero is very similar. The order of degene- 
racy at each — E, can never be greater than two. 

If all U,=0, there are at least » bound states, namely all {Æ,, gi(r)}. 
Therefore if v=n —1, we are sure that at least one of Us is not zero. 

Finally suppose that U(k)=0, i.e. U(r) = Sapir). (2.31) reduces to: 


c, U# 
E +X 


(2.32) S Me 


Here, however, one can never get more that n negative roots. This is because 
(0) > — e. The bound states wave-functions are linear combinations of Pi(r). 
Besides, the scattering wave function p(k, r) is identical to (k?,r). This is 
clear from (2.1). The only change concerns the bound states. Remark that 
(2.27) gives 6°= 0. 
The equation (2.32) can have sometimes one positive root. This occurs 

when ¢=+1 and D(0)>0 This positive root corresponds, as we shall see 
later, to a positive energy degenerate state. 
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2°5. Positive energy degeneracy. — Consider again the equation (2.22). We 
have resolved it by imposing to y(k,7) the asymptotic form A sin (kr+6) 
However, (2.22) can admit solutions with vanishing asymptotic form. on 
pose that x(k, 7) is such a solution. y(k, 7) must vanish identically if U(r) = 
because the local potentials considered here admit no such solutions. (2.22 4 
gives immediately: 


U; 
2.36 NS ey 
(2.33) x(k) =e Bik Ny 
‘ se titi) 
(2.34) vk, p) = ue ii 


By replacing these values into (2.20), we are left with: 


Le È U? 
(2.35) D(k ee P| dein) a Peay. 


We see that positive energies at which degenerate solutions can exist are given 
by positive roots of the eigenvalue equation (2.31). However, this condition 
is not sufficient to insure the existence of y(k, 7). We must verify that y(k, 1) 
has really a vanishing asymptotic form. (k,r) is given by (dropping the 
normalization constant N): 


U, U(p 
“(ky =e> 3 2, PH) — Paden" ee p(p*, 7) . 


At infinity, only the integral contributes. This integral is evaluated as follows: 


—p,r).—}{—-p)}p,r)}= 


po L la 
P| deo") 37 HP = P | dow") 


0 


p? — A 


1 (pap U®) a ie io 
ere. r~ — ; exp [ipr] = 
alia pr") ri] jp) pe — 
Li U(k 
Fie 
Therefore, 7(k, 7) has a vanishing asymptotic form if and only if U(k) —0. 


We conclude that: 
Theorem: The necessary and sufficient conditions to have a degenerate 


state at the positive energy E,— ke are: 


py [exp [ikr] + exp [— ikr] exp [— 270,]] . 


(2.39) (ko) = 9, 
(2.40) U(ky) = 0. 
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The breakdown of the famous Levinson theorem is intimately related to 


the existence of such states (91°). 
In the special case considered in the end of Sect. 24, if (2.32) admits a 
one obtains such a degenerate state. For the corresponding 


positive root has 
and therefore vanishes at in- 


wave function is a linear combination of gi(r) 


finity. 


9°6. Some properties of Oy. — 6° is given by: 


co 


= 2k 7 ° F(p Seau: 
van aan + Pee? Lp © ln 
i iti 


F(k) has a constant sign, the same as €. Equation (2.41) shows that tg ò? 
goes through zero at k=k, if and only if: 


O(k?) = 0. 


These conditions, the same as (2.39), (2.40), are sufficient to insure that tg 0, 
woes effectively through zero at k=%k,. This is shown as follows: : 
Due to (2.3) and (2.4), q(k?, r) is a well behaved function, i.e. it is a con- 
tinuous function of r and has a continuous derivative 09/0k. It is also an 
entire function of H=k?. Then, (2.5) and (2.6) insure that U(k) is also a 
well behaved function. Its derivative with respect to k, dU/dk, exists and 
is continuous everywhere except, may be, at the origin. Near each zero, k 
of U(k), we can therefore write: ee: 


(2,42) U(k) = (k — ko) U(ko) ; | T (ke) |< 00. 


k—>Ky 


Then, F(k) has at least a double zero at k= ky: 


(2.43) F(k) = e(k — ko)? F(ko), 0 < Fk) < co- 


kaka 


() N. Levinson: Dans. Mat. Fys. Med., 25, 9 (1949). 
(10) A. MARTIN: Nuovo Cimento, 7, 607 (1958). 
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This last equation shows that 


°F) | Ff F(p) 
Gr pr? | eh (| (p? — RP ar) 


exist and the latter expression is the derivative of the former with respect 
to ki. It follows that: 


ap Ci U; i eF(p) | 
2.44 — = — 2k P aa ; d —i 0. 
moe di, =. CEA ca i 
0 


Therefore, in the neighbourhood of each k,, one can write: 


— 2k — ko} FK 


2.45 
“Sead 2 (= “E n a ko) 


aa A>O. 


We infer that tg 6” goes effectively through zero at k= k,. Incidentally, one 
sees also that: 


(2.46) = 0e 


One has therefore the two following theorems: 


Theorem a): A necessary and sufficient condition for tg 0, to go through 
zero at k= %, is that there exists a degenerate state at E, = ks. 


Theorem bd): 6° goes once and only once through each na. 


Most often, the zeros of the numerator and of the denominator of (2.41) 
do not coincide and the sign changes of tg 0, correspond to zeros of the 
denominator. In this case, dò, remains between nz and (n-+1)z. 


27. The generalized Levinson theorem. — As we shall see later, the value 
of 6°(0) — 0,(00) is of great importance when one attempts to solve the integral 
equation connecting U(k) to 62(k). It can be deduced from the Levinson 
theorem (°) generalized by A. MARTIN(?°) for non-local interactions. 

According to this theorem (*): 


6(0) — 6(co) = (v+0)m 
6,(0) — d,(00) = na. 


(*) We disregard the possibility of a resonance at zero energy. 
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tive. 
g=n—l1 and o= 0. This is in contrast to 


6,,(0) — 6,(00) is always non-negative. 


KH. CHADAN 


o is the number of degenerate states at 
positive energies; v and ” respectively 
are the number of bound states of the 
total and local interactions. By sub- 
tracting these two equations, we find: 


Corta vicn sake 520) — dt(co) = @ — 17 + on. 
Fig. 3. 

Generally, 62(0) — d,(00) is non-nega- 

The only permissible negative value is —7; and occurs when #9, 

the purely separable case in which 


nt0- (8541, Ysa 


Fig. 5. 


Fig. 3, 4 and 5 show some typical cases of variations of d?(k) against the 


momentum È. 


3. — Derivation of the integral equation. 


We have to determine U(r), or its integral transform {U(k), U,}, from, 
V(r), Ô(k) and the discrete spectrum. This is done as follows: 
From V(r), we calculate the phase-shift, ô,(k). Then, we determine 6°(k) by 


(3.1) Ok) = Ô(k) — d (4) 
and therefore: 
(3.2) g(k) = — = te x(k). 


dk) must have all properties previously established. 
e, the sign of the separable interaction, is given by the sign of g(k) near 


infinity : 
eg(k) — ASO, LUE 
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g(k) and e being determined, we have the following integral equation connecting 
{U(k), Us} to g(k): 


°F (p) ap ec UE 
(3.3) re L P| HI Ti. 
(k) = g(k) [1 + | pi De nl 
De a 
FD == - U2(k 
BE 0 pe 0) 


This integral equation by itself is not sufficient to determine F(k) and the 
constant U,. Binding energies are also needed. This fact is in agreement with 
general theorems about the connection between the S-matrix and the inter- 
action. There are therefore two steps to obtain U(k) and U,, 1=1,2,..., n. 

We begin by supposing that the constants {U,} are known. In other 
words, we suppose that the inhomogeneous part of (3.3) is given. Some of 
the {U,} may be zero. This can be seen by simple examination of the discreet 
spectrum (see Sect. 2). Solving the equation (3.3) for F(k), we find a solution 
depending on n constants U,. Putting this solution F(k; U,) into the eigenvaule 
equation (2.31), we obtain a set of equations which determine {U,}. As we shall 
see, our method cannot be applied in the particular case v —n—1 and o=0. 

Before solving equation (3.3), we establish some theorems connecting the 
sign of g(k) near the origin and near infinity with the number of bound states 
in order that the integral equation can have admissible solutions, è.e. solutions 
with constant sign. These theorems are based on the discussion of (2.4) and 
the equation 

eF(k) 


g(k) = — due) * 


Theorem 1: If g(0)g(oc)>0, the number of bound states must be 
equal to n. 

Proof: There are two cases, according to the sign of g(00). 

If g(co) >0, we have e=1 and D(0)< 0. This shows that v= n. 

If g(co) <0, e is equal to — 1 and then (0) > 0. This requires also y=”. 

Theorem 2. If g(c) > 0 and g(0)< 0, there is no admissible solution 
unless y= n —1. 

Proof: g(co)>0 implies that «=1, and therefore » cannot be greater 
than n. Besides, g(0)< 0 imposes to ®(0) to be positive. This shows that 
ven — 1. 

Theorem 3: If g(0)>0 and g(co) <9, one must have v=n+1. 

Proof: e is equal to — 1 and D(0)< 0. The number of bound states 
must therefore be equal to n+1, according to previous discussions. 


Theorems 1, 2 and 3 give necessary conditions for the existence of solutions. 
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4. — Solution of the integral equation. 


As in (2), we make the following changes of variable and of function: 


e=h+1, 
Pe) ,,. 9@) 
Equation (3.3) then becomes: 
ec, Ui 1 (hye) 
= -+- : - dy, 
cae De Ii à Pere i 
1 


mg (x) g(x) 


(42) hiv) = ME Ê in Hes sin 62 (x) exp [— id (æ)] . 


The new integral equation is of the same type as equation (2.1) of OMNES ("!), 


with òf=—ny and 


We can therefore apply Muskhelishvili’s method generalized by OMNES. 


Inhomogeneous equation. — We define a function H(Z) of the complex va- 


riable Z by (1°): 


1 [ht 
(4.3) Ha) | ee dy. 
(4.4) h(y) = exp [in(y)] sin 7(y) . 


H(Z) is holomorphic in the complex plane cut from 1 to oo along the real 
axis. Obviously, the solution of (4.1) is: 


(4.5) p(x) = Mx) + 2iH(x*) . 
Now, H(Z) is put into the form: 
(4.6) H(Z) = QZ)D(Z), 


(1) R. OmnÈs: Nuovo Cimento, 8, 316 (1958), equation (2.1). 

(*) To avoid confusion, we use the notations n and À (x) instead of 6 and f(x) used 
by Omnès. 

(12) We suppose that 6; is chosen so that dul co) = 0. 

(18) Strictly speaking, there are some indications in favour of separable interac- 
tions. See, for example, Lévy’s repulsive core: M. M. Livy: Phys. Rev., 88, 725 (1952). 
Also the footnote (6) of Y. YAMAGuCHrS paper: Phys. Rev., 95, 1628 (1954). 


3868 


ON THE CONNECTION BETWEEN THE S-MATRIX ETC. 907 


where 

(4.7) QZ) = exp[u(Z)], 

(4.8) uz) = 2 UO. 
a) y—Z 


1 


@(Z) must be holomorphic in the complex cut plane. It must also go to zero 
in the limit |Z|— oo. From (4.1), one deduces that ®(Z) must satisfy the 
following functional equation: 


(4.9) D(xt) — D(a-)= — = A(x) [exp [— u(a*)]— exp [—u(æ hi; 


u(x+) = tim u(æ + te) = (x) + in(x), 


Si pf 


TT y — 


1 


A particular solution of this equation is: 


À He ec U? È 
(4.10) D(Z) =} er [—(Z)] 402) —1 +5 poe a PL il], 
tk _i ° nly) x 
sat 8) =2 [ay 
1 
We therefore deduce, from (4.5), (4.6) and (4.10) that: 
#004 ‘ 
(4.12) p(x) = exp [u(a*)]]1 oe Pat exp [— vi]. 


One must check that (4.12) is a solution of the integral equation (4.1). 
This is done in the same way as in (?) by assuming that n(0) — (co) < 0 
This means that the solution (4.12) is not valid if v=n—1 and o=0. In 
fact g(x) is no longer finite at æ—1 in this last case. ; 

Homogeneous equation. — The homogeneous equation corresponding to (4.1) is: 


roue dy 


(4.13) Pole) = lag Sweets 
1 


The same method as used for the inhomogeneous equation leads to: 


exp [u(at)] > - Le = if--4{0 — n(oo)) = — Na<0, 
(4.14) po(r)= oe geet) 


0 otherwise ; 
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b,, are arbitrary constants. Therefore, the general solution of (4.1) is: 
a b, ec; Ui ses 

(4.15) (x) = explu(a*)] È ro > TI 2 yay: pes exp [ 1) 


Taking into account the changes of variable and of function, we deduce the 


general solution OL (33): 


SOIT 
y 


7 1 2 
(4.16) F(k) = 22 sin 6% (k) exp [ax(k?)] f + > or «adi exp [— ul ; 
This solution F(k) depends on N+n parameters, namely {U,} and {bm} 
It must keep a constant sign. Therefore, the changes of sign of sin 0, must 
coincide with the changes of sign of the expression between brackets in (4.16). 
This requirement gives us a set of relations between the parameters. These 
relations, added to those deduced from the eigenvalue equation, permit us to 


find the value of the parameters. 

The particular case v=n—1 and o=0 must 
been able to find the solution of the integral equation in this case. 

In all other cases (v > or o > 0), it can be easily seen that there are at 


be excepted. We have not 


least as many equations as parameters. 
The question of the uniqueness of {b,,} and {U,} has not been clarified yet. 


Note added in proof. 


We shall give in a forthcoming paper the solution of the integral equation (4.1) 
in the particular case v = n—1, 6 = 0, as well as the conditions of existence and 


uniqueness for the constants {Uj}. 


RIASSUNTO (*) 


Il problema di trovare dalla matrice S una interazione fenomenologica separabile 
è stato risolto da M. Gourdin e A. Martin che hanno dato formule esatte ed esplicite 
che mettono in relazione l'interazione con lo spostamento di fase e l’energia di legame 
dello stato legato (1:2). Generalizziamo il loro metodo al caso in cui l'interazione è la 
sovrapposizione di un’interazione non locale separabile e una locale, ambedue centrali (3). 
Mostriamo come si possa determinare la parte separabile per mezzo della matrice S. 
Eccettuate poche modificazioni, il metodo e i risultati sono analoghi a quelli di (*:°). 
Il risultato principale è che, data la parte locale dell'interazione, lo spostamento di 
fase e le energie di legame degli stati legati, è possibile determinare esattamente la 
parte separabile dell’interazione se lo spostamento di fase ha un andamento regolare. 


(*) Traduzione a cura della Redazione. 
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(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Polarization of 5 -Particles in *Li-Decay. 


L. KeszrtHELYI and J. ZIMANYI 


Central Research Institute for Physics, Department of Atomic Physics - Budapes 


(ricevuto il 10 Ottobre 1958) 


In their excellent experiment GoLp- 
HABER et al. have measured the helicity 
of the neutrino produced in K-capture (1). 
The helicity of the antineutrino produced 
in 8 -decay, however, can be determined 
only if two joint experiments (£--v an- 
gular correlation and @ -polarization) are 
carried out together. Therefore it seems 
worth while to perform these two expe- 
riments on the same [ -decay. 

MAyER-LEIBNITZ et al. (?) have meas- 
ured the angular correlation between 87 
and v-particles produced in *Li B -decay. 
Since the polarization of { -particles in 
SLi decay has not been measured yet, 
we have carried out this polarization 
experiment. 

The half life of SLi is 0.858, the 
maximum f-energy Ena = 12 MeV (3). 
In order to determine the 8 -polarization 
it has been found convenient, because 
of the high B--energy, to measure the 
circular polarization of the external 
bremsstrahlung of ( -rays (*). 


(1) M. GOLDHABER, L. GRODZINS and A. W. 
SUNYAR: Phys. Rev., 109, 1015 (1958). 

(2) K. H. LAUTERJUNG, B. SCHIMMER and 
H. MAIER-LEIBNITZ: Zeits. f. Phys., 150, 657 
(1958). 

(3) F. AJZENBERG and T. LAURITSEN: Rev. 
Mod. Phys., 27, 77 (1955). 

(4) M. GOLDHABER, L. GRODZINS and A. W. 
SUNYAR: Phys. Rev., 106, 826 (1957). 
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Radioactive 8Li nuclei were produced 
by bombarding a natural Li target o 
about 1mm thickness making use o : 
the reaction ?Li(d, pli. The ~ 650 keV 
deuteron beam was produced by mean s 
of a Cockcroft-Walton generator. The 
8 -rays of SLi initiated bremsstrahlung 
in a Pb absorber. The circular polari- 
zation of bremsstrahlung was measured 
in the usual way (4). 

In order to eliminate the direct 
+ radiation produced by the nuclear 
reactions, the d beam was periodically 
pulsed. An activation of 28 was fol- 
lowed by a measurement of 1s. The 
direction of the magnetic field was al- 
tered per 15s. 

In the Li+d reaction also fast neu- 
trons are produced in a great number 
producing activities rather disturbing the 
measurement. We were particularly 
annoyed by the radioactivities of ?°F 
nuclei produced in Nal(Tl) crystal by 
(n, «) reaction and that of 1$N produced 
in the cooling water and in other ma- 
terial (e.g. the glass envelope of the 
multiplier) containing oxygen, by (n, p) 
reaction. Therefore, a CsI(T1) crystal was 
used instead of Nal as y-ray counter 
and the target was cooled by paraffine 
oil. 

An integral discriminator connected. 
to the detector was applied for meas- 
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uring y-rays of higher energy than 
~ 5 MeV. 

The effect of the magnetic field on 
detectors was studied by means of y-rays 
of Co, Na and nuclear reactions 
7Li(p, y)*Be. It has been found that the 
number of pulses remained constant 
independently of the direction of the 
magnetic field within the error of meas- 
urement amounting to +0.15%. 

Let us denote by + the number of 
pulses in the case when the magnetic 
field in the measuring magnet is oriented 
towards the target (the spins of the 
electrons in the iron are oriented parallel 
to the momentum of the incident y-ra- 
diation) and n— will be the number of 
pulses in the case of an inverse magnetic 
field. In our experiment the value of 


n_— N+ 


n+ ne 


A 


was found to be 


Kite He (3.4 2 08)-10-. 

Since the disturbing activity of !6N 
could not be completely eliminated, this 
value has to be corrected. By analysing 
the decay curve it has been found that 
of the total number of pulses a fraction 
of 72% was due to the bremsstrahlung 
of the @ -particles produced in *Li-decay. 
The correction being taken into account 


ce (4.86 Eli) 10-48. 


For completely left circularly polarized 
y-rays we would obtain with the satura- 
tion magnetization of the analysing 
magnet of 7 cm effective length, applied 


L. KESZTHELYI and J. ZIMANYI 


in our experiment 


, = + 5610 


%100% 


Thus it follows that the average polari- 
zation of the bremsstrahlung produced 
by  -particles in SLi decay, is 


P,= (86 4 20)%> 


for the energy range investigated in our 
experiment. The y-rays are left handed 
circularly polarized (their spin is anti- 
parallel to their momenta). It follows 
from this fact that the 8--particles of 
sLi are polarized in backward sense. It 
is rather complicated to determine 
accurately the degree of 8 -polarisation 
from the value of Py. (The shape of 
the 8 -spectrum and that of the brems- 
strahlung, as well as the energy depend- 
ence of the polarization of bremsstrahlung 
y-rays and the energy dependence of 
counter efficiency and so on, have. to 
be taken into account). 

In our measurement the sign of 
polarization is unambiguously obtained 
and the degree of polarization of electrons 
is found to be in agreement with the 
relation P, = v/c. 

By comparing our experiment with 
that of MAYER-LEIBNITZ et al. it follows 
that the spin of the antineutrino pro- 
duced in 8 -decay is parallel to its 
momentum in accordance with the result 
obtained by GOLDHABER et al. @). 

* OK * 

We have to thank our colleagues 
Dr. J. Erò and T. Nagy for their help 
in planning the experiment and L. SZABÔ 
for building the electronic equipment 

and G. Scumipt for his help in carrying 
out the measurements. 
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The Selection Rules on Hyperon Decays. 


K. NAKAGAWA and H. UMEZAWA 


Department of Physics, University of Tokyo - Tokyo 


(ricevuto il 21 Ottobre 1958) 


Recently we have published a short note (1) with the same title. There, use has 
been made of the observed branching ratio 


(1) Rs = 


the ratio of the life times 


{2) R's = 


and the ratio of the anisotropy coefficients 
(3) clou = Ale 


The symbol x denotes the anisotropy coefficients for the X-decay processes; a for 
X++p+-7° and «+ for X++n+x*. It was concluded in the previous note that the 
rule |AZ|= } does not hold as far as the one-to-one law and the invariance under 
the time reversal (T-invariance) are assumed. The one-to-one law says that the 
coupling constants for the parity conserving and non-conserving terms are of a 
same magnitude. ; 

We now come back to the same problem, because recent experiments have not 
always supported (3); the situation for the anisotropy coefficient now is quite 
obscure (*). 

Assuming the rule |AZ|= 3 and 7-invariance, the X-decay interaction is of 
the form 
— Lt, 


1—7T 
> (A+ Bys)VuWzidu®zi + (Fe 


(4) H= yx 


Je + B'ys)vufbsti Ou Pt he. 


Notations here are as follows, 


Yrs Pr3 


(1) K. NAKAGAWA and H. UMEZAWA: Nuovo Cimento, 8, 945 (1958). 
(*) The authors express their sincere thanks to Prof. F. S. Crawrorp for information on 


the experimental results under subject. 
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1 ; | CD pe 
| Pot = 7/2 (Py, È PTE Prt = V2 (Gri + Pra) » 
| pro = Yx3 Pr = Pas» 
0 0 0 0 0 a 0 ‘—i 0. 
fel 0 0 —i f=) 2 0 0 4 =\|i 0 0:35 
0 v 0 —1 0 0 0 0) 0 
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Here, the constants A, B, A' and B' are real. The one-to-one law requires that 


f AB, 
6) | AB", 


This together with experimental results (1) and (2) leads us to the relations 


AJB = — A'/B’, 


wi | 4224”. 


— 


It should be noted that A/B = 4'/B' was obtained in our previous letter where (3) 


was assumed. The relations (5) and (6) lead us to the interaction 


1 — 


OT uv. Ta A MTS È 
(7) = JPN (ht ers) %u PEs Ou Pat ES YN Ti (1—eys) Vp ati Ou Pa > 


2 


where g is a real coupling constant and e= 41 
calculated to be 


The anisotropy coefficients are 

Oy (+> p+r) 

a,(it>n+nt)=0, 
( 


(Sone) 0. 


6, 


(8) 


The ratios Ry, and Ry, are shown in the table. It should be noted that the one- 
to-one law depends on the expression for the interaction; it does on the definition 
of the coupling constants. Indeed, the one-to-one law would be lost if it would be 
formulated in terms of the coupling constants for the y,yx*p ete. 


This is a situation different from that in our previous note where the relation (3) 
was assumed. 


Tape I. — Branching ratios, the ratios of the life times and the anisotropy coefficients. 


A'JA | BJA =e | B'JA"= — € Rs Rs % CAS o 
I 1 — ] 0.88 0.44 — 1 0 0 
1 — 1 1 0.88 0.44 1 0 0 
— 1 Il — 1 1.28 0.75 — 1 0 0 
— I — 1 1 0.75 1 0 0 
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Inelastic Scattering of Electrons from ‘Li and #N. 


P. CAZZOLA 


Istituto di Fisica dell’ Università - Padova 


CC. HOGLLA 


Istituto di Fisica dell Università - Parma 


(ricevuto il 15 Novembre 1958) 


The inelastic scattering of electrons from 1?C has been calculated by MORPURGO (!) 
making use of the shell model in its most simple form, both in LS and JJ coupling. 
The agreement of these calculations with the experimental results is rather good: 
the very strong increase of the ratio inelastic to elastic scattering with the energy 
is given correctly and also the absolute value of the inelastic cross section differs 
at most by a factor of 3 from the experimental value. 

It is then reasonable to ask whether an analogous situation is true also in other 
nuclei. We shall give here the results of calculations similar to the above mentioned 
ones for the case of *Li and !4N; unfortunately no experimental data are as yet 
available on these nuclei; but we may reasonably hope that the experiments will 
be performed soon, and these calculations may be also useful in getting some idea 
of what one has to expect. 

The reasons why we have chosen the above nuclei are: 1) the lower levels are 
sufficiently distant from each other so that the elastically and inelastically scat- 
tered electrons may be energetically resolved; 2) the shell model structure of the 
nuclei in question is rather simple; 3) the spin and parities of the levels involved 
are known. 

Concerning the question of the resolution (point (1) above) we would like to 
mention that, in principle, also the cross section for excitation of rather close levels 
may be measured provided one takes the coincidences between the scattered elec- 
trons and the y-rays from the de-excitation of the level in question. Such a technique 
is certainly much more complicated, but, once established, might give the pos- 
sibility of getting by inelastic electron scattering very important information on 
a large variety of nuclei. 


(®) G. MoRPURGO: Nuovo Cimento, 3, 430 (1956). 


~~ 
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Coming back to the present case we have, for the ‘Li, considered the scattering 
from the ground state (J=1*, T=0) to the first excited state (7=3*, T=0) at 
2.18 MeV. We have used in this case only the LS coupling — in which the two 
states in question are 39 and 3D respectively — because it is well known that 
the JJ coupling does not give correctly the order of the levels in SLi. 

For the 14N we have considered the scattering between the ground state and the 
state at 3.95 MeV which both have J=1* and T—0. In this case the calculations 
have been performed both in LS and JJ coupling. The calculations are very similar 
to the ones performed by MORPURGO for the 12C; we have neglected recoil effects 
and the corrections arising from the form factor of the nucleons. Unless explicitly 
mentioned, the notation used in the following expressions will be the same as in 
ref. (1). The quantity to be calculated is, of course, the square of the inelastic form 


factor defined as: 


where J is the spin of the ground state and 


re =| et exp [iq:r|dt- 


. 2 . . . . . 
The calculation of Fi, is most easily done choosing a system of coordinates with 


the + axis in the direction of q, and using the expansion of exp [ig] in series of 
spherical harmonics. 

We list below the normalized wave functions which we have used, confining 
ourselves to the p nucleons and writing only the angular and spin part; 1 is proton 
and 2 is neutron; the isotopic spin is obviously 0 for all the wave functions 
written down. 


6Li (only LS coupling): 


pera state = SOU V E (pipa + Pi Ps — PYP2) » 
pesca stato = CO VE (2prp2 + Pi Pa + Pi pz) + other terms . 
Here M = M'=—1, 0, 1; are the triplet spin functions; %).= 2 if M'=— +1 and 


Ly = V 5 st N 0; the other terms in ya a8 well as the other values of M’ 
do not intervene in the calculation, with our choice of the 2 axis. 


14N (LS coupling) 
M PERSE [ope e 
cond sie = VO VI (pi pi + Di DI — PIPA) + 
M ’ a en — ro pere PALI 
vi nia VI (2pîp: + pips + Pi pi) + other terms , 


with the same meaning of the symbols as above except that: 1) a bar indicates 
a hole; 2) % = Va if M'= +1 and &,,— aft at M'=T0 
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915. 
MN (JJ coupling): 
” Pitt PRE ME 
Weround state — . 
VE (PE.-4 Pha + Pha DE. -3)> M = 0, 
ya “he VITE VEDI 4 Phas + VE PR ui Piz 3 tn =P}, M= TF1, 
excited state — Ad re È 
VE{— VEDRA Pins + VIDI D+ np), M=0, 


where we have used n and p instead of 1 and 2 to indicate neutron and proton 
and we have symmetrized the wave function with respect to n, p interchange in 
order to have a total T= 0. 

The expressions for Fi, are summarized in the table below; we have also given 


the expressions for F2 (compare for this the formula (129) of the review paper 
by HoFsTADTER (?)). 


TABLE I 
Nucleus | Coupling| Transition Fa Fà, 
; = ni 3 2\2 2 2 \2 2 
LI (0-=-4.18) MeV) — — Le) exp nr pesi exp af 
| 72270 \ 2» 18» 2y 
| 

| leglorg=a q? 

14N LS 0+3.95) MeV; — —|—|e __ 

| Medi 2 a(S) sn ( 1 
i | 5 q? g 
] — — —)exp | —— 
11/0 ge 42 y 2y 

LUNI JJ (0+3.95) MeV| — —|—] exp |—— 

| 72 72\» v 

| 


It appears that the inelastic to elastic ratio at 90° should be one in both cases 


around 230 MeV. The corresponding cross section is small but should nevertheless 
be measurable. 


The authors wish to express their gratitude to Prof. G. MorpurGo for sug- 
gesting this problem and for heipful discussions. 


¢) R. HOFsTADTER: Ann. Rev. Nucl. Sci., 7, 231 (1957). 
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LIBR RICH RECENSIONI 


L. D. LaAnDAU and E. M. LIFSAITZ — 
Statistical Physics, Pergamon 
Press, London, 1958; pagg. 484, 
80 scellini. | 


Si tratta del quinto volume del Corso 
di Fisica Teorica, il grande sforzo di 
questi autori di coprire in un trattato 
tanta parte di fisica. Una edizione in 
inglese era gia stata pubblicata nel 1938, 
ma la presente edizione vuol essere una 
completa rielaborazione piuttosto che un 
rifacimento della precedente. 

Indubbiamente si tratta di un libro 
notevole, che non dovrà mancare à 
nessuna biblioteca di fisica ed a cui gli 
specialisti dovranno continuamente rife- 
rirsi. Ma non si tratta di un libro facile 
e, tanto meno, da consigliare agli stu- 


denti. Il carattere un poco enciclopedico 
ne rende difficile lo studio sistematieo, 


‘anche a causa delle forti discontinuità 


di livello degli argomenti trattati. 

Dopo una parte generale di termo- 
dinamica statistica, che fa subito uso 
di concetti quantistici, seguono nume- 
rose le applicazioni che vanno dalle più 
tradizionali alle più moderne, e che in 
parte sono anche contributi originali 
degli autori. Così la teoria del liquido 
di Fermi e di Bose, le discontinuità nel 
calore specifico ed altri paragrafi ancora, 
sono molto preziosi perchè ci permettono 
una conoscenza diretta e concisa dei 
lavori originali dei due illustri autori e 
della scuola russa contemporanea. 
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